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THE MAGELLANIC CLOUDS. XV. ON THE TILT OF THE LARGE 
CLOUD 


By VirGiIntA McKisppen NAIL AND HARLOW SHAPLEY 
HARVARD COLLEGE OBSERVATORY 
Communicated June 15, 1955 


The preceding communication of this series! (paper XIV) presented evidence, 
derived from the photometric study of 321 cepheid variable stars, on the inclination 
of the bar of the Large Cloud to the plane perpendicular to the line of sight. No 
positive indication of tilt was found, but the method is not sensitive. The scatter 
in apparent median magnitude of cepheids of all ordinary period lengths was found 
to be so large that, whatever the difference may be in the distances to the ends of the 
bar, it is concealed in the magnitude dispersions (see Figure 1 of Paper XIV). The 
positive contribution of Paper XIV lies in the discussion of the variable stars and of 
a certain peculiarity in the frequency of their periods. 

To test the possibility that variable stars might reveal the inclination of the as- 
sumed central plane of the Large Cloud as a whole, we have carried through a 
special study of variable stars in two regions, G and H, which are 4°8 apart on op- 
posite sides of the bar. This investigation extends an earlier discussion of the me- 
dian magnitudes in various sections of the Cloud.? It bears on the indication of a 
conspicuous tilt derived by de Vaucouleurs, who based his determination in part 
on the shape of luminosity contours over the Cloud and in part on some of the var- 
iable stars for which Harvard observers had published periods and apparent median 
magnitudes. The elongated contours were believed to indicate the tilt of a flat- 
tened system, and the variables were used to detect which side of the Cloud was 
nearer the observer.* 

The magnitudes available to de Vaucouleurs were based on various secondary 
sequences. They were not as homogeneous as those we can now provide. For the 
present investigation of tilt we have based the two secondary sequences for regions 
G and H directly on a common master-sequence, using the flat field ADH plates for 
the intercomparisons. 

Fifty-six new periods were derived to supplement earlier work in these regions, 
and twice that many variables were measured for the determination of type. The 
corrections to the median magnitudes already published, while not large, are suf- 
ficient to modify the de Vaucouieurs result, so far as cepheid variables are concerned. 
Two series of Bruce plates, centered, respectively, on regions G and H, were used 
for the measures of periods and magnitudes. 

For the present account of our work most of the details of procedure can be 
omitted. Table 1 and Figure 1 present the summarized data. To supplement the 
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TABLE 1 
Test FoR TILT OF THE LARGE CLOUD 
Mean Median : Mean Harvard 
No. of Period Period De Error Reprint 
"al (Days) (Days (Mag. (Mag.) No. 
6.700 4.299 +0.04 This paper 
7.083 4.189 +0 .05 This paper 
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Fic. 1.—Systematic deviations from a mean period-magnitude curve for eleven regions in the 
Large Magellanic Cloud 
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evidence from regions G and H, we have incorporated in Table | the similar mate- 
rial for the bar from the preceding paper of this series and some earlier results for 
three other regions.” 

For each cepheid for which we have determined the period, P, and the photo- 
graphic light curve, we compute the apparent median magnitude from the period- 
magnitude relation 


mM, = 17.14 — 2.08 log P. 


Then for each region we compute the mean value of the deviations of the observed 
median magnitudes 7) from the foregoing period-magnitude curve, that is, we com- 
pute the mean of mo — m,. When for a given area the mean deviation is positive, 
most of its median magnitudes fall below the period-magnitude curve, which sug- 
gests that, on the average, the variables of that area are more distant; when the 
deviation is negative, the variables are perhaps nearer. But the mean-square 
errors of the means (sixth column of Table 1) are sufficiently large to advise caution 
in assuming that the southwest part of the Cloud is nearer. The difference of two- 
tenths of a magnitude between the mean systematic deviations of region G and 
region H could also represent merely the average difference in space absorption for 
the two regions. 

Ignoring the possible difference in absorption, assuming that the magnitude 
standards are properly homogeneous and that the data are sufficiently numerous 
for this test, and taking the distance of the Cloud as 170,000 light years, we can 
compute that region G appears to be 17,000 light years nearer than region H. If 
we interpret this difference as the result of a tilt of a flattened system, the inclina- 
tion to the tangent plane is about 50°. 

The evidence from regions E and F, however, lends no support to this deduction. 
The difference is better attributed to differential absorption, to uncertainties in the 
magnitudes, and to the intrinsic wide spread of cepheids about the mean period- 
magnitude relation. 

In Figure 1 the data for the bar of the Cloud have been divided into six parts. 
All six sections are more positive than any region outside the bar. Probably this 
result indicates, as mentioned in the preceding paper of this series, that more than 
average space absorption is concentrated among the stars of the bar. It should be 
noted in this connection that region D borders on the great 30 Doradus complex of 
bright and dark nebulosity and that the positive deviation of the region doubtless 
owes something to this location. 

In the course of the study of fields G and H, a considerable number of variables 
of types other than classical cepheids were found. In Table 2 these variables are 
listed (as in Table 1 of Paper XIV), and in Table 3 (as in Table 2 of Paper XIV) new 
variables heretofore unannounced are given, with positions, log periods, median 
magnitudes, and amplitudes for the cepheids and with positions, types, and maxi- 
mum magnitudes for noncepheid variables. 

Figure 2 shows the period-magnitude relation for regions G and H. As found 
along the bar of the Cloud, a considerable scatter is also found here in the median 
magnitudes; but these variables, which are a couple of degrees from the center of 
the bar, do not fall below the standard period-magnitude curve as do those of the bar. 
Apparently these two fields are fairly clear of dust, although both of them contain 
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TABLE 2 


NONCEPHEID VARIABLES 


HV ———-M aAGNITUDES——-—— PERIOD HV -—— MAGNITUDES 
No. Max. Min. 1 Min. 2 (Days) No. Max. Min. 
Eclipsing Irregular 
2659 16.0 16.6 pee ts 2627 15.0 16.5 
2731 15.2 15.9 pe 2665 14.9 15.5 
5936 14.7 15.8 ee 2700 14.3 16.3 
12411 14.8 15.7 : 5532 14.7 15.3 
12413 15.8 16.4 + , 5554 13.5 14.4 
12414 15.7 16.2 5889 15.0 18.5 
12440 14.65 15.75 15.05 1.92149 5891 14.0 14.9 
5914 14.4 15.4 
Long Period 
1001 14.1 16.8 542.1 5933 14.3 15.5 
1004 13.8 16.3 548 3 5949 14.8 15.4 
12439 16.3 18.5 261.8 5961 14.1 14.9 
5964 13.8 14.6 
Short Period 
12425 15.1 15.6 5966 13.4 14.2 
5967 15.8 16.7 
Irregular 
986 14.5 15.5 6003 15.1 15.6 
996 14.1 16.0 12407 14.4 15:7 
2243 13.5 14.5 12420 13.8 14.7 
2255 14.0 16.4 12437 14.7 15.8 
2289 14.0 14.8 12441 13.6 14.6 
2595 13.5 14.6 
TABLE 3 
New VARIABLES IN ReGions G AND H oF THE LARGE CLOUD 
CO-ORDINATES 
HV No. X sd LoG or MEDIAN Amp. 
6 a PeRIOD MAGNITUDE (Maa.) 
12406 1818 3600 0.737 15.70 0.80 
12407 1818 4650 Irregular 14.40 ns 
12408 1902 3552 0.498 15.75 1.10 
12409 2148 3384 0.677 15.30 0.70 
12410 2286 2160 0.760 15.30 0.90 
12411 2400 4656 Eelipsing 14.80 : 
12412 2400 5238 0.582 15.55 1.30 
12413 2820 4734 Keclipsing 15.80 
12414 2832 2532 Eclipsing 15.70 
12415 3012 882 0.663 15.50 1.40 
12416 3198 2478 0.594 15.85 0.70 
12417 3210 1848 0.563 16.05 0.70 
12418 3330 2028 0.481 15.95 0.90 
12419 3426 5328 0.507 16.10 1.00 
12420 3720 5400 Irregular 13.80 
12421 3738 864 0.729 15.60 1.00 
12422 3738 1482 0.633 15.90 1.00 
12423 3882 546 0.528 15.90 1.00 
12424 4068 3990 0.620 15.90 1.00 
12425 4278 2928 Short 15.10 
12426 4350 4674 0.406 15.70 1.00 
12427 4452 2778 0.718 15.50 1.00 
12428 4560 870 0.849 15.50 1.00 
12429 5028 3678 0.595 15.70 1.00 
12430 5196 888 0.670 15.65 1.10 
12431 5202 3636 0.488 16.00 1.00 
12432 5244 4362 0.630 16.00 1.00 
12433 5412 1296 0.874 15.45 0.90 
12434 5634 4410 0.448 15.75 1.30 
12435 13752 18612 bi 16.35 0.90 
12436 13968 18762 0.636 15.85 1.30 
12437 14838 18198 Trregular 14.70 
12438 14880 18426 er 16.00 1.00 
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TABLE 3 (Continued) 
NEw VARIABLES IN ReGions G aNb H of THE LARGE CLouD 


CO-ORDINATES 


HV No. X Y LoG or MEDIAN Amp. 
a ’ PERIOD MAGNITUDE (MaG.) 
12439 * 14886 13508 Long 16.30 
12440 14892 17628 Keclipsing 14.65 
124417 L5815 17956 [Irregular 13.60 bias 
12442 16068 18522 15.90 0.60 
12433 L6806 18060 16.10 0.60 
* HV 12439 = susp. 43 (Harvard Ann., Vol. 90, No. 1, 1933). 
+ HV 12441 = susp. 47 (Harvard Ann., Vol. 90, No. 1, 1933). 


some bright nebulosity. Each area covers about one and a half square degrees, 
centered as follows: 

G: X = 3,810’, y = 2.880": 
H: X = 15,570’, Y = 16,140’. 
The two regions are much alike in the distribution of periods over the interval from 
two to thirty days and, as shown in Table 1, in their mean and median periods. 
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Fia. 2.—Period-magnitude diagram for regions G (©) and H (@) 


They do differ, however, in the amount of recognized bright nebulosity. Region 
H is much more nebulous, and that condition may well be responsible for the fact 
that 18 of its variables (Table 2) are irregular supergiants, whereas only 7 irregular 
variables are found in region G. Three long-period variables are found in region 


H and none in region G. 
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Note on the Small Cloud.—It has been suggested that the Small Magellanic Cloud 
is also a quasi-barred spiral like the Large Cloud and that its elongation is an indica- 
tion of its tilt. The variable stars provide no supporting evidence. For 211 ceph- 
eid variables in five areas outside the central body of the Small Cloud, we have 
assembled the mean systematic deviations of the observed median magnitudes from 
the adopted period-magnitude curve, as shown in the accompanying tabulation. 


No. of Systematic 
Area Position Variables Deviation 
A Southwest 37 0700 + 0"03 
B Northeast - 55 —- 038+ .03 
C Southeast 41 —- O08 + .04 
D East 4] + O12 .04 
E Northwest 37 —0.17 + 0.038 


Summary.—The shape of the Large Cloud of Magellan, whether flattened or 
roughly spherical, remains undetermined. So far the evidence favors a chaotic 
form, but future studies with radio telescopes and possibly with stellar spectro- 
scopes may clarify the matter. Our examination of the magnitudes of several 
hundred variable stars (38 of which were newly discovered in the course of the work) 
gives no positive evidence that these stars lie along a plane; but the criterion is not 
sensitive. A similar examination for the Small Cloud, involving 211 cepheids, is 
also inconclusive. 

A by-product of the work is the evidence that a considerable amount of dusty 
obscuration affects the bar of the Large Cloud. Another by-product is an indication 
that the scatter in median magnitudes along the period-magnitude curves is not 
due to photometric inaccuracies but rather is attributable to other factors, including 
thickness of the Clouds in the line of sight, localized obscurations, and inherent 
dispersions in the median luminosities of classical cepheids. 

' These PRocEEDINGs, 41, 185-190, 1955; Harvard Reprints, No. 407. 

* These PROCEEDINGS, 38, 281-289, 1952; Harvard Reprints, Nos. 359 and 360. 

8 Observatory, 74, 26, 1954. 


RECONSTITUTION OF ACTIVE TOBACCO MOSAIC VIRUS FROM 
ITS INACTIVE PROTEIN AND NUCLEIC ACID COMPONENTS* 


By H. FRAENKEL-CONRAT AND RoBLEY C. WILLIAMS 
VIRUS LABORATORY, UNIVERSITY OF CALIFORNIA, BERKELEY 
Communicated June 17, 1955 


Much recent evidence from chemical, physicochemical, electron microscopical, 
and X-ray studies has resulted in a definite concept of the structure of the tobacco 
mosaic virus (TMV) particle.'~> It appears that about 2,800 protein subunits of 
a molecular weight near 18,000 are arranged in a helical manner to form a rod with 
a hollow core. The nucleic acid is believed to occur as strands in the core. Elec- 
tron micrographs which support this concept have been obtained of the virus at. 
various stages of disaggregation.*-> A protein isolated from infected plants has been 
found to reaggregate—first to short pieces of the presumed helix lying on end and 
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resembling disks with central holes and then to much longer, but inactive, rods of 
the diameter of the virus yet free from nucleic acid. It has now been possible to 
achieve the co-aggregation of inactive virus protein subunits and inactive virus 
nucleic acid to give nucleoprotein rods which appear to be infective. 

Preparation of Protein and Nucleic Acid Components.—TMV was dialyzed against 
pH 10-10.5 glycine buffer (0.01 M) or pH 10.5 carbonate-bicarbonate (0.1 M) at 
3° C. for 48-72 hours. Undegraded virus was separated by cold ultracentrifuga- 
tion, and the supernatant was brought to 0.4 saturation with ammonium sulfate. 
The protein alone precipitates (optical ‘density 2¢9 »,/optical densityogo ,, = R = 
0.65), leaving only nucleic acid (R = 2.0) in the supernatant; if this separation 
is not clean, longer alkali treatment is necessary. The protein moiety is precipi- 
tated twice more with 0.25—0.35 saturated ammonium sulfate, dialyzed, brought to 
pH 7.0-8.0 with NaOH, and finally again freed from heavy particles, such as un- 
degraded virus, by ultracentrifugation. The protein gives a water-clear solution 
at pH 7; the masked —SH group isstill present. The spectrum resembles that of a 
mixture of tryptophan, tyrosine, cysteine, and phenylalanine, simulating the com- 
position of the protein, although the minimum (at 250 my) is not quite as low 
(max./min. = 2.4 versus 2.9) (Fig. 1); P analyses (0.01—0.03 per cent) indicate 
removal of about 95-98 per cent of the nucleic acid. Evidence for the absence of 
detectable virus particles will be discussed below. 

The nucleic acid fractions from such alkali-degraded TMV are not as effective 
for reconstitution as that obtained by the detergent method.? A virus solution 
(1 per cent) containing | per cent sodium dodecy! sulfate is adjusted to pH 8.5 and 
held at 40° for 16-20 hours. Following this treatment, ammonium sulfate is 
added to 0.35 saturation, and the protein precipitate is separated by centrifugation. 
When the supernate is refrigerated, from 60 to 90 per cent of the nucleic acid precip- 
itates (R = 2.0) and is centrifuged off the next day, it is further purified by re- 
peated resolution in ice water and precipitation with two volumes of cold ethanol 
and a few drops of 3 M pH 5acetate. The nucleic acid solution is finally subjected 
to cold ultracentrifugation to remove any traces of virus. The virtual absence of 
protein is indicated by a minimum in O.D. near 230 my (max./min. = 3.0) (Fig. 1). 

In a few preliminary ultracentrifuge experiments, kindly performed by Dr. 
Howard K. Schachman, the nucleic acid preparations exhibited one principal 
boundary with a sedimentation coefficient of about 8 8, resembling the preparation 
of Cohen and Stanley.’ The protein, in pH 9 borate (0.01 M) or in 0.01 M NaCl 
at that pH, exhibited largely a single boundary with a sedimentation coefficient of 
4.5 8. At lower pH values larger components, presumably due to aggregation, 
were observed. 

Reconstitution of Active Virus.—For reconstitution of the virus, 1 ml. of ap- 
proximately 1 per cent protein solution is mixed with 0.1 ml. of a 1 per cent nucleic 
acid solution. Opalescence appears after addition of a suitable buffer; 0.01 ml. 
of pH 6 acetate (3 M) has given the best results, but phosphate (pH 6.3 and 7.0) 
and pH 6 ammonium acetate have also been used successfully. The samples are 
held at 3° for at least 24 hours. They may then be directly diluted and assayed. 
More often they were ultracentrifuged, the pellets redissolved in water, traces of 
insoluble material separated by centrifugation, and aliquots of the opalescent 
supernatants diluted for spectrophotometry. Most of the protein and 40-60 
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Fig. 1.—Ultraviolet absorption spectrum of purified virus pro- 
tein and nucleic acid. 





per cent of the nucleic acid*® were in the pellet, and the spectrum was that expected 
for a nucleoprotein such as TMV (max./min. = about 1.17). On the basis of the 
O.D. of TMV (0.27 at 260 my for an 0.01 per cent solution), approximate concen- 
trations were calculated, and the sample was diluted to a range of 100 to 10 ug/ml. 
for assay. For pellets composed of protein only (absorption maximum at 280 my), 
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calculations were based on the O.D. of the virus proteins (0.13 for an 0.01 per cent 
solution at that wave length). Assays were performed on groups of 8-12 plants 
(Nicotiana glutinosa), distributing each of an equal number of samples (including 
at least one standard and often a solvent blank) over 10 equivalent half-leaves. 
It was indeed quite surprising to find that the reconstituted nucleoprotein prep- 
arations produced local lesions (from 2 to 30 per half-leaf), when tested at 10-100 
ug/ml., which were indistinguishable in appearance from those of TMV at 0.1 
ug/ml. In contrast, no activity was observed when the protein or nucleic acid 
alone was tested up to 1,800 and 52 ug/ml., respectively. When 0.005 per cent 
TMV was held for several hours in 5 per cent protein or 0.5 per cent nucleic acid, it 
was found fully active upon dilution and assay. This strongly suggests that the 
noninfectivity of the two components cannot be attributed to inhibition phenomena 


(Table 1). 


TABLE 1 
Activity OF THREE PREPARATIONS OF RECONSTITUTED VIRUS AND OF TMV art 
DIFFERENT Assay LEVELS; ABSENCE OF INHIBITORS FROM COMPONENTS 


scat eanares -AssaYs——_>_—_- 
MATERIAL TESTED ug./ MI. Lesions/ Half-Leaf* 
Reconstituted virus: 
25,000-rpm. pellet 100 27 +54 31 
120+ 
30 18 - 
25 35 — 
10 9 7 14f 
40,000-rpm. pellet 100 223 — 
TMV stock preparation 0.25 27 
15 
0.1 9 (6.1-12.1)f 
0.02 5 
TMV-protein added at 1000:1 0.1 10 
TMV-RNA added at 100:1 0.1 5 
13 


* Each figure is the average of 10 equivalent half-leaves. The 3 columns in top half of table refer to 


separate preparations. 
+ This is one of the preparations which appeared to increase in activity during storage. 


t Average and range of 9 assays (10 half-leaves each) of one TMV preparation over 3 months. 


When each of the two components was diluted for assay to 0.2 and 0.02 mg./ml., 
respectively, prior to mixing, no activity was generated. Also, when the complete 
reaction mixture (5-10 mg./ml.) was diluted to 0.1 mg./ml. one minute after ad- 
dition of the buffer which triggers the aggregation, no activity was obtained. 
About one hour at room temperature appears to be required for the formation of 
any active rods. These experiments represent convincing control data; they also 
seem to prove that we are dealing with a definite chemical reaction and to exclude 
most other interpretations. The observation that the nucleic acid gradually loses 
most of its activity during several weeks of storage at 3° (in aqueous solution, pH 
5) also serves as a control experiment, indicating the crucial nature of its physical 
state. Treatment of the nucleic acid (600 ug.) with ribonuclease (0.2 ug) in 10-4 
M magnesium sulfate rendered it unable to combine with the protein to produce 
active rods; also, substitution of other nucleic acids (DNA from thymus and RNA 
from turnip yellow mosaic virus") yielded almost nucleic acid-free pellets which 


were inactive (Table 2). 
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A few preliminary experiments on the properties of the reconstituted virus in- 
dicated somewhat greater lability to alkali (pH 9) than that of TMV. Ribonu- 
clease (0.2 ug for 5 mg.) also appeared to decrease the activity. A fractionation 
of the nucleoprotein material by centrifugation at 25,000 rpm. for 30 minutes 
yielded a pellet of considerably higher infectivity than the residual half of the 
material, pelleted at 40,000 rpm. (1 hour). Material so purified seemed to increase 
in activity upon storage in 0.3-0.5 per cent solution at 3°. It also appeared to 
become more resistant to alkali and to ribonuclease. 

Electron Microscopic Studies —Electron microscopy was used qualitatively to 
study the shapes and sizes of the components involved in the formation of the 
virus rods and quantitatively! to assess the degree of homogeneity of the suspen- 
sions as well as to obtain counts and length distributions of the reconstituted rods. 
In the protein, the only visible particles, aggregates presumably of the original 
4.55 material, appeared as disks about 5-15 my thick and with central holes. The 
diameter of the disk was 15 mu and that of the central hole about 4 my (Fig. 2). 


TABLE 2 


E-rFEcT OF REACTION CONDITIONS ON REGENERATION OF VIRUS ACTIVITY 


> — ASSAYS 
REACTION CONDITIONS* ug-/MI. Lesions/ Half-Leaf 

Protein (1 per cent) + RNA (0.1 per cent): 

1 minute, pH 6 100 0.6 

24 hours 100 10.2, 13.1f 

96 hours 100 13.0 
Protein (0.01 per cent) + RNA (0.001 per cent): 

24 hours, pH 6 100 0.6, 0.2t 

96 hours 100 0.1, 0.2T 
Protein + RNase-digested RNA 100 0.0 
Protein + TYMV-RNA 100 0.3 
Protein t 500; 625; 1800 0.0; 0.0; 0.0 
RNA? 30; 52; 230 0.0; 0.0; 0.3 


* Protein and RNA stand for the preparations isolated from TMV by the me »thods described in the 
text. TYMV-RNA is the ribonucleic acid isolated from turnip yellow mosaic virus. RNase is ribo- 
nuclease. 

t Reassayed dilute solution after one week. 


t Different preparations, assayed about 10 times, generally on separate plants (see n. 9); in the 
customary half-leaf assay method, up to 0.5 lesions were occasionally obtained with the same prepara- 
tions. 


The appearance of these perforated disks is identical with that of partially poly- 
merized X-protein.’ In the nucleic acid solution only occasional poorly defined 
fibrils could be detected. The reconstituted rods appeared to be identical in 
shape and size with intact TMV, except for a greater randomness of length (Fig. 3). 

Preparations of the purified protein and nucleic acid were examined for any 
electron microscopic evidence of contamination with intact TMV. In a specific 
instance the protein was sprayed upon the specimen screens at a concentration of 
0.1 per cent and the RNA at 0.007 per cent. These are tenfold and sevenfold 
greater concentrations of the two components than the highest levels used in the 
assays of the reconstituted virus (100 ug./ml.). The protein preparation showed 
no rods of a length even as great as 100 my in six droplet patterns each of which had 
a volume of approximately 3 X 10-*° ml. The RNA material also showed no rods 
in six droplets of approximately this average volume. From these figures it can 
be calculated that the constituent protein and RNA solutions contribute fewer 
than 5 X 10’ typical TMV particles per milliliter to the reconstituted virus in the 
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highest concentration assayed. Additional counts made on the standard TMV 
used as an infectivity control showed that approximately 7 x 10% particles per 
milliliter, of typical length, are required to produce the 10 lesions per leaf custom- 





Fic. 2.—Electron micrograph of the TMV protein used in the reconstitution ex- 
periments. The particles are characteristically disk-shaped, with central holes. 
X 120,000. 





Fic. 3.—Particles of the reconstituted tobacco mosaic virus. Their morphology 
is identical with that of normal TMV, except for a greater proportion of short par- 
ticles. The longest rod in this field is about 300 my long. 60,000. 


arily obtained. Inasmuch as several preparations of the reconstituted virus gave 
about twice that number of lesions at the 100 ug./ml. level, we can conclude that 
the starting materials failed by at least a factor of 30 to contain enough contam- 
inating TMV particles to account for the final infectivity. This conclusion rests 
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upon the supposition that the specific infectivity of any TMV particles that might 
remain as contaminants in the protein and RNA preparations is the same as that 
possessed by the untreated control TMV. 

Counts and length distributions were secured for the particles of both TMV and 
the most active reconstituted virus, the 25,000-rpm. pellet described above, in order 
to ascertain the relative biological efficiency of the two materials in terms of numbers 
of particles per milliliter per lesions per leaf. The results of the counts and length 
distributions are given in Table 3. From this table it appears that about one- 
tenth of the total particles in the reconstituted virus, representing about one-third 
of the total mass, were of length ca. 300 muy, the length of the monomer of normal 
TMV. These particles were therefore only about 3 per cent as infective as the 
particles of similar length in the control TMV. It thus appears that polymeriza- 
tion of the protein by nucleic acid to form ~300-mu rods is fairly frequent but that 
only a fraction of the rods is reconstituted with the structural faithfulness necessary 
for infectivity. 


TABLE 3 
ELECTRON Microscope Counts or PARTICLES OF TMV Anpb REcoNSTITUTED VirRUS AT ASSAY 
CONCENTRATION 
Assay 


Concentration Lesions Total Particles/ M1. of Particles/ M1. 
Preparation (ug./ M1.) Half-Leaf Particles/ Ml. Length 290-310 My 10 Lesions 


Control TMV 0.1 10 1.6 X 16 7.0 X 108 7 x 106 
Reconstituted virus 
(25,000 rpm. ) 10 10 2.2 & 102 2.0 K 10” 2 X< 10” 


The length distribution of the rods of this pelleted material exhibited a reasonably 
random character for lengths less than about 260 my, but a highly uniform length 
between 290 and 310 my for the 10 per cent of the rods falling in this range. Only 
3 per cent of the particles were of lengths greater than 310 mu. This is in contrast 
to aggregates of X-protein, which show a complete randomness of distribution of 
lengths. 

To ascertain whether the nucleic acid was localized in the center of the newly 
formed virus rods, Dr. R. Hart applied his technique of detergent treatment, followed 
by electron microscopical analysis.> He found the reconstituted rods appreciably 
more labile to SDS than was standard TMV, but after SDS treatment for 10 
seconds many rods were partially degraded and showed a central strand of material 
protruding from the ends, as does standard TMV after 60 seconds of reaction 
(Fig. 4). These strands disappeared when ribonuclease was added. 

Summary.—The preparation from TMV of protein and RNA fractions which 
tend to recombine at about pH 6 to form a nucleoprotein carrying virus activity 
(0.1—-1 per cent of that of TMV) is described. 

An electron microscopic search revealed no TMV rods in either of the two com- 
ponent solutions at a concentration level thirty fold to three hundred fold greater 
than those at which the reconstituted virus was assayed. In the latter, on the 
other hand, up to about one-third of the material consisted of rods of the typical di- 


ameter and length of TMV, many, if not all, containing a nucleic acid core. 

The concentration, time, and pH dependence of virus regeneration is that of a 
typical chemical reaction. Freshly prepared RNA is required for appreciable 
reaction; degraded RNA, or nucleic acids from other sources, are inactive. 
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No inhibition was observed if know small amounts of TMV were added to con- 
centrated solutions of each of the components and subsequently diluted and 
assayed. 





Fic. 4.—Particles of reconstituted TMV treated briefly with hot detergent. 
Rods are seen with strands of nucleic acid projecting from their ends. This appear- 
ance is identical with that found for normal TMV more severely treated. 55,000. 


The evidence thus seems reasonably complete that, under the conditions de- 
scribed, TMV nucleic acid enters into combination with TMV protein subunits 
and favors aggregation to rods, some of which are of sufficient length and structural 
integration to carry infectivity. 


The capable assistance of Mrs. B. Singer and Mr. J. Toby is gratefully ac- 
knowledged. 
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ENZYME FORMATION IN PROTOPLASTS OF BACILLUS 
MEGATERIUM* 


By O. E. LANDMAN AND S. SPIEGELMAN 
DEPARTMENT OF BACTERIOLOGY, UNIVERSITY OF ILLINOIS, URBANA, ILLINOIS 


Communicated by M. M. Rhoades, July 22, 1955 


INTRODUCTION 


Yesults of the last two years * have led rather compellingly to the conclusion 
that the induced synthesis of enzyme occurs via a mechanism which condenses 
amino acids on a template. The possibility that amino acids are “‘activated”’ prior 
to adsorption on the template seems likely from recent preliminary reports.* 4 
In any event, no evidence for a stable intermediate, smaller than the fully formed 
enzyme, has yet emerged. If one accepts the template hypothesis, it is evident 
that further understanding of the mechanism of enzyme formation will not be 
achieved through attempts at analysis of intermediate stages of protein synthesis, 
since none exist. This view forces concentration on experiments designed to eluci- 
date the nature and functioning of the template. Data obtained? * with intact 
cells synthesizing enzyme suggest that the template is RNA. However, the 
distance between the observations and the derived conclusions in such experiments 
is too great for certainty. The need is clear for a subcellular fraction possessing 
enzyme-forming ability and amenable to more directly interpretable dissection of 
its components and their function. The recent work of Gale and Folkes®* 7 on 
ruptured cell preparations of Staphylococcus aureus, and of Zamecnik and Keller® on 
a liver microsome fraction, indicates that the attainment of such systems is close 
at hand. 

Weibull’s® observation on ‘“‘protoplast’’ formation resulting from lysozyme treat- 
ment of Bacillus megaterium cells in hypertonic medium suggested that an- 
other analyzable subcellular system was available. Protoplasts are morpho- 
logically quite dissimilar from the rod-shaped cells from which they are derived, 
being spherical and smaller. Each rod usually yields between 2 and 3 spheres. 
Lester” and Beljanski!! were able to show that lysozyme treated B. megatertum 
could incorporate labeled amino acids. The question remained whether protoplasts 
could be induced to synthesize enzyme. This potentiality has been realized. It 
has been found possible to devise a properly supplemented stabilizing medium, 
which permits protoplasts of B. megaterium to form beta-galactosidase at rates 
comparable to those of intact cells. Simultaneously, Wiame and his collaborators 
(personal communication) succeeded in obtaining formation of arabinokinase in 
lysozyme treated preparations of Bacillus subtilis. 
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It is the major purpose of the present paper to describe the methods of obtaining 
such preparations from B. megaterium and their enzyme-synthesizing properties. 


MATERIALS AND METHODS 


a) Strain and Conditions of Growth. The B. megaterium strain used is a lactose- 
positive mutant derived from Weibull’s® strain KM by selective enrichment in a 
medium containing lactose as the primary carbon source. It has a tendency to 
revert, particularly in dense liquid cultures, to a lactose-negative rough variant, 
insensitive to lysozyme. Stability is achieved by keeping stock slants on minimal 
(1 per cent) lactose medium. The minimal medium has the following composition: 
KNO;, 1 X 10-!' M; K.SO,, 1 X 10-? M; Na,HPO;, 1 X 10-' M; MgSO, 5 xX 
10-* M; MnSO,, 1 XK 10-4 M; Fe.(SO,);,5 K 10-°M; ZnSO, 5 X 10° M; and 
aspartic acid, histidine, valine, cysteine, all 50 ug./ml. 

It was empirically found that the growth stage of the cells employed is crucial 
to reproducible success in obtaining enzyme formation in protoplasts. The follow- 
ing procedure has consistently yielded active preparations: A suitable volume of 
2 per cent Difco peptone is inoculated at a level of about 1/200 that of a stationary 
phase culture, using cells from stock slants or young 2 per cent peptone cultures. 
The culture is incubated at 30° C. with shaking to provide vigorous aeration for 
12-15 hours. By this time it will have reached the stationary phase and have an 
optical density (O.D.¢60) of 0.55—0.64, using filter 66 of the Klett-Summerson color- 
imeter and the medium as a blank. The culture is then “rejuvenated” by dilution 
with 2 per cent peptone to an O.D.¢69 of between 0.12 and 0.16 (2.3 XK 107 cells/ml.) 
and reincubated until an O.D.¢¢0 of 0.20—-0.22 is reached. It is essential that growth 
beyond this extent be avoided. The cells are then harvested and washed with 
0.2 M Naes/HPO, adjusted to pH 7.8. They are then suspended (adjusted to an 
O.D.¢60 of 0.70) in the media described below and converted to protoplasts with the 
aid of lysozyme (Armour and Company) at a level of 200 ug./ml. The density of 
protoplasts in the usual induction experiment is such that 1 ml. contains 1 mg. of 
protein. 

b) Enzyme Assay.—The enzyme being followed is a beta-galactosidase, which 
hydrolyzes the synthetic chromogenic substrate ortho-nitrophenyl-8-p-galactoside 
(ONPG). Its activity can therefore be readily ascertained by procedures similar 
to those employed" for the assay of Escherichia coli lactase. 

A study of the purified B. megater‘um enzyme and its induction has been made 
and will be reported separately. The properties of the enzyme are such as to lead 
to the following assay mixture for the attainment of maximal activity; ONPG, 
1 X 10-*M; NasHPO,, 10~! M, adjusted to pH 7.84; MnCh, 1 XK 10-*M; glucose, 
0.8 M. The enzyme reaction is allowed to proceed at 40° C., the color due to the 
ortho-nitrophenol released being followed at suitable intervals by reading in a Klett- 
Summerson colorimeter, using filter 42a. Activities are reported as millimicro- 
moles of substrate hydrolyzed per minute. Assays are always carried out on 
lysed preparations, since neither protoplasts nor cells exhibit their full enzyme con- 
tent when intact. Dilution of enzyme leads to marked inactivation, which, how- 
ever, can be avoided by including 0.1 M methionine in the diluent. 

c) Analytical Methods.—The procedures of Ogur and Rosen'* were employed to 
separate the nucleic acids and protein. DNA was estimated by the Dische™ 
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reaction, and the RNA spectrophotometrically. Protein was analyzed by the Folin 
reagent. 

d) Enzymes.—The following enzymes have been employed in the present study: 
lipase (Nutritional Biochemical Company); 2X crystallized trypsin, crystalline 
ribonuclease, crystalline desoxyribonuclease (all from Worthington Biochemical 
Sales Company); and crystalline lysozyme (Armour and Company). 


EXPERIMENTAL RESULTS 


a) The Stabilizing Medium for Protoplasts.—A medium suitable for the purposes 
in mind must fulfill two principal requirements. Protoplasts suspended in it should 
remain physically intact and physiologically functional for a period sufficient to 
perform an induction. Further, the medium must permit the synthesis of the 
enzyme being followed. 

The retention of beta-galactosidase was employed as an indication of physical 
intactness in the search for a utilizable medium. Fully induced cells, obtained by 
growth in a 2 per cent lactose-peptone medium, were converted into protoplasts 
under various conditions, and the distribution of enzyme between the protoplasts 
and the supernatants was determined at various intervals. Protoplasts were spun 
down by centrifugation in a Serval SS1 run at top speed for two minutes. This 
relatively high speed does not injure the protoplasts. 

The validity of released beta-galactosidase as an indicator of the extent of physi- 
cal destruction was justified in separate correlative experiments. These established 
that the proportionate loss in the enzyme from the sedimented protoplasts was equal 
to that observed for total protein and DNA. Further, direct microscopic counts 
of the remaining protoplasts correlated well with what would have been predicted 
from the measurements of enzyme and the other parameters noted. 

Using the enzyme as the test indicator, it was found that 0.2 M sucrose and 7.5 
per cent polyethylene glycol® failed to stabilize our protoplasts. All the enzyme 
activity was found in the supernatant, even after short incubation periods. On 
the other hand, the following compounds, at the concentrations indicated, were 
effective stabilizing agents: sucrose above 0.6 M, KCl and NaCl at and above 
0.8 M, and NasHPO, at and above 0.5 M. 

The effects of these various agents on enzyme-forming ability of intact cells were 
examined. On the basis of these experiments sucrose was eliminated as a possi- 
bility by virtue of its strong inhibitory effect. The data obtained suggested that 
the sodium phosphate buffer was the most desirable. 

b) Induction Medium.—An examination was next undertaken to discover sup- 
plements necessary to achieve rates of enzyme formation in protoplasts comparable 
to those of intact cells held under the same conditions. As a result of these studies, 
the following induction medium was eventually evolved: 0.5 M NasHPO,, adjusted 
to pH 7.84; enzymatic hydrolyzate of casein, 2 per cent; MnCl, 1 X 10-4 M; 
adenosine triphosphate (ATP), 1 mg./ml.; hexose diphosphate (HDP), 6 mg./ml.; 
lactose, 6 X 107? M. 

Of the components listed, the hydrolyzate’ HDP and lactose are mandatory. 
The omission of any one of them results in the complete absence of enzyme syn- 
thesis. The reasons for the amino acids and lactose requirements are evident, 
since the latter is the inducer and the former are needed for the formation of 
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new enzyme molecules. The HDP is apparently necessary for the energy-generat- 
ing system. 

In actual experiments with protoplasts the lactose is omitted during the period 
of protoplast formation. This avoids including in the experiment enzyme formation 
by cells or stages intermediate between intact cells and protoplasts. At 
30° C. and a level of 200 ug/ml. of lysozyme, virtually complete conversion to 
protoplasts is attained in 30 minutes. The inducer and other agents are 
introduced at the end of this period. The progress of the lysozyme treatment 
can be easily followed by direct microscopic observation and count in an hemocytom- 
eter chamber. To avoid foaming during this and all subsequent incubations, 
aeration is continuously maintained by means of a device which rotates the tubes 
in a plane inclined to the horizontal. 

The 0.5 M concentration level of the buffer is 
necessary for prolonged experiments with the pro- 
toplasts. It was found possible to lower the con- — 
centration to 0.3 M and still obtain sufficient sta- 
bilization to permit enzyme synthesis. Figure | 
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compares enzyme formation in cells and proto- 
plasts in the induction medium adjusted to 0.5 M 
and 0.3 M of the NasHPO,g. Except for length of 
the lag period, enzyme synthesis in intact cell is 
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essentially the same at both concentrations. Pro- ti 

toplasts, on the other hand, exhibit a difference re- poe: 

flected by the cessation of enzyme formation in 300 

the 0.3 M buffer after 3.7 hours of induction. At 200 

this point it can be shown that 50 per cent or 100 

more of the enzyme formed has been released into 05 

the medium, indicating the onset of physical dis- HOURS 

solution of the protoplasts. In the 0.5 M buffer thease be oclinanliieneineen: 
enzyme synthesis continues for considerably pended in 0.5 M and 0.3 M induc- 


longer periods tion medium. Activities are ex- 
; f ‘ 5 pressed in terms of muM of sub- 
In dealing with protoplasts, there are a few useful strate hydrolyzed per minute 


facts which should be mentioned. A continual en- per milligram of protein. 

ergy source is absolutely essential for the mainte- 

ance of enzyme-forming ability. Manipulations in the cold, absence of HDP, and of 
aeration are therefore to be avoided. Complete disappearance of the ability to 
synthesize enzyme is followed by even brief exposures to a medium containing amino 
acids but no HDP. Even the omission of the amino acids leads to severe irrever- 
sible loss in enzyme-forming ability. 

Finally, if it is desirable to wash and resuspend protoplasts subsequent to some 
treatment (e.g., with ribonuclease), resuspension should be in supernatant obtained 
from a freshly prepared protoplast suspension if maximal enzyme synthesis is to 
be observed. The nature of the responsible factor(s) here is now under investi- 


gation. ' 

c) Cells versus Protoplasts.—Qualitatively, the metabolic properties of intact 
cells versus protoplasts are remarkably similar, provided that the comparison is 
made under identical conditions. Thus the need for relatively elevated levels 
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of inducer and external supply of amino acids is not unique to protoplasts but is 
rather referable to the hypertonic medium. Surprisingly, this holds true even for 
the HDP requirement, as is shown in Table 1. It is evident that in 0.5 M buffer 
enzyme synthesis in intact cells is impossible if the HDP is omitted, a condition 
which disappears at the lower buffer levels. The utilization of HDP evidently 
cannot be a means of distinguishing intact cells of B. megateritum from protoplasts 
under hypertonic conditions. This is a situation markedly dissimilar from that 
reported® for S. aureus, in which ruptured cells can use HDP and intact cells cannot. 
It is when one examines responses to various enzymes that striking differences 
between protoplasts and cells begin to emerge. This is clearly exhibited in Table 
2 in the case of lipase and trypsin. Intact cells are completely insensitive to the 
enzymes whereas the synthetic ability of protoplasts is completely abolished. 
These results illustrate a point worthy of the attention of those concerned with 
performing and interpreting experiments with subcellular fractions. Based 
simply on the observation recorded with lipase, one might perhaps be led to con- 
TABLE 1 
Errect oF HyPpERTONIC MEDIUM ON HEXOSE-DIPHOSPHATE REQUIREMENT FOR ENZYME 
SYNTHESIS IN INTACT CELLS 
(Enzyme Activities Expressed as muM of ONPG Hydrolyzed per Milliliter per Minute) 


Enzyme Enzyme 

Synthesized Synthesized 

Buffer HDP after 90 Min. Buffer HDP after 90 Min. 

Concentration (Mg./M1.) Incubation Concentration (Mg./ML1.) Incubation 
0.5 M 0 0 0.05 M 0 966 
0.5 M 6 294 0.05 M 6 1,080 


TABLE 2 
IxrFFECT OF TRYPSIN AND LIPASE ON ENZYME FORMATION IN CELLS AND 
PROTOPLASTS * 


Enzyme Present Cells Protoplasts 
None 1,090 620 
Trypsin (100 ug. /ml. ) 1,140 0 
Lipase (100 ywg./ml.) 1,020 0 


* Cells or protoplasts were suspended in inducer-free induction medium (0.5 M) and 

incubated with the indicated enzyme for 1 hour, subsequent to which inducer (0.06 M 

lactose) was added. The enzyme formed in the next 2 hours is recorded in terms of 

muM ONPG hydrolyzed per milliliter. 
jecture that a lipid is a key component of the enzyme-forming mechanism. How- 
ever, the fact is that the loss of enzyme-synthesizing ability is a simple consequence 
of physical dissolution of the protoplast. After incubation with either lipase or 
trypsin at the levels indicated, few protoplasts can be recovered. 

It is thus important in any given case to demonstrate that an observed inhibition 
of enzyme synthesis is not the result of a generalized destruction. This caution is 
also relevant to experiments involving ribonuclease (RNase) and desoxyribonu- 
clease (DNsae). Lysis of protoplasts by RNase has been observed under certain 
conditions both in our own laboratory and by Brenner (personal communication). 
To be interpretable, experiments must be accompanied by evidence that the en- 
zyme treatment has resulted in a selective removal of the homologous compound. 

An extensive examination has been made of the effects of RNase and DNase 
on enzyme synthesis in intact cells and protoplasts. The details will be reported 
in extenso in a subsequent publication. We may here summarize the principal 
findings. RNase at levels of 1 or 2 mg./ml. does not suppress enzyme formation 





aetna 


ee 





Hy 


0.5 M buffer, and resuspended in 0.5 M induction 
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in intact cells suspended in either hypertonic (0.5 M) or normal (0.05 /) induction 
media. Indeed, frequent stimulations are observed. Protoplasts suspended in 
0.5 M buffer induction media gave variable responses to RNase. However, con- 
sistent inhibitions are observed if the exposure to RNase is carried out in 0.32 M 
induction medium. A typical experiment is shown in Figure 2. Here protoplasts 
were incubated aerobically in 0.32 M inducer-free induction mixture with and with- 
out | mg./ml. of RNase for 1 hour. They were then spun down, washed with 

















mixture containing inducer, as is described in the [ 

last paragraph of section b of “Experimental Re- ae i 

sults.”” Samples removed after resuspension 

showed no destruction of protoplasts by the cri- 400}- 

teria noted above. It is evident from Figure 2 3 

that exposure to the RNase leads to drastic in- f eoil: 

hibition of enzyme-forming ability. ‘The apparent 5 

absence of the lag period seen in Figure 2 is due Z 

to the fact that the ‘‘zero”’ time sample is taken "1 

after the RNase pretreatment. The protoplasts 

have thus had 1.5 hours to adjust to the hyper- 1ooF 

tonic medium. Analyses of protoplasts subjected 

to incubation with RNase show that 80-90 per , F 

cent of the RNA can be removed by the enzyme, , ‘ yours * . 

with no significant loss of DNA or protein. Fic. 2.—The effect of pretreat- 
DNase did not inhibit enzyme formation in ae "cae of eaten 

intact cells. On the contrary, in analogy with to form beta-galactosidase. Ac- 

the observations of Lester’ and Beljanski,"! strong weal of tak stiade sede ae 

stimulations were frequently obtained. Two re- minute per milligram of protein. 


sults were encountered on treating protoplasts with 

DNase (200 y/ml.). One was the selective removal of DNA in amounts ranging 
between 20-50 per cent of the original content. The other involved a concomitant 
loss of about 30 per cent of the RNA, with little or no loss of protein in either case. 
Inhibitions of enzyme synthesis were observed only when both nucleic acids were 
extracted and in these cases the inhibitions were consistent and severe. 

The results support the obvious and frequently suggested conclusion that RNA, 
derived from DNA, is a key component of the enzyme-forming mechanism. A 
more detailed discussion of the significance of these observations is deferred until 
it can be done with the aid of reconstitution experiments using the system described, 
and made possible by the results reported here. 


SUMMARY 


A procedure and stabilizing media have been described which permit the study 
of enzyme synthesis in protoplasts of B. megaterium. Enzyme-formation properties 
have been compared with those of intact cells under the same conditions. The 
protoplasts differ mainly in being susceptible to enzymatic resolution. 

Treatment with either RNase or DNase can result in the selective removal of the 
homologous substrates. When 50 per cent or more of the RNA is thus removed in 
non-induced protoplasts, loss of enzyme-synthesizing capacity results. 
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THE DISTRIBUTION OF PARENTAL PHOSPHORUS ATOMS AMONG 
BACTERIOPHAGE PROGEN Y* 


By Gunther S. Stent and Niels K. Jernet 


VIRUS LABORATORY, UNIVERSITY OF CALIFORNIA, BERKELEY, AND CALIFORNIA INSTITUTE OF 
TECHNOLOGY, PASADENA 


Communicated by W. M. Stanley, July 15, 1955. 


Approximately half the deoxyribonucleic acid (DNA) contained in a population 
of T2, T4, or T6 bacteriophage particles reappears among the descendants ultimately 
issuing from phage-infected bacterial host cells. This transfer is not due to the 
reincarnation of entire, intact parental DNA units in progeny guise, since at least 
half the DNA of each of those descendant particles which harbor the transferred 
atoms must be of nonparental origin.* For an understanding of the mechanisms 
involved in the reproduction of the hereditary structures of the bacteriophage, it is 
desirable to know the distribution of the parental atoms among the progeny pop- 
ulation, i.e., the extent to which the atomic identity of the parental DNA has been 
conserved or destroyed. It is the purpose of this communication to present the 
results of experiments which indicate that most of the transferred phosphorus atoms 
of the parental DNA are distributed over at least 8 but no more than 25 of the prog- 
eny. A more detailed description of these results will be presented elsewhere. 

The basis of these experiments is that the bacteriophages lose their infectivity 
upon decay of radiophosphorus P* incorporated in their DNA, the rate of inactiva- 


2 


tion being proportional to the number of P*? atoms per particle.* The fraction of 
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P* disintegrations which inactivate the T2 or T4 particles in which they occur is 
0.10 at 4° C., this “efficiency of killing” having been established for “nonparental”’ 
radiophosphorus atoms, i.e., for those assimilated into the phage DNA from the 
phosphorylated constituents of host cell or growth medium.* * Since neither the 
transferred phosphorus atoms nor those whose decay leads to inactivation appear 
to reside in any “special fraction” of the bacteriophage DNA,*°* it would seem 
reasonable that the decay of transferred P*? atoms should similarly inactivate the 
progeny particles harboring them. We have adopted this at present unprovable 
assumption and have endeavored to detect the presence of parental P*? atoms in the 
descendant phages by observing the lethal effects of the decay of these atoms on 
the progeny population. 

Survival of the Infectivity of Progeny of a Highly Radioactive Parent.—Most of the 
phosphorus atoms transferred by P**-labeled T2 or T4 parents appear among the 
first fifty or so progeny particles to mature after the termination of the eclipse 
period; the last hundred or so progeny particles to mature just before lysis of the 
infected cell receive very little or none of the parental atoms.* ® 7 An experiment 
was carried out to examine whether the parental DNA has been dispersed so effee- 
tively that its atoms are widely distributed among the individual early progeny 


particles. 

A culture of strain B of Escherichia coli was infected at a multiplicity of 0.2 
phage per cell with a 2-hour-old stock of T2 labeled with P*® at 2,100 mc/mg 
(or 2,800 atoms of P** per particle). The infected bacteria were freed from un- 
adsorbed phages by centrifugation and incubated in broth at 37° C. for 14, 16, 
21, and 30 minutes, when potassium cyanide and chloroform were added to aliquots 
to induce premature lysis.» * The lysates, which corresponded to yields of 0.11, 
2.2, 35, and 230 mature intracellular T2 particles per infected bacterium, were stored 
at 4° C. and assayed from day to day for their infective titer. 

The result of this experiment is presented in Figure 1, where the logarithm of the 
titer of the surviving infective phages has been plotted against the fraction of P** 
atoms decayed by the time of assay. It is seen, first of all, that the highly radio- 
active parental T2 stock loses its infectivity rapidly in the exponential manner 
characteristic of inactivation by P* decay,* 4 so that less than 1 per cent of the 
initial titer remains after 8 hours. The progeny populations do not, however, 
appear to be subject to any significant loss of titer, none of them being inactivated at 
a rate superior to 0.2 per cent of that of the rate of inactivation of their parent. It 
appears, therefore, that the majority of the individuals among even the earliest progeny 
do not contain more than 0.2 per cent of the parental phosphorus. Since the average 
transfer to the first intracellular progeny has been estimated to b= 2 per cent parental 
phosphorus per particle,'® it would follow that most of the transferred phosphorus 
resides in a minority of the progeny population. The first progeny, therefore, are 
highly heterogeneous with respect to amount of parental phosphorus possessed by 
each, 

Survival of the Transferability of the Transferred Phosphorus.—A second ex- 
periment was designed to probe further into the nature of the heterogeneous dis- 
tribution of transferred parental atoms among progeny bacteriophages. This 
experiment is based on the notion that in the case of a population of P*-labeled 
bacteriophages of which a minority of the particles contain most of the P**, it should 
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be possible to determine the rate of inactivation of the radioactive minority through 
P*? decay by measuring at various times the fraction of P*? which the entire pop- 
ulation can transfer in single infection to its progeny. Initially, when all the radio- 
active phages are still infective, the global P** transfer to the progeny should be 
somewhere near 50 per cent; at later stages, when most of the radioactive phages 
have been inactivated by decay of some of the P*? atoms they contain, but when all 
the nonradioactive phages are still alive, the transfer of P*? to the progeny should 
be nearly zero. The rate at which the transferability of the P** atoms contained by 
the heterogeneous population decreases with the fraction of P**? atoms decayed 
should, then, in principle, be a reflection of the specific P*? content of the radio- 
active minority. 
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Fig. 1.—The infective titer as a function of the fraction of P*? atoms decayed of a parent stock 





of T2 labeled with P*? at 2,100 mce/mg (left-hand panel) and of its progeny collected at stages of 
intracellular growth corresponding to various burst sizes (right-hand panel). The stippled lines 
indicate survival curves corresponding to rates of inactivation constituting various percentages 
of the rate of inactivation of the radioactive parent. 

A control experiment was carried out to test the practicability of this idea. 
Strain B of EL. coli was infected at a multiplicity of 0.1 phage per cell with a mixed 
population of phages consisting of 97 per cent nonradioactive T4 particles and 
3 per cent purified T4 particles labeled with P*? at 57 me/mg (or 75 atoms of 
P*? per particle). The infected bacteria were separated from any unadsorbed 
radioactivity or phages by repeated centrifuging and washing in salt-poor broth, 
and the radioactivity adsorbed to the bacterial cells was counted. The infected 
bacteria were incubated in salt-poor broth for 27 minutes at 37° C., at which time 
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potassium cyanide was added to induce lysis,® yielding an average burst of 100-180 
phages per infected cell. Bacterial debris and unlysed bacteria were then removed 
by centrifugation, and a second culture of bacteria was infected with an aliquot of 
the supernatant fluid, containing the progeny phages. After allowing adsorption 
of the progeny, the infected bacteria were centrifuged and the radioactivity con- 
tained in this last bacterial pellet counted. The ratio of radioactive counts in this 
last pellet to those initially adsorbed to the first bacterial suspension represents the 
fraction of the P* transferred from parent to progeny. The mixed population of 
parental phages was stored at 4° C., and this transfer experiment was repeated on 
several successive days, while an ever decreasing fraction of the radioactive minority 
population of T4 still remained alive. 

The results of this experiment are presented in Figure 2, a, where the logarithms 
of the percentage of the surviving particles of the radioactive T4 population and of 
the percentage of P*? transferred to the progeny by the mixed population of radio- 
active and nonradioactive parental phages are plotted against the fraction of the 
P*®? atoms which have decayed by the day each transfer experiment was carried 
out. It is apparent that the transferability of the P** does, in fact, decrease con- 
tinuously from its initial value near 40 per cent to a value below 10 per cent at a 
rate similar to that at which the infectivity of the radioactive T4 particles them- 
selves is lost. The two inactivation curves do not remain parallel after the transfer 
value has dropped below 15 per cent, because, due to certain experimental com- 
plications, minimum transfer values between 2 and 5 per cent are usually en- 
countered even when multiplication or maturation of the adsorbed radioactive 
phage has been prevented.?. This minimum may increase further when the radio- 
active parents have been inactivated by agents which prejudice the injection of 
the parental DNA into the host cell.'' The control demonstrates, however, that 
the rate of loss of transferability of the P® contained in a phage population is an 
adequate reflection of the rate of inactivation, and hence of the specific P** content, of the 
particles that possess the radiophosphorus. 

It has been shown that the first-generation progeny possessing P*? transferred to 
them from a randomly labeled parental population again transfer about 30-50 per 
cent of these P*? atoms to their own, second-generation progeny. Hence it would 
follow that the rate of loss of the transferability of the parental radiophosphorus 
possessed by the first-generation progeny should reveal the amount of the parental 
P*? carried by individual progeny particles. 

A suspension of E. coli B was infected at a multiplicity of 0.1 phage per cell 
with a purified stock of T4 labeled with P** at 200 me/mg (or 260 atoms of P* 
per particle). The infected bacteria were permitted to lyse to yield a burst of 100 
first generation radioactive progeny per cell. Analysis of these first-generation 
progeny showed that they contained 43 per cent of the parental P**. The lysate 
was stored at 4° C., and a transfer experiment similar to the control experiment just 
described was carried out on successive days in order to determine the amount of 
P*. transferable from the first-generation to the second-generation progeny. 

The results of this experiment are presented in Figure 2, b, where the logarithms 
of the percentage of the surviving particles of the radioactive parental T4 popu- 
lation and of the percentage of the P*? transferred from the first-generation to the 
second-generation progeny are plotted against the fraction of P*? atoms decayed. 
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It is seen that the rate of loss of transferability amounts to no more than 6 and no 

less than 2 per cent of the rate of death of the parental population. It may be 

concluded, therefore, that those particles of the first-generation progeny which f 
contain most of the transferred P*? atoms appear to possess between 2 and 6 per 

cent of the phosphorus of their progenitor. 
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Fic. 2.—a: The percentage of the initial infectivity and of the 
transferable radiophosphorus of a stock of T4 labeled with P*%? at 
57 mc/mg as a function of the fraction of P*? atoms decayed. 06: 
The percentage of the initial infectivity of a parent stock of T4 
labeled with P*? at 200 mc/mg and the percentage of the parental 
phosphorus contained in the first-generation progeny transferable 
to the second-generation progeny as a function of the fraction of 
the P*? atoms decayed. The stippled lines indicate survial curves 
corresponding to rates of inactivation constituting various percent- 
ages of the rate of inactivation of the radioactive parent. 





Discussion.—The survival of the infectivity of progeny of a highly radioactive 
parent, presented in Figure 1, indicates an upper limit to the dispersion of the pa- 
rental phosphorus atoms, for the accuracy of titration of the infective centers in 
that experiment would appear to be such that the initial inactivation of no more 
than a third of any of the progeny populations should have escaped our notice. 
The bulk of the transferred phosphorus can, therefore, reside in no more than a 
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third of the first 50 or so progeny, that is to say, in no more than about 15-20 
particles. These particles would then carry at least 3 times the average of 2 per cent, 
or 6 per cent parental phosphorus each. The survival of the transferability of the 
transferred phosphorus, presented in Figure 2, 6, indicates a similar dispersion, 
for the 50 per cent of the parental atoms which are transferred to the descend- 
ants appear in progeny particles each of which seems to contain between 2 and 6 
per cent of the parental phosphorus. From this it would follow that the number 
of phages which carry the bulk of the transferred phosphorus atoms must be be- 
tween 50/2 and 50/6, 1.e., between 25 and 8 particles. The fact that the parental 
atoms are distributed over at least 8 particles indicates that the parental bacteri- 
ophage DNA experiences a certain dissociation in the course of its reduplication 
within the host cell. It should be noted that this dissociation proceeds in single 
infection, i.e., under conditions in which it is certain that the progeny are, in fact, 
the true offspring of the individuals whose radioactive atoms they bear. Previous 
experiments, by which evidence for a dispersion of the parental DNA had been 
adduced, were invariably carried out in multiple infection, where it was possible, 
and sometimes certain, that some of the transferred atoms had not actually been 
derived from that parent which was responsible for the issue.2 '' '* The quanti- 
tative inferences drawn from the present experiments, however, depend for their 
part on the as yet unverified assumption that the “efficiency of killing” by decay 
of transferred P*? atoms is also 0.10. The significance of the present findings for 
various proposals concerning the mechanism of DNA replication will be consid- 
ered elsewhere. 

Summary.—The descendants of highly P**-labeled T2 or T4 bacteriophages are 
inactivated by radiophosphorus decay in a manner which indicates that the bulk 
of the atoms transferred from the deoxyribonucleic acid of the infecting parental 
virus particle to that of its progeny is distributed over no less than 8 and no more 
than 25 progeny particles. 

* This investigation was supported by grants to the Virus Laboratory from the National Cancer 
Institute of the National Institutes of Health, United States Public Health Service, and from the 
Rockefeller Foundation; and to the California Institute of Technology from the National Founda- 
tion for Infantile Paralysis. 
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STUDIES OF MOLECULAR INTERACTION IN ISOIONIC: PROTEIN 
SOLUTIONS BY LIGHT-SCATTERING* 
By Serce N. TimasHerr, Howarp M. Dinrzis,t Joun G. KrrKwoop, 
AND BERNARD D. COLEMAN 
STERLING CHEMISTRY LABORATORY, YALE UNIVERSITY, NEW HAVEN, CONNECTICUT 
Communicated August 16, 1955 
The technique of light-scattering has been extensively employed in recent years 
in physicochemical investigations of proteins in solution." ? It is useful not only 
for the determination of molecular weights of macromolecules but also for the 
measurement of their thermodynamic interaction. The excess turbidity 7, arising 
from composition fluctuations, in a solution of two components is related to the 
chemical potential us and molecular weight / of the solute in the following manner: 
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where C2 is the solute concentration in grams per milliliter, uo“ is the excess chemi- 
cal potential of the solute, and H is the optical constant. From the intercept of 
HC:,/r, plotted as a function of C2, the solute molecular weight WM, may be deter- 
mined, and, from the slope of the curve, the derivative (Oye /OC2)7,, measuring 
the departure of the solution from ideal behavior may be calculated. Attractive 
forces between the solute molecules lead to negative values of this derivative, 
whereas repulsive forces lead to positive values. 

Kirkwood and Shumaker‘ have shown that fluctuations in charge and charge 
configuration of protein molecules can lead to a supplementary attractive force 
between them, in excess of the usual] electrostatic force associated with their average 
charge distributions. Although it is not, in general, possible to distinguish be- 
tween these two types of force by thermodynamic measurements, that part of the 
fluctuation force arising from total charge fluctuations of the molecules of a salt- 
free isoionic solution may be isolated from effects due to all other intermolecular 
forces, since it gives rise to a term in the excess chemical potential proportional to 
the square root of the solute concentration, while all other intermolecular forces, 
both van der Waals forces and electrostatic forces associated with permanent and 
fluctuating multipoles, contribute only terms proportional to the first and higher 
powers of the concentration. Thus the theory of Kirkwood and Shumaker leads 
to the following power series in C,'”* for the function HC./r in salt-free isoionic pro- 
tein solution, for which the average molecular charge is zero: 
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where (2°),, is the mean-square charge of a protein molecule in protonic units 
e, N is Avogadro’s number, k is Boltzmann’s constant, 7’ is the thermodynamic tem- 
perature, and D is the dielectric constant of the solvent. The virial coefficient B 
reflects the combined effect of all types of intermolecular forces.> According to 
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equation (2), the function HC,/r7 should approach linearity in C.'’? at high dilu- 
tions, and the charge fluctuation (z*),, may be determined from the limiting slope 
of HC,/r, plotted as a function of C,'”’. 

It is the purpose of the present investigation to provide an experimental test of 
the predictions of the theory of Kirkwood and Shumaker for salt-free isoionie pro- 
tein solutions and to employ equation (2) for the determination of the total charge 
fluctuation of a representative protein. The protein used was Armour crystalline 
bovine plasma albumin (BSA). An 8-10 per cent solution of this protein was made 
up in redistilled water and was deionized by passing it over an ion-exchange col- 
umn, according to the method of Dintzis and Oncley.6 The conductivity of the 
resulting solution corresponded exactly to its pH. This material was then clari- 
fied for light-scattering by centrifuging it for three hours in a Spinco Model L cen- 
trifuge at 40,000 rpm and then passing it, under one pound of nitrogen pressure, 
through a sintered glass filter of ultrafine porosity, especially designed for light- 
scattering work.’ The water used as diluting solvent was purified by redistilling 
it in an all-Pyrex still. . 

The light-scattering measurements were carried out in a Brice-Phoenix pho- 
tometer, with the proper modifications for use with small cells. The measurements 
were done in 3-ml. square glass cells held in a specially designed cell-holder. In each 
set of measurements the turbidity of the solvent was first determined, and then a 
series of small increments of concentrated stock protein solution were added from 
a Gilmont ultramicroburet. Each series of measurements was repeated several 
times. Points obtained in overlapping dilution series indicated that this procedure 
was free of errors resulting from protein denaturation or from the introduction of 
small amounts of dust during addition and mixing. 

The value for the refractive increment of the protein used throughout this study 
was 0.1953. The results obtained with one preparation of isoionic Armour BSA in 
salt-free aqueous solution are shown by the filled circles in Figure 1. In this 
figure the light-scattering data are plotted in the usual manner as a function of 
concentration. The data, which cover a concentration range of 0.005—1.8 per cent, 
reveal a marked curvature in the dilute region. This curvature is most pronounced 
in the region below a concentration of 0.2 per cent, which, in this case, is deter- 
mined by a set of 40 points. 

When plotted against the square root of the concentration, as is shown in Figure 
2, the points fall very nearly on a straight line over the four-hundred-fold concen- 
tration range studied in these experiments. The curve is well represented by 
equation (2), with the following coefficients determined by the method of least 


squares: 


HC, 


T 


= 1.27 X 10> — (8.35 X& 10-7)C." — (2.46 * 10-8) C., (3) 


the standard deviation being +0.024 X 10-°. 

The intercept of this curve corresponds to a molecular weight of 78,900 without 
the use of the customary correction for depolarization.’ This value is in good agree- 
ment with values previously reported for BSA from light-scattering measurements.’ 
An ultracentrifugal analysis of our sample of albumin revealed the presence of a few 
per cent of a heavier component, which should be sufficient to account for the fact 
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Fic. 1.—Light-scattering of Armour bovine serum albumin; usual plot. Filled circles: isoionic 
salt-free aqueous solution; open circles: solution in 0.001 M NaCl. 
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that our value of the average molecular weight was somewhat higher than that ob- 
tained by other methods. Results obtained in a control run done in 0.001 M 
NaCl, which covered the same concentration range as the salt-free measurements, 
are given by the open circles in Figure 1. In this case the data are linear in con- 
centration over the entire concentration range studied, as would be expected for a 
protein solution in the presence of salt. Extrapolation to zero concentration re- 
sults in a value of 83,300 for the apparent molecular weight. The difference in 
intercepts between the experiments on salt-free solution and those in 0.001 M 
NaCl could be accounted for largely in terms of the difference of refractive incre- 
ments of BSA in the presence and in the absence of salt,’ since in our measurements 
the same value of dn /dcz was used in the two cases. A small contribution would also 
be expected from the thermodynamic interaction constant between the protein and 
the salt. 

The possibilities that the data for the salt-free isoionic solutions were influenced 
by instrumental artifacts or by loss of protein due to surface denaturation or ad- 
sorption were eliminated by control experiments which will be described in a subse- 
quent paper. The fact that the measurements in 0.001 M NaCl yielded values of 
HC,./r, linear in the protein concentration, which did not exhibit a sharp rise at high 
dilutions is good evidence for the absence of artifacts. The possibility that the up- 
ward curvature in the salt-free solutions was due to progressive ionization of the 
protein as its self-buffering power diminished with decreasing concentration was 
eliminated by control experiments in 10-° HCl. Results identical with those 
in water alone were obtained. 

The theoretical possibility that our results should be attributed to a concentra- 
tion-dependent aggregation arising from specific short-range attractive forces may 
be rejected, since such an effect would always be linear in the protein concentration 
at high dilutions. In fact, the small value of the virial coefficient B of equations 
(2) and (3) shows that there is rather a close balance between short-range attractive 
forces and repulsive forces, the latter making a positive co-volume contribution to 
B which nearly cancels the negative contribution of the former. 

In a recent paper Dandliker' has also reported on light-scattering measurements 
on isoionie solutions of BSA. From his data it is not possible to decide whether 
HC./r approaches proportionality to the square root or to the first power of pro- 
tein concentration at infinite dilution. However, Dandliker’s measurements cover 
a smaller concentration range than our own, and the region of high dilution is cov- 
ered by fewer points, with a larger random error, than in our experiments. 

Our experiments provide convincing evidence of the validity of the charge- 
fluctuation mechanism of molecular interaction in protein solutions proposed by 
Kirkwood and Shumaker. In the absence of a means of estimating the relative 
effect. of other short-range forces, quantitative interpretation is possible only for 
the long-range component which gives rise to the term proportional to the square 
root of concentration. From the coefficient of this term in equation (3) we are 
able to calculate a value of 3.5 protonic units for (z%),,’*, the root-mean-square 
charge fluctuation of BSA in isoionic solutions. This value is in excellent agree- 
ment with the value 3.4 which has been calculated from titration data obtained for 
this protein by Tanford.!! 

In conclusion, we remark that further studies on Armour BSA, human serum 
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mercaptalbumin, and bovine serum mercaptalbumin carried out in this laboratory 
are in complete agreement with the results presented here and may be quantita- 
tively interpreted by the charge-fluctuation concept. These studies will be reported 
in a separate communication. 


* Contribution No. 1316 from the Sterling Chemistry Laboratory, Yale University. This 
investigation was supported in part by the Office of Naval Research, Contract No. Nonr 659(00). 
It was presented in part at the 125th Meeting of the American Chemical Society, Kansas City, 
March, 1954, and at the 126th Meeting, New York, September, 1954. 
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A NOTE ON THE ORIGIN OF LIFE 
By Martrynas YCAS 
QM RESEARCH AND DEVELOPMENT CENTER, NATICK, MASSACHUSETTS 
Communicated by H.C. Urey, July 26, 1955 


Oparin'! made a significant contribution toward the solution of the problem of the 
origin of life by suggesting that the prebiological environment contained a great 
abundance of organic compounds, whose existence was possible because of the re- 
ducing nature of the atmosphere at that time. The gap between living organisms 
and simple inorganic compounds is thus considerably reduced. Since Oparin made 
his suggestion, Urey* has brought forward geochemical and astronomical arguments 
in favor of the thesis that the original atmosphere of the earth was of a reducing 
nature. Recently it has been shown experimentally* that passage of electrical 
discharges through such an atmosphere produces organic compounds in great 
variety and abundance. The original speculations of Oparin thus appear to have 
been well founded. 

There is, nevertheless, a large gap between the simplest living organism and an 
aqueous solution of organic compounds. I propose here that this gap can be further 
narrowed. 

A biochemical system can be regarded as a system of catalysts regulating the 
transformation of other compounds so as to make available to the system energy 
and matter for its further increase and maintenance. Since the catalysts cannot be 
produced without metabolism, and metabolism cannot proceed without the cata- 
lysts, the origin of such a system has naturally been difficult to comprehend. 
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If one accepts the arguments mentioned above as to the prebiological environ- 
ment of the earth, it will be noted that the transformation of matter would have 
formed an “ecological cycle.”” Under the influence of the energy of light or elec- 
trical discharges, simple compounds (methane, ammonia, etc.) of the original 
atmosphere form a great variety of organic compounds in solution in the ocean. 
These compounds are metastable and tend, at rates determined by the activation 
energies of the processes involved, to disintegrate into the compounds from which 
they originally arose. While in solution in the ocean, the organic compounds will 
interact, forming numerous molecular species. On a two-dimensional figure the 
molecular transformations can be represented as a reticulum, within which cycles 
can be arbitrarily delineated (heavy lines in Fig. 1). 

As has been shown, particularly by the long researches of von Langenbeck,‘ 
primary valence-bond catalysis is a common phenomenon in organic chemistry. 
Hence it is to be expected that many of the compounds in solution will catalyse 
some reaction represented by the reticulum. Now consider the following case: a 
product of cycle A catalyzes a rate-limiting reaction in cycle B, and conversely a 
product of cycle B catalyzes a rate- 
limiting reaction in cycle A. In 
such a case the quantity of matter 
passing through these two cycles 5 ie on a 
will increase at the expense of other 














reaction cycles. A kind of “natural etnies 
selection” based on reaction rates 
will result. Eventually (and auto- ?? ?? 
matically), a system of interlocking aa | 
cycles will be selected which oper- a 
ates at the maximal rate. 

Such a system of self-generating — 
catalysts meets the minimal defini- 
tion of life. In this initial stage, 














however, there are no discrete organisms, and there exists only one living 
thing, the metabolizing ocean. The further evolution of this system presum- 
ably led to the production of catalysts of a high molecular weight and peptide 
nature. These may have agglomerated, and eventually the system would have 
become delineated into small discrete masses, the individual organisms. A better 
understanding of the biosynthesis of proteins and the nature of enzymatic catalysis 
may eventually make it possible to speculate in more detail on the nature of the 
above process. 

This proposal has the advantage that it is not necessary to invoke a series of 
highly improbable events to account for the origin of life. Furthermore, since the 
metabolism of such a primitive system is primarily determined by the catalytic 
properties of compounds rather than by their initial concentrations, metabolic 
systems of a similar kind might arise in any large solution of organic compounds, 
relatively independently of the exact initial composition of the system, given a 
suitably coupled energy source. If bodies similar to the earth occur elsewhere in 
the universe, the occurrence of life with very similar properties may therefore be 
cosmologically a common event. 
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1 A. I. Oparin, The Origin of Life (New York: Dover Publications, 1953). 
2H. C. Urey, these PRocEEDINGs, 38, 351, 1952. 
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ON THE STRUCTURE OF THE STRESS-STRAIN RELATIONS 
By T. Y. THomas 
GRADUATE INSTITUTE FOR MATHEMATICS AND MECHANICS, BLOOMINGTON, INDIANA 
Communicated J aly 22, 1955 


1. Jntroduction.—In a recent publication Truesdell! has emphasized a difficulty 
associated with the usual application of a system of differential stress-strain relations 
of the form 


AT ap 


dt = Fis(mu, ey) eR eee, B (1) 


in which the left member is the total time derivative of the components 7, of the 
stress tensor and the quantities F,, are the components of a tensor invariant of the 
stress tensor and the rate of strain tensor whose components have been denoted by 
nas The difficulty in question arises when one attempts to apply relations (1) toa 
body in rigid motion, and it stems, mathematically speaking, from the fact that the 
right members of (1) are of tensor character under the group of point transforma- 
tions relating moving rectangular co-ordinate systems, while the left members are 
not of this character relative to such transformations. In meeting this difficulty 
Truesdell devised a procedure which led to the addition of other terms to the left 
members of (1), and this provided a set of stress-strain relations of correct tensorial 
character. Later, Noll' obtained stress-strain relations similar to those obtained 
by Truesdell on the basis of the requirement of invariance under the above group of 
transformations. 

While this work must be regarded as constituting an outstanding contribution to 
the theory of the stress-strain relations, we feel that further clarification of the 
difficulty associated with the form (1) of these relations is desirable, and we have 
therefore been prompted to conduct our own investigation of this problem.? We 
have shown in this paper that the difficulty in question can be removed merely by a 
suitable interpretation of relations (1), and hence these relations are, in fact, when 
suitably interpreted, the stress-strain relations for a dynamically sound theory of 
the deformations of a continuous medium. 

2. Interpretation of the Noninvariant Stress-Strain Relations.—The local motion 
of the medium in the neighborhood of a point P can be resolved into an instanta- 
neous rigid rotation about an axis through P, a translation, and a motion resulting 
in distortion or actual change of volume. The components of the local motion in- 
volving rotation and translation will not produce structural changes in the medium, 
i.e., changes inherently connected with variations in the stress and strain rates. It 
appears natural, therefore, to view the stress-strain relations from the standpoint 
of an observer in a moving co-ordinate system relative to which the local rotation 
vanishes. Whether or not one chooses this system so that the local translation 
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likewise vanishes appears immaterial, since the stress tensor and the other quanti- 
ties used in the construction of the stress-strain relations will not be affected by 
translations. 

To develop the idea involved in the above paragraph precisely, let us denote by 
x* the co-ordinates of an arbitrary rectangular co-ordinate system and by y* the 
co-ordinates of a rectangular system which may be moving relative to the x system. 
Then between the x* and the y* we will have relations of the form 


co k 5 
y* = Agp(t)xX + b(t); AakLam = Rm; (2) 


in which the coefficients a,, and 6, are functions of the time ¢ alone. Now select 
an arbitrary point P of the medium, and choose the y system so that relative to it 
the instantaneous rotation of the medium vanishes at the point P. When the y co- 
ordinate system is chosen in this way, it will be called a kinematically preferred sys- 
tem associated with the point P. It will be assumed that this kinematically preferred 
system is selected so that its origin lies at the point P. 

We now make the following assumption. Relations (1) are valid at a point P of the 
medium only when they are referred to a kinematically preferred system associated with 
the point P and evaluated at the origin of this system, 1.e., at the point P. The 
problem of expressing relations (1) in the x co-ordinate system must now be con- 
sidered. When this problem is solved, it will be seen that our interpretation of 
equations (1) as equations valid only at the origin of a kinematically preferred sys- 
tem leads to a set of stress-strain relations which is invariant under the group of 
co-ordinate transformations defined by (2). 

3. Auxiliary Formulas.—The components of particle velocity in the z and y sys- 
tems related by (2) are given by », = dx*/dt and w, = dy*/dt, respectively. By 
differentiation of (2) with respect to the time ¢, we obtain the relations between 
these components, namely, 


es ‘ 
Wa = AgkVk + TD apt + Ruy (3) 
where the dot denotes the ordinary time derivative. We now differentiate (3) with 
. B . . . . 
respect to the co-ordinate y® and then eliminate the derivatives dx"/dy* by means 
» 7 k B . n > 4 out '* ° 

of the relations 02° /Oy" = ag,, which can be derived from (2). This gives 

Wap = GaklpmVk,m + Gapdgp- (4) 


Now, when we interchange the indices a,8 in (4) and then add corresponding mem- 
bers of the resulting equations and the equations (4), we find that the terms involv- 
ing the @’s cancel, so that we obtain 


Nap = QakIpmEekmy (9) 

where 
Nap = , 2(Wa B + Wea); €km = 2Vz m + Um k)- 
It follows from (5) that the quantities e€,,, or are the components of a tensor 
| km Nap 

under the group of transformations defined by (2). Similarly, when we inter- 
change the indices a,8 in (4) and subtract, the resulting equations can be given the 
form 


Ger = Waptpe — PikIais (6) 
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where 
Vas = '/2(Wap — Wea); biz = ‘Jolie — Ve) 
ap / 2\%a,B B,a/} tk /2\%RR k,i/- 
Hence quantities of the type ¢; do not constitute the components of a tensor under 
A tk 
the group (2). 

Denoting by o,3 the components of the stress tensor relative to the x system and 
transforming these quantities into the components 7, of this tensor with respect 
to the y system, we have 

Tap = File A ge. (7) 
Let us now take the total derivative of (7) with respect to the time, to obtain 


drag Mes doiz 


it Agidge + TixLaiApe + Fi~Aaid pp. 
tf 


Eliminating the quantities d,; from these relations by the substitution (6), we 
vasily find that the resulting equations can be written in the form 


Tap a FikVaiA ge, (8) 
corresponding to (7), where 
a drag 
Tap = é a VayT By Aj VByT ar (9) 
dt 
? doix 
a= it — DimTkm — PkmTim- (10) 
¢ 


The quantities 7,, and ¢ given by (9) and (10) are therefore related by the tensor 
transformation law under the group of transformations defined by (2). 

4. Invariant Form of the Stress-Strain Relations.—Let us now consider the co- 
ordinates y* in section 3 to be the co-ordinates of a kinematically preferred sys- 
tem. Then the condition that the rotation vanish is given by pag = 0 at the 
origin of the y system. Hence the left members of (1) can be replaced by the cor- 
responding quantities 7, on account of (9). But, since the quantities 7,, are of 
tensor character by (8) and the right members of (1) are of tensor character by hy- 
pothesis under the transformations of the group defined by (2), it follows that the 
relations 


6a = Fys(€u,- oe er eee (1 1) 


are valid in the a system. These relations represent the stress-strain relations (1) 
relative to the x system and are, moreover, invariant under the transformation 
group (2). 

5. The Absolute Time Derivative-—The reduction of the quantities 7,, defined 
by (9) to the total time derivatives dr,,/dt at the origin of a kinematically pre- 
ferred system is exactly analogous to the reduction of the absolute derivative of a 
quantity along a curve in a Riemann or general affinely connected space at the origin 
of a system of normal co-ordinates in this space. The situation here appears also to 
be connected, although perhaps not quite so closely with the reduction of the com- 
ponents of the covariant derivative to partial derivatives at the origin of normal 
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co-ordinates. Adopting the notation commonly used in the tensor analysis, let us 


now write? 
Dox, _ (44) (12) 
Dt dt] 2 


where the quantities in the right members of these equations are evaluated at the 
origin of a kinematically preferred system. Jn connection with equations (12), we 
assume that the kinematically preferred system is chosen so that its axes have the direc- 
tions of the corresponding axes of the x co-ordinates system at the time t under considera- 
tion. ‘The quantities Do;,/Dt defined in the x system by (12) will now be the com- 
ponents of a tensor under the group (2) which we will call the absolute time derivative 
of the tensor «. It follows immediately that 


Dox . doi, 
Dt dt 


= PimTkm — DiemTim 


since the left and right members of these relations are (a) of tensor character under 
the group (2) and (6) equal at the origin of the y system. Hence the stress-strain 
relations (11) can be written 


Doug si 
Dt 


Fs(en, «++» €33) Oily 5. 933). (13) 


On the basis of relationship (12), it appears that the absolute time derivative 
gives the natural modification of the total time derivative of a tensor when one 
wishes to obtain relations of invariant character under the group (2). In fact, by 
the mere substitution of the absolute time derivative of the stress tensor for the 
total time derivative of this tensor in (1), we obtain a set of invariant and dynami- 
cally correct stress-strain relations. 

6. Final Remark.—If we make the substitutions 


’ “a 2 ’ ) 1 ’ ’ —_ 
Vag = 2(Ve,8 + UB .a) + 2Vq.B tes V8.) = €a8 + Pap 


for the derivatives v, 3 in the left members of the stress-strain or constitutive equa- 
tions derived by Truesdell,‘ we find that we obtain the left members of our relations 
(11) or (13), plus additional terms. However, these additional terms are of the 
nature of terms that can appear in the right members of the relations. It follows, 
therefore, that our relations (11) or (13) and the relations of Truesdell are equiva- 
lent. On the other hand, relations (11) and (13) appear to be strictly similar to 
those obtained by Noll, apart from the fact that the latter are expressed in a form 
involving both covariant and contravariant indices. 

‘©, Truesdell, “The Simplest Rate Theory of Pure Elasticity,” Communs. Pure and Appl. 
Math., 8, 123-132, 1955. Also, W. Noll, “On the Continuity of the Solid and Fluid States,” 
J. Rat. Mech. and Anal., 4, 3-81, 1955. References to earlier work are given in these two papers. 

2 Prepared for the Applied Mathematics Branch, Naval Research Laboratory, Washington, 
D.C. This paper was written while the author was visiting the University of California, Numeri- 
cal Analysis Research, Los Angeles, California. 

3 For the analogous definition of tensors and invariants by means of normal co-ordinates see, 
for example, T. Y. Thomas, The Differential Invariants of Generalized Spaces (London: Cam- 
bridge University Press, 1934), pp. 96-107. The introduction here of relation (12) was discussed 
at the Naval Research Laboratory with Dr. Ericksen, who pointed out to me in this connection 
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that he had previously employed co-ordinate systems in which the rotation vanished, as an aid 
in the construction of invariant relations. See R.S. Rivlin and J. L. Ericksen, ‘‘Stress-Deforma- 
tion Relations for Isotropic Materials,” J. Rat. Mech. and Anal., 4, 323-425, 1955. 

4 Op. cit., p. 124. 

© Op. cit., p. 27. 


COMBINED ELASTIC AND PRANDTL-REUSS STRESS-STRAIN 
RELATIONS 


By T. Y. THomas 


GRADUATE INSTITUTE FOR MATHEMATICS AND MECHANICS, INDIANA UNIVERSITY, 
BLOOMINGTON, INDIANA 


Communicated July 29, 1955 


1. /ntroduction.—The concept of stress-stain relations valid for both the elastic 
and the plastic domain has certain theoretical advantages which appear to com- 
pensate for the fact that such relations are necessarily somewhat more complicated 
than the corresponding relations for either of these domains treated separately.! 
For example, there is the possibility of a more gradual transition from the elastic to 
the plastic state corresponding to the actual behavior of materials; also, the use of 
a single system of equations eliminates the troublesome matter of matching the 
elastic and the plastic solutions, which arises in the process of unloading when dis- 
tinct systems of equations are assumed to govern the elastic and plastic behavior of a 
body.” 

The discussion in this paper is based on the assumption of a set of stress-strain 


relations of the form 


Dogs 


aon Fap(€i, - - +) €335 Ou, » + +» O33), (1) 


in which the F’s are the components of a tensor invariant of the stress tensor o and 
the rate-of-strain tensor e; the quantities Do,,/Dt are the components of the abso- 
lute time derivative of the stress tensor.’ Briefly stated, we have attempted to con- 
struct the tensor invariant F so as to obtain the utmost simplicity consistent with the 
usual assumptions, such as isotropy of the material, homogeneity, and the quadratic- 
yield condition of von Mises—considered here as the extreme case of an inequality, 
since relations (1) are intended to embrace both the elastic and the plastic domain 
as commonly conceived. While this set of stress-strain relations is of differential 
character it reduces approximately on integration, as shown in section 9, to the 
usual algebraic relations of classical elasticity theory for small displacements about 
the natural or unstressed state of the body. In the extreme case of yield, the stress- 
strain relations of this theory become identical with the condition of incompressibil- 
ity and the Prandtl-Reuss relations (sec. 8) which have been taken as basic in the 
recent book by Prager and Hodge.‘ 

2. Definitions.—Rectangular co-ordinates will be employed exclusively, so that 
there will be no distinction between the covariant and contravariant indices of ten- 
sorial quantities, and covariant differentiation, denoted as usual by a comma, will be 
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identical with partial differentiation relative to these co-ordinates. The use of such 
systems involves, of course, the implicit assumption that a rectangular system can be 
found for which the equations of the theory will be valid and will lead to results 
which approximate the results of experiment. An essential part of the solution of 
any problem is therefore the selection of the proper rectangular co-ordinate system. 
The velocity is defined by its components v, = dx“*/dt, i.e., the components of 
velocity are the total time derivatives of the co-ordinates x* of the material particles. 
We also denote by o,3 the components of the symmetric stress tensor. Other quan- 
tities which will enter the discussion can now conveniently be defined as follows: 


. ot * , 
Stress deviation: ogg = Gag — '/s0104s, 
Rate-of-strain tensor: €.3 = '/2(¥a.g + Usa) 

7 . . . * 
Rate-of-strain deviation: €gg = €ag — '/s€:%as; 


where it is understood that repeated indices, unless the contrary is stated, are to be 
summed over the values 1, 2, 3, and that 6,, denotes the usual Kronecker delta. 

3. Dynamical Relations—The following dynamical equations must be satisfied 
in a continuous medium which is not subjected to body forces, namely, 


dp 


7 + aa = 0 (equation of continuity), (2) 
( 


die 


1 (equations of motion), (3) 
al 


Tup,8 — P 
where p denotes the density. Since the first of these relations is a consequence of 
the assumption of the conservation of mass, it follows, by putting p = 1, that we ob- 
tain M%q@ = 0 or €,4 = 0 as the condition for incompressible flow. 

Equations (2) and (3) above are invariant under the group of orthogonal co- 
ordinate transformations which relate the preferred co-ordinate systems of our Eu- 
clidean metric space. These equations are also invariant under the group of uni- 
form translations, in conformity with the principle that it is not possible from dynam- 
ical considerations to distinguish between two reference frames which are moving 
relative to each other with uniform velocity. 

4. Simplification of the General Stress-Strain Relations.—As the first step in the 
process of simplifying relations (1), we assume that these relations have the form® 


Dou 


Dt = AS €iz + Bis Ciky (4) 


where A and B are tensors which can be taken without loss of generality to be sym- 
metric in their contravariant and also in their covariant indices. One could assume 
the coefficients A‘, and B‘%, in these relations to be functions strictly of the co- 
ordinates «* and the time ¢, in which case the right-hand members of the relations 
will be linear and homogeneous in the quantities ¢, and o. But certain of the 
following requirements will be satisfied with greater facility by sacrificing this 
linearity and allowing these coefficients to be expressed as tensor invariants of ¢ and 
a. We shall find, however, that this can be done in such a way that only simple 
quadratic expressions in the quantities €,, and o.3 are involved in the components 


of the tensors A and B. 
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The above assumption (sec. 1) that the right-hand members of relations (1) are 
the components of a tensor invariant of the stress and rate-of-strain tesors implies 
that the material is isotropic. Now, as generally understood, the requirement of 
isotropy means that at any point x and at any time? the properties of the material 
are independent of direction; this is expressed mathematically by the condition 
that the stress-strain relations retain their strict form under proper orthogonal co- 
ordinate transformations. It follows that A and B must be tsotropic tensors. As is 
well known, the components of the symmetric tensors A and B must therefore have 
the form 


ik ik k is kk 
A ye = hd b a8 + u(6,'5¢ + 5354 i 
ik k . k is ke 
Bye = £6" 548 oo ¢ (6,'68 + 5354 ), 
where A, uw, £& and ¢ are scalars. Substituting these expressions in equation (4), we 
now have 


Dogs 


Di = (Néibag + 2ueag) + (Ecidag + 26048). (5) 


We shall wish in particular to consider the single equation obtained from equation 
(5) by summing on the indices a and 8. We shall likewise need to deal with the set 
of equations which result when this latter equation is combined with equation (5) in 
the manner indicated in section 2 for the formation of a deviation tensor. We thus 
obtain the following relations: 


Doi; ‘ 
= (SA + Quen + (3& + 2W)ou, (6) 
Dt 
Dox ) * ae De * (7 
= 2 af - ap? | 
Dt to eine | 


which are readily seen to be equivalent to the system (5). 

5. The Fundamental Inequality.—It will be assumed that there exists a (positive) 
material constant K, regarded here as one of the quantities characterizing the be- 
havior of the material, such that for all flows we have 


* * r 
TapFap S K. (8) 


A point x of the material at which the equality holds in relation (8) will be called a 
yield point at the time ¢ in question. A region of the material over which the equal- 
ity holds will be said to be in a state of yield. The strict equality in condition (8) 
above is the well-known yield condtion® of von Mises. 
For a region in a state of yield it follows from equation (7) and the definition of the 
absolute time derivative’ that 
. dors ; Doers Pe PS ee hae 
Tag i = Cag Di = Queaplas + 2Caplasg = Y. 
Hence we can write 
* * 
Ser <_< (for 28% ap = K). (9) 
im K 
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Now ¢/u is a scalar, and it can therefore be considered quite generally to be a func- 
tion of the co-ordinates x“ and the time ¢. However, from the standpoint of the 
structure of the theory it is advantageous to represent this scalar by writing 


- 


* * * * 
ic GE ap7 ap + bo agF ap + C, (10) 
m 
where the quantities a, b, c, are scalar functions of position and time. When the 
yield condition is satisfied, the right-hand members of equations (9) and (10) must 
be equal. Assuming this equality to be of the nature of an identity, it follows that 


a = —1/K and bK +c = 0. Hence equation (10) takes the form 
¢ €agd ap a aa? = 
= — —y({l- — }. (11) 
m kK ( K ) 


in which y now appears as a scalar function of the co-ordinates x«* and the time /. 
When yield occurs, this relation reduces to relation (9) above. 
We can now use equation (11) to eliminate the quantity ¢ from the right-hand 


members of equation (7). Hence 


Doig J * ik T ik ( casi) on 5 
— ae 7f{1 — e 12 
Dt b ) €ag | K + Y K Tag { ( ) 


6. Condition of Incompressibility.—It is generally considered that one can as- 
sume with sufficient accuracy that volume changes do not occur during plastic flow, 
i.e., when the material is in a state of yield (sec. 5). This condition of incompressi- 
bility can be realized automatically by choosing the coefficients in equation (6) to be 
suitable scalar invariants.’ 

Let us write equation (6) in the form 


Dei: 
A = €; + Boi, (13) 
Dt 
where 
a ] ; a 3& + 2¢ 
— BA + 2p’ Sede ns. 


We now assume that A — 0 and B > 0 as o2goag > K at any point x inde- 
pendently of the time ¢. It follows immediately from this assumption and from 
equation (13) that the condition of incompressibility ¢;; = 0 is satisfied for plastic 
flow.® 

Corresponding to the procedure in section 5, let us represent the scalars A and B, 
as can be done without loss of generality, by expressions of the form 


* * * * 
A = a+ boggeas} B= c¢ + doggeag, 


in which a, b, c, and d are functions of the co-ordinates x“ and the time ¢. Then, for 
plastic flow, we have a = —bK and c = —dK by the above assumption. Hence A 
and B become 


A = —bK + bongo ag; B seas —dk + dong ap- 
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When we substitute these expressions in equation (13), this equation can be given 
the form® 


Dei: hei: 
_ 2 es + koi, ( 14) 
Dt l — (¢p%ae/K) 
where h and /: are scalar functions of the 2“ and ¢t. Equation (14) reduces to the 
condition of incompressibility ¢;; = 0 for plastic flow. 


7. Homogeneous Material.—The material will be said to be homogeneous if the 
stress-strain relations are left strictly unchanged by co-ordinate translations; this 
condition is sometimes expressed by saying that the properties of the material will 
appear the same to observers at different points. Assuming that this invariance re- 
quirement applies independently of the time i, it follows that the scalar functions y, 
y, h, and k in the stress-strain relations (12) and (14) reduce to material constants for 
homogeneous material. 

The determination of the original scalars \, £, and ¢ in terms of the above material 

* 


° r . . * 
constants uw, y, A, and k, the yield constant K, and the invariants €,3,0,, and 


. * " » . > . Re 
To3% a3 can be obtained from equation (11) and the following two relations: 


ar 
h 


7. * > , 
— . (ongtog K) 


3\ + Qu = 3¢ + 2 =k. 
8. Stress-Strain Relations for Plastic Flow.—In the case of plastic flow the stress- 
strain relations (12) and (14) become 


+ 
Dos 9..." 2u "es ud 
= ERO jb OC ass 15) 
tk? tk? aps 
K 


€x = 0; Seay 2ueag — 
where u and K are material constants. There are only four independent relations 
in the second set of equations (15), (a) owing to the fact that these equations are 
identically satisfied when we sum on the indices @ and 6 and (6) because the com- 
bination o.gD 008 Dt vanishes identically in view of the yield condition. Hence 
equations (2), (3), and (15), when combined with the yield condition, constitute 
a set of 1 + 3+ 5+ 1 = 10 equations for the determination of the 10 dependent 
variables p, vg and ogg. 

9. Small Stresses —Let us now consider the situation for ordinary elastic dis- 
placements of isotropic and homogeneous material where the stresses are small in 
comparison with the stresses which produce yield. It can evidently be assumed here 
that terms in equations (12) and (14) containing the squares of components e. 
are negligibly small in comparison with the value of the yield contant K. Assum- 
ing also that the terms which contain y and k in equations (12) and (14) can likewise 
be neglected and that the absolute time derivative of the stress tensor can be ap- 
proximated by the total time derivative for the usual elastic problems involving 
small stress, it follows that the stress-strain relations reduce to 

dong a tig doi: : me 
dt ~MEabi in (3A + 2y) ex, 
where yu is a material constant and the quantity \ can be taken to have the constant 
value (hk — 2u)/3. These relations can readily be shown to be equivalent to the 


relations 
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doag - 


A €iiO a8 = 2 ME ap: (16) 
dt 


We assume that the co-ordinate system is chosen so that, relative to it, the body is 
at rest in its unstrained state. Then du,/dt = v,, where the v, are the velocity com- 
ponents and the u, are the components of particle displacements from the un- 
strained position. Hence 


OUg # 
Nt * Ua BU B a Va- ( I ( ) 
c 
But, 
OUa,p OU 
at a ot = (Vq 7 Va,kVR) Bs ( 18) 
B 


on account of equation (17). Expanding the right-hand member of equation (18), 
it now follows readily that 


dug g ? 
= Vap — Ya,kYk, Bs 
dt a a 
or 
de l 
ap 
= = €ap — 5 (VakVe,p + Ug 2Ve.a); (19) 


where the e,, are the components of the ordinary elastic-strain tensor. 

It is easily seen that the terms quadratic in the derivatives v, , in the right-hand 
members of equation (19) can be neglected in comparison with the terms egg, since 
the velocities and their derivatives can be assumed to be arbitrarily small for suf- 
ficiently small values of the stresses. Hence from equation (19) we have 

de ag 


= €.2 (approximately). (20) 
dt 


Using equation (20), equation (16) now becomes 


doag m de; buabil deog d a Bini ; 
= f Z = ( Cii pe Z e : 
dt Fie vi dt dt uae ning 
Hence, by integration of these equations, we have 
Tap ne Al iiOas 2ulas + ( ‘ab (21 ) 


where the C,, are constants associated with the trajectory of a moving particle. 
But, if we suppose that initially the body is in its undisturbed state, with the stresses 
and displacements equal to zero, then C,, = 0 for each trajectory. Hence, under 
the above assumptions, considered to be realized for sufficiently small stresses, we 
see that the differential stress-strain relations (12) and (14) reduce to the well-known 
algebraic stress-strain relations of classical elasticity theory. 

' Prepared for the Applied Mathematics Branch, Naval Research Laboratory, Washington, 
D.C. This paper was written while visiting the University of California, Numerical Analysis Re- 
search, Los Angeles, California. 
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2 We have not seen this question adequately discussed anywhere in the literature, and, as far as 
we are aware, it has never been completely resolved. 

3 For a discussion of the absolute time derivative of the stress tensor and related matters see 
T. Y. Thomas, “On the Structure of the Stress-Strain Relations,” these PROCEEDINGS, 41, 716— 
720, 1955. See also C. Truesdell, ‘“The Simplest Rate Theory of Pure Elasticity,’’ Communs. Pure 
and Applied Math., 8, 123-132, 1955, and W. Noll, “On the Continuity of the Solid and Fluid 
States,” J. Rat. Mech. Anal., 4, 3-81, 1955, where references to previous work are to be found. 

4W. Prager and P. G. Hodge, The Theory of Perfectly Plastic Solids (New York: John Wiley & 
Sons, 1951), p. 29. The second set of eqs. (15) constitutes the invariant form of the Prandtl- 
Reuss stress-strain relations, and these relations are strictly equivalent to the noninvariant form of 
the Prandtl-Reuss relations given by Prager and Hodge in view of the interpretation of these latter 
relations proposed by Thomas, op. cit., sec. 2. 

5 The right-hand members of eq. (4) do not represent the most general symmetric tensor in- 
variant F. In fact, Rivlin and Ericksen have shown that the most general such invariant is given 
by a linear homogeneous combination of the Kronecker delta and a set of polynomial tensor in- 
variants of ¢ and o which may be said to constitute a complete set; for the three-dimensional case 
under consideration this complete set consists of seven invariants exclusive of the Kronecker delta. 
See R. 8S. Rivlin and J. L. Ericksen, ‘‘Stress-Deformation Relations for Isotropic Materials,”’ J. 
Rat. Mech. Anal., 4, 421, 1955. This provides the opportunity for further generalization. Here 
we are concerned primarily with the construction of relations of the simplest possible formal 
character consistent with the requirements of this theory. 

6 A derivation of the yield condition from physical or thermodynamic considerations would of 
course be very much worthwhile. In this connection it may be mentioned that a discussion of 
yield conditions based on a variational procedure, which may possibly have some bearing on this 
physical problem, has recently been carried out by us. See ‘‘Determination of the Plastic Yield 
Condition as a Variational Problem,’”’ these PRocEEDINGs, 40, 322-331, 1954. From a purely 
mathematical standpoint it appears that the yield condition is associated with a breakdown of the 
general independence of the stress-strain relations. See ‘Interdependence of the Yield Condition 
and the Stress-Strain Relations for Plastic Flow,’ these PRocEEDINGs, 40, 593-597, 1954. 

7 Instead of the assumptions made in this section, we might take the viewpoint that the co- 
efficients 3\ + 2u and 3& + 2¢ in eq. (6) are at most functions of the x® and ¢, i.e., that these co- 
efficients are not actual scalar invariants of eand a. By the homogeneity condition of sec. 7 it then 
follows that 3A + 2u and 3& 4- 2¢ are material constants. This modification of our assumptions 
would allow for compressibility effects during plastic flow. 

8 It is assumed implicitly here that the quanity o;; and its derivative Do;;/Do do not become in- 
finite. 

® Dr. G. R. Irwin has suggested that the quantity 1 — Cub xp /K in the denominator of the right- 
hand member of eq. (14) be replaced by a power of this quantity in order to take account of the 
differences in speed at which different materials approach the yield state. A corresponding general- 
ization of eq. (14) would also result by allowing the scalars A and B in sec. 6 to contain polynomial 
invariants of the third and higher degrees in the components Tap: We have already mentioned 
the possibility of generalization of relations (12) in n. 5. In view of these possibilities, it would 
appear that our concept of the combined stress-strain relations for the elastic and plastic states 
could be made sufficiently general to account satisfactorily for the observed behavior of most 
materials. However, we feel that, before such generalizations are made, one should explore the 
consequences of the simpler stress-strain relations in the text and that extensions of these relations 
should be based, if possible, on variational principles or thermodynamic considerations. 
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CHROMOSOME STRUCTURE AND THE MECHANISM OF CROSSING 
OVER* 


By R. P. Levinet 
BIOLOGICAL LABORATORIES, HARVARD UNIVERSITY, CAMBRIDGE, MASSACHUSETTS 
Communicated by Ernst Mayr, August 25, 1955 


Chromosome structure and the mechanism by which genes recombine through 
crossing over are intimately related. A knowledge of one should shed light on the 
other. Further, an understanding of the mechanism of crossing over should 
contribute to our conception of the delineation of genetic units within the chromo- 
somes. The principal hypotheses proposed to explain the mechanism of crossing 
over! utilize certain mechanical features of chromosome structure to explain the 
way in which crossing-over occurs. One drawback of these hypotheses lies in 
their not being amenable to strict genetic and cytological proof. 

It has been demonstrated® * that the spontaneous chromosome breakage in 
Tradescantia can be increased approximately seventeen times under calcium- 
deficient conditions. Other cation deficiencies do not have this effect. Further- 
more, it has been shown‘ that the salivary gland chromosomes of Drosophila 
melanogaster can be broken in vitro by treatment with a chelating agent. Because 
calcium and magnesium are major cationic constituents of nuclei,> Mazia‘ has 
suggested that they play a fundamental role in the structure of chromosomes. He 
has proposed that chromosomes are divided into particulate units linked through 
divalent cations, either calcium or magnesium. This model for chromosome 
structure has been extended by Steffensen,* who suggests that deoxyribonucleic 
acid (DNA) molecules, after the structure proposed by Watson and Crick,® 
can be linked in a linear fashion through their terminal phosphate groups, these 
groups being capable of chelating calcium ions. Data to be presented here suggest 
that cations, particularly calcium, may play a major role in chromosome stability, 
perhaps at the sites where calcium ions are presumed to be bound into the structure 
of the chromosome. The validity of a hypothesis considering such a function for 
divalent cations can be tested by attempts to alter, in a specific manner, condi- 
tions to which the structure of chromosomes is believed to be sensitive, and to 
measure the effects genetically, by crossover analysis, and cytologically, by ex- 
amination of the chromosomes themselves. The purpose of this paper is to 
present in summary the results of genetic experiments which support the hypothesis 
mentioned above. Though most of the evidence to be presented is indirect or perhaps 
circumstantial, it points in such a direction as to support the hypothesis under 
consideration. 

Three sorts of experiments on crossing over in D. melanogaster will be mentioned 
here. All of them have been carried out in this laboratory, and certain of them 
have been done independently by Dr. Dale M. Steffensen, of the Brookhaven Na- 
tional Laboratory. Since space does not permit the full development of experi- 
mental methods or the statistical analysis of the results, a detailed presentation 
will be published elsewhere. 

The first genetic experiment is summarized in Table 1. It concerns the effect 
on crossing over of the feeding of calcium chloride to adult females carrying the 


727 











728 GENETICS: R. P. LEVINE Proc. N. A. S. 


X-chromosome markers y ct ras f (yellow body color, cut wings, raspberry eye color, 
and forked bristles). The usual corn meal, molasses, agar, and yeast medium 
was prepared in a 0.1 M calcium chloride solution. The females were permitted 
to remain on this medium for four days, after which they were transferred to 
normal media. Thereafter they were transferred to fresh media every two days for 
a period of sixteen days. Appropriate statistical tests have demonstrated that the 
differences between the various treatments and their controls in this and in the other 
genetic experiments considered in this paper need not be attributed to viability 
differences. In all cases the significance of the results was tested by variance 
analysis, following a procedure similar to that of Levine and Levine.’ Feeding 
of excess calcium results, as shown by a heterogeneity chi square (Table 1), ina 
significant decrease in crossing over between the control and treated series from 
the sixth to the fourteenth day after feeding (see last column, italicized figures). 


TABLE 1 
CRrossING OvER IN FEMALES FED AS ADULTS ON THE UsuAL MEpIUM 
SUPPLEMENTED WITH 0.1  Caucrum CHLORIDE 


Two-Day 


EGG-LAYING No. Per Cent Crossing OveR - 
PERIODS FLIES y-ct ct-ras ras-f Total x? P 
(T* 1,278 14.63 10.48 25.19 50.30) ws 
LAC 1) 433 17.94 10.00 24.91 52.85/ a 
ae » 1847 16.40 10.98 25.77 53.15 : 
Paes jee wae 1 See oper © <9: 0.8 
eee 2°162 16.42 9.30 24.74 50.46 
. 719 Be ius Be an 6So.UeOe 
iT 1771 16.02 10.66 22.18 48.86 a 
4 ic 37 0 ie Ee See ee... ee 
ae > 1846 13.87 9.26 22.64 45.77 a ie 
6 , 1,988 16.30 10.77 23.25 50. 32/ met al 
Beet 054 13.52 10.48 22.64 46.64 va 
pitas Li? eo ss s2 2. = — =o 
ae 5 1179 -12.38 8.73 22.99 44.10) : 
( ( 273 44té«i=S (ss sB.s—«i‘amkwgz CW; «88-5? << 9.0001 
ee. 864 14.94 7.76 25.47 48.17 ; 
wit | 1,011 14.88 oe Ba ea ORM 8 
PRES 624 15.54 11.38 25.63 52.55 — 
' C 870 17.46 9.99 25.17 52.62/ 
iT 357 1820 10.92 24.37 53.49 m 
eee ca BD uo mame. ne °* os 


*‘T”’ represents the treated series; ‘‘C’’ represents the control series. 


It has been noted in small-scale experiments that crossing over is not altered with 
the feeding of magnesium chloride but can be increased with feeding potassium 
chloride. 

The second experiment is one in which larvae were raised on the usual medium, 
which was made 0.01 / for the chelating agent ethylenediamine tetra-acetic acid 
(EDTA). Crossing over was studied in the X chromosome marked with y ct ras f. 
The results are summarized in Table 2. EDTA has the effect of causing signifi- 
cantly increased crossing over during the early egg-laying periods. These results 
are similar to those of Steffensen.2 The feeding of EDTA to adults has not brought 
about significant changes in crossing over.’ 

The final experiments in this series have been an extension of work previously 
reported? in which decreased crossing over was observed between the second- 
chromosome loci black and cinnabar in D. melanogaster as a result of exposing 
pupae to extreme water loss at low relative humidity. Data have been obtained 














Se 





Tlie Ss alle 


—— 


Vou. 41, 1955 GENETICS: R. P. LEVINE 729 


for other loci on the second chromosome as well as for the X and third chromo- 
somes. In each case crossing over was significantly decreased in the first egg-laying 
period when the female pupae had been maintained over concentrated sulfuric 
acid. The results of one experiment are given in Table 3 for crossing over between 
the second-chromosome loci a px or (arc wings, plexus wings, and orange eye color). 
In addition to significantly reduced crossing over, it has been observed that the 
nuclei of the primary oécytes and nurse cells of desiccated animals have been 
reduced about eight times in volume when compared with control animals.” 


TABLE 2 


CROSSING OVER IN FEMALES FED AS LARVAE ON THE UsuAL Mrepium 
SUPPLEMENTED WITH 0.01 MW EDTA 


Two-Day 
Eh GG-LAYING No. Per CENT CROSSING OVER 
PERIODS FLIES y-ct ct-ras ras-f Total x? | 
I 549 17.67 12.02 24.96 54.65 P 
; Aas ao eee a > ( 52 
;' ( 5,158 19.15 11.38 25.47 56.00 — 6 
{3 2,166 18.33 13.29 26.31 57.93 af 
2 ' Fale bi 4 , = rik ¢ 0.006 
¥ ( 8,897 19.06 11.39 24.27 54.72 7.29 any 
a i 3,649 20.80 12.49 24.83 58.12 ‘ e 
» ‘ pa ‘ or 9 -e@ 79? 2 10 0 15 
( 9,874 19.99 11.48 25.31 56.73) 
I 2,580 21.42 12.24 25.46 59.12 wha : 
- bajo Te ecy ‘ < 0.00 
: ( 6,499 18.67 11.05 24.55 54.27 17.38 oes 
: I 2,803 18.99 9.92 25.58 54.49 yee 
‘ . pad : pas - md 9 59 0.002 
6,207 16.22 10.74 24 01 50.97! saa 
. “i 2,114 19.30 10.83 24.02 54.15 ~ or 
6 5) 140 17.66 10.01 25.27 52.94/ 0.84 0.35 
2 T 2,034 18.29 10.27 23.94 52.50 3 r 
d 1954 16.81 10.96 25.03 52.80; °%8 9.8 
. 4 1,621 19.24 9.3 24.42 52.97 ‘ . 
5 C 1,358 16.03 10.27 25.81 52.11 0.34 0.60 
‘y 1,295 17.07 11.05 26.34 4. 46 B: = 
+8] ‘ algo Pallig sie ‘ nage ime a 0.1 
( 3,468 17.76 10.98 25.05 53.79 15 0.70 
T 1, 293 18.71 10.83 25.06 54.60 : 
10 C 3/001 1756 10.64 24.59 52.79) 1.13 0.30 
TABLE 3 
CROSSING OVER IN FEMALES TREATED AS PUPAE TO ZERO PER CENT 
RELATIVE Humipity at 25° C. 
Two-Day 
EGG-LAYING No Per Cent Crossinc Over 
PERIODS FLIES a-px pxr-or Total x? P 
‘ 3,849 1.92 +31 5.23 
‘ oe re ss ; 0.00€ 
( 4, 253 3.34 5.65 8.89 10.84 < 0.0001 
a} 3,767 0.93 1.91 5.84 P 
9 i ‘ 
2 ( 1317 1.27 5.49 6 76 3.02 0.08 
T 3,281 0.85 5.06 5.91 9 . 
P C 3,970 1.18 5.04 6.22 0.24 0.60 
T 3,430 0.55 5.80 6.35 biviad ‘ 
. C 1,393 0.39 $96 5.35 3.43 0.06 


If one interprets crossing over as resulting from genetie exchanges at the sites 
where, as Steffensen suggests, calcium is bound, then one would predict that the 
amount of crossing over or stability of the chromosome should be directly related 
to the amount of calcium present. The data reported here support such a pre- 
diction. Thus, when the amount of calcium is increased by feeding calcium 
chloride to adult flies, crossing over is reduced. The removal of water from the 
primary oécytes by desiccation, with a resultant visible effect on their nuclear 
structure, could result in increased calcium concentration and reduced crossing- 
over. On the other hand, a calcium deficiency, as might arise with the feeding of 
EDTA, results in increased crossing over. 
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Though Steffensen’s model suggests that the crossovers will occur at the sites 
where calcium might be bound into the chromosome, it must be emphasized at 
this point that generalized changes in ionic strength or composition could affect 
crossing over in a manner other than that predicted by the model. Regardless 
of this, the summarized available experimental data as presented here suggest 
that crossing over can be modified by ionic changes according to predictions as 
to their effects on chromosome structure. 

Summary.—The following facts are interpreted as related to a mechanism of 
crossing-over in which it is assumed that genetic exchanges will occur as the result 
of altering cationic conditions (particularly with respect to calcium) to which 
chromosomes may be sensitive: 

1. Decreased crossing over in D. melanogaster when excess calcium is fed to 
young adult females. 

2. Increased crossing over in D. melanogaster with the feeding of EDTA during 
the larval stage. 

3. Decreased crossing over in D. melanogaster following loss of water from the 
nuclei of the ovaries and a possible concomitant increase in ionic concentration. 

4. Increased chromosome breakage in T'radescantia with calcium deficiency 
(Steffensen? *). 

5. In vitro chromosome breakage in D. melanogaster with EDTA (Mazia‘). 

These facts are discussed in terms of a model for chromosome structure, as 
proposed by Steffensen, in which it is assumed that genetic exchanges will occur 
at sites where calcium ions are presumed to be bound into the structure of the 
chromosome. 

* Supported in part by an Institutional Research Grant from the American Cancer Society 
and by research grant No. RG-4000 (C), from the National Institutes of Health, United States 
Public Health Service. 

+ The author wishes to acknowledge the kindness of Dr. Dale M. Steffensen in permitting him 
to quote certain unpublished work. 
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THE GULF STREAM AS AN INERTIAL BOUNDARY LAYER 
By J. G. CHARNEY 
INSTITUTE FOR ADVANCED STUDY* 
Communicated June 8, 1955 


Introduction.—Notable advances have been made in our understanding of the 
ocean circulations in the past few years. The general westward intensification 
of the wind-driven circulation was shown by Stommel! to be a consequence of the 
variation of the vertical component of the earth’s angular velocity. Owing to a 
tendency toward conservation of absolute vertical vorticity, a fluid column moving 
poleward in the Northern Hemisphere acquires a relative clockwise spin, and a 
fluid column moving Equatorward a counterclockwise spin. As a result, the clock- 
wise circulation of the ocean is intensified in the west and reduced in the east. 
Using actual wind stresses, Munk? was able to account for the principal features of 
the mean circulation of the open ocean. The detailed structure of the western 
current remained, however, essentially unexplained. The eddy-frictional mecha- 
nism employed by Munk to provide the necessary energy dissipation for the 
entire ocean was inapplicable to the individual boundary current. In order to 
account for the observed width of the current, Munk was forced to postulate an 
eddy viscosity so large that the eddy sizes were themselves comparable to the 
width. Such eddies are observed in the meanderings of the boundary current 
after it leaves the continental slope but are then cut-off branches of the current itself 
and therefore can have nothing to do with its structure. 

In an analysis of transient wind-driven circulations the author* found that narrow 
intense coastal currents are generated as purely inertial reactions to wind-produced 
shoreward transports of water and raised the question whether the Gulf Stream, 
the Kuroshio, and the Agulhas are not themselves such inertial reactions, requiring 
for their generation neither friction nor local wind stresses.‘ It is the purpose of 
the present article to show that the main features of the western currents in their 
generating areas can be accounted for by treating the motions as boundary-layer 
flows in which only pressure and inertial forces act. Since more is known of the 
detailed structure of the Gulf Stream than of the other currents, the discussion is 
confined to it. Specifically, an explanation is advanced for the Florida Current, 
that portion of the Gulf Stream system in which most of the flow is generated and 
in which the current remains parallel to the coast. After the current leaves the 
continental slope, it may no longer be treated as a boundary layer nor even as a 
steady-state flow. 

Construction of the Model.—The solution to the ocean circulation problem pre- 
sented by Munk is deemed to be valid at the outer or seaward edge of the Florida 
Current. Here, as Munk has shown, the flow is essentially determined by the 
balance of the pressure, Coriolis, and wind-stress forces. As the coast is approached, 
the velocity gradients increase and both the frictional forces and the nonlinear field 
accelerations become large. One may satisfy the boundary condition of vanishing 
normal velocity either by invoking friction or by taking into account the nonlinear 
field accelerations.’ Since the frictional forces are important only at distances from 
the coast small in comparison with the width of the boundary current (unless ex- 
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cessively large eddy viscosities are assumed), whereas the field accelerations are 
large across the entire current, the latter course is adopted. 

An essential property of the Florida Current is its intense baroclinicity. An 
isotherm at the center of the mean thermocline, say the 10° C. isotherm, rises 
from a depth of 900 m or so off South Carolina to nearly the surface at the coast. 
Accompanying the large horizontal density gradients are strong density currents 
which diminish rapidly to zero belowthe 10° isotherm. We adopt as the simplest sys- 
tem embodying the above properties an ocean consisting of two homogeneous incom- 
pressible layers. From a mean of actual soundings we take the density difference 
to be 2.0 X 10-* gm cm~*, and for comparison of theory with observation we 
identify the interface with the 10° isotherm. Disregarding seasonal or intra- 
seasonal variations, we suppose the motion to be steady and to be confined to the 
upper layer. The implied assumption, which is made for simplicity and for lack 
of good evidence to the contrary, is that there is little frictional momentum transfer 
across the thermocline and that the deep waters have been brought to rest through 
the action of bottom friction. 

The Munk theory prescribes the total volume transported by each section of the 
Florida Current. These transports will be regarded as the basic data of our prob- 
lem. Since, however, their computed values depend on a very imperfectly known 
wind-stress distribution over the entire ocean, we shall instead use observed values 
for the purpose of testing our theory. 

Although the coast line in the Florida Straits-Cape Hatteras generating region 
is crescent-shaped, the curvature effects induced in the stream are small. For 
example, that part of the vertical component of relative vorticity which is due to 
curvature is found to be negligible in comparison to the part which is due to shear. 
The motion may therefore be described in a rectangular co-ordinate system in which 
the coast line coincides with one of the axes. We place the y-axis at the coast line, 
or, more precisely, at the edge of the continental shelf, and the x-axis at right angles 
and pointing seaward. With slight approximation we may assume that the y-axis 
points north and the x-axis east. 

Derivation of the Governing Equations.—Referring now to the upper layer, we let 
p be the pressure, p the density, u the z-component of velocity, v the y-component 
of velocity, and h the depth. The Coriolis parameter is denoted by f. Following 
Rossby,® we ignore the unimportant kinematic effects of the earth’s curvature by 
assuming the motion to be planar, but retain the important dynamic effects asso- 
ciated with the latitudinal variability of f by allowing f to be a function of y. The 
equations of motion and continuity for the upper layer then take the simple forms 


o o re) 
u E + y all fv= — (*), (1) 
Ox oe Ox \p 
ov 1 ov ay ) (”) (2) 
u v u=— =], 2 
Ox ae Oy \p ) 
re) re) 
a (hu) + ae (hv) = (0. (3) 


Here, in accordance with our assumptions, no frictional forces are included. 
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If the velocity is to be zero in the lower layer, the pressure must be constant in a 
horizontal plane. Denoting quantities in the lower layer by primes, we have, 
hydrostatically, 


p’ = gph + gp’h’ + constant, 
p = gp(h+h’) + constant, 


in horizontal planes, whence 
) p F 
f + constant = g ( ~ ) h=g'h. (4) 
p p 


Cross-differentiation of the equations of motion and combination with the equation 
of continuity yield the vorticity equation, 


(¢(+ f+ : («+f (Cc+f ( wan *) 
Uu = ie eg . ) = (¢ ) {2 v 
Ox ' OY , ; Ox oy/’ 


a) (‘ + ‘) ioe (' + ‘) z 
U + » =. @, (3) 
Or\ h oy\ A 


where ¢ is the relative vertical component of vorticity 0v/Ox — Ou/dy. This equa- 
tion states that the potential vorticity (¢ + f)/h is conserved along a streamline. 
From the equation of continuity we may introduce the volume transport stream 


or 


function y defined by 


0 
hu= —- ¥ hv = 7 (6) 
OV Ov 


and equation (5) may be written in the form 


(2G) . 2 (1 2H) | ry (7) 
+ = F(p), 
h Lox\h Ox Oy\h Oy 


where F is a function of y to be determined. 

A second equation relating y and fA is the Bernoulli equation. Multiplying 
equation (1) by u and equation (2) by v and adding, we obtain, with the substitu- 
tion of p/p = g’h + constant, 


0 fu? + v? ; 0/u- + vv ; “ Q 
u a r + g’h}+ ie 9 + g’h} = 0, (8) 


l 1 oy\? 1 ow\? 
. ( ‘h = G - 9) 
2 (7 =4 7 (, =) | Tg (Y) ( 


Here G is another function of ¥ to be determined. 
The inertial boundary layer is the region in which the nonlinear field accelera- 


~ 


or 


tions are comparable in magnitude with the Coriolis and pressure forces per unit 
mass. To estimate the orders of magnitude of the terms in equations (7) and (9), 
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we introduce the characteristic magnitudes of u and v, U and V, the characteristic 
width of the current «, and the characteristic length of the current L into the con- 
tinuity equation (3) and the vorticity equation (5). From the first we obtain 
U/V ~ e/L, and from the second V/Le ~ df/dy = 8. Here the symbol “~”’ de- 
notes equality in order of magnitude. Combining the two relations, we get e ~ 
U/s. But U is the characteristic externally prescribed lateral current velocity at 
the outer edge of the Florida Current, i.e., U ~ 10 cm sec™!; L is of the order 10° 
em; and 8 ~ 10-8 em-! sec~!. Hence « ~ (10/10-!%)'? ~ 107 em; V ~ LU /e 
~ 10? em sec~!; and within the Florida Current |v| > |u|, |dv/O x + Ou/dy| ~ 
VL/Ue~ L?/2 ~ 10? > 1. Equations (7) and (9) may therefore be replaced by 


1 /Ov : 
: = F(y) 10 
h es + ‘) (y (10) 





and 


py? 
+ gh = G(y), (11) 
respectively, v being given by the second of equations (6). 
Evaluation of F() and G(); Determination of Empirical Functions.—To evaluate 
F(w) and G(y), we introduce the outer boundary conditions and utilize the fact that 
the highest-order terms in the equations of motion (1) and (2) as well as in equa- 





tions (5) and (8) become negligible at the outer boundary (2 = ©). Denoting 
quantities at this boundary by a bar, we obtain, from equations (2), (4), and (6), 4 
1 dy gq’ dh 
ee Se eS (12) 
hdy f dy 


and, from equations (10) and (11), 


f 
= = F(y), (13) 
h 

g'h = G(y). (14) 


At the outer boundary the quantity y is a prescribed function of y. Hence integra- 
tion of equation (12) gives h as a function of y and therefore also of Y. But, at the 
outer limit, both f// in equation (13) and g’h in equation (14) are also functions of 
y and therefore of y. In this way the functions F(y) and G(y) are determined 
along the outer boundary and consequently at every point in the interior region 
connected with the outer boundary by a streamline. Since in the present formula- 
tion of the problem all interior points are so connected, the problem of the deter- 
mination of F(y) and G(y) is solved in principle. The actual forms of the functions 
will depend on the particular form taken for J(y). 
The parabolic function 


Y= h — vy — yw)? (15) 


corresponds in Munk’s theory to the requirement that the transport into the Gulf 
Stream change sign approximately where the mean zonal wind changes sign. 
The function fits the observed distribution of Y quite well for y < yo, yo being the 
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point at whieh the transport reaches its maximum value Y. To evaluate the con- 
stants y. i+, we take y = Oat the coast, so that Y becomes the volume transport 
of the cuirent, and we take the zero point of y to be midway between the Florida 
Straits and Cape Hatteras (y = yo), a distance of 700 km measured along the coast 
from each. The geostrophically calculated transport in the Florida Straits is 
approximately 30 X 10° m* sec~', and the increase from there to Cape Hatteras is 
approximately 50 XK 10° m* see~!.75 Hence % = 80 & 108 m* sec~', and y has 
the value 50 X 10® + (14 X 10°)? = 2.55 K 10 m sec!. 
It is sufficiently accurate to take f a linear function of y. Setting 


f = fo + Bly — yo), (16) 


we have fy = 0.84 K 10-4 sec—' and 6 = 1.8 X 107—'! m~! see—!. 
With the substitution of equations (15) and (16), equation (12) may be integrated 


directly, to give for h 


») e 

j.9 i. 9 2YJ0 9 ty8 ‘ - 

h? = hy? — ov Cm gh? es | Jom Yo)’, (17) 
] 3g 

or, in view of equation (15), 

: ‘ Me. ‘ tyB (do — P\"” 

h?=h? — — (W—y)+: = (= “) (18) 
y 3g Y 


This expression for / is now substituted in equation (14) and the bar removed from 
vy. We obtain for G(y) 


Sia 2a "cage PS ls 
Gr(y) = 2 hy’ ee : nn — y) + pe - (Wo = y) : . (19) 
g 3g'y 


a 


The function f(y) may be obtained by a similar process, or, more simply, as follows: 
We differentiate equation (14) with respect to y along the outer boundary thus: 


dG dG dy ; dh dy 


= ¢ . 
dy dy dy d dy dy 
and substitute from equation (12) 
dy f dy 


dy g'h dh’ 
The result is 


Gif 
Ta shied Saag (20) 


which holds for all y and therefore for all y.° 

Solution of Equations.—Equations (10) and (11) are solved as follows: Observ- 
ing that y enters only parametrically, we replace x as independent variable by y. 
Then dv/dx = (Ov/dp) (OW /Ox) = hv(dv/OwW), and equation (10) becomes 


) 


d 
hv a + f = hF(¥). 
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§ 
We eliminate v by utilizing equation (11) differentiated with respect to y: 


Ov ig! ; Oh a dG 
"Oy 7 Oy dy’ 
The final equation is 
oh ; 1 /dG 
na J - (' — r) == (), (21) 
Op gh g' \dy 


Its solution, subject to the boundary condition h = h(y) at ¥ = Py), is 


which is recognized as a form of the statement that v is geostrophic. Indeed, the 
above relation could have been derived directly from equations (1) and (6) by 
applying the boundary-layer approximations which were justified in the derivation 
of equations (10) and (11). We have merely shown here that these approximations 
have been used consistently. 

The velocity v is obtained as a function of y and y from equations (i1), (19), and 
(22), and 2 is obtained as a function of ¥ and y by solving the equation 

én (23) 

hv . 

which requires only a numerical quadrature. The boundary condition here is x = 
Oat y = 0. 

Evaluation of he and the Phenomenon of Separation.—Before the integrations can 
be performed, it is necessary to determine the constant of integration in equation 
(19). This quantity is the depth of the upper layer at the latitude of maximum 
current transport. A theoretical value may be obtained by the following heuristic 
argument: We observe first that the flow can no longer be regarded as an inertial 
boundary flow at latitudes greater than yo. This is so because fluid columns which 
enter the stream south of yo and leave it to the north conserve both the stream 
function and the potential vorticity; but the potential vorticity at the outer bound- 
ary increases monotonically with y through yo, whereas the stream function in- 
creases south of yo and decreases north. We must therefore drop the assumption 
that ¢ < f in the expression (¢ + f)/A for the potential vorticity at the outer 
boundary; but this means that the current, defined as the region in which the field 
accelerations are appreciable, must leave the coast at y = yo and go out to sea. 
On the assumption that the current remains coherent because of rotational stability, 
it follows that the shoreward edge must do likewise. Since the interface separating 
the upper. from the lower homogeneous layer also separates the region of motion 
from the region of no motion, we have the result that the interface must rise to 
the surface of the ocean at x = 0, y = yo and thenceforth extend out to sea as a 
front dividing the Gulf Stream from the slope water. That the interface cannot 
come to the surface at the coast for y < yo may be seen from equation (22). Setting 
y = Oand differentiating with respect to y, we obtain, with the aid of equation (12), 
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Aeoas 7 ik 26 p 
Rs ¥ <0. 
dy g 


which shows that / has to decrease monotonically as long as the current continues 


to flow along the coast. 
In view of the foregoing analysis, we may state that h = 0 at x = 0, y = yo. 
Substituting y = Oandh = Oat y = yo in equation (22), we obtain for 4» the value 


= hobo 
hy = ( - ) = $20 m-. 
g 


which compares well with the observed mean value of 900 m given by Iselin'® for 
the depth of the 10° isotherm at the offshore edge of the Gulf Stream in the vi- 
cinity of Cape Hatteras. 

We note that in our theory the phenomenon of separation does not require the 
presence of a protrusion of the coast line such as actually exists in both the Gulf 
Stream and the Kuroshio case. It would be of interest to see whether this con- 
clusion is verified by model experiment. It may be shown that an analogous phe- 
nomenon can be expected to occur in the case of the more easily modeled ocean 


consisting of a single homogeneous layer. 
Comparison of Theory with Observation.—The results of the integrations are repre- 


sented in Figures l(a, b) and 2. Figure l(a) depicts contours of the depth h of the 


upper of the two homogeneous layers at intervals of 100 m, together with the vol- 
ume-transport streamlines y = constant drawn at intervals of ten million cubic 


meters per second. Since the x-scale of motion is very small compared to the y- 
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Fig. 1.—(a) Calculated depth h and volume transport stream function y of the upper layer. 
The contours of h (sclid lines) are drawn for each 100 m, and the streamlines for each ten million 
cubic meters per second. (b) Schematic chart of the observed mean depths of the 10° C. iso- 
therm. 
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scale, the former has been exaggerated in the diagrams by the factor 5. For com- 
parison with observation the depth contours of the 10° isotherm taken from Iselin’s 
article,’ together with his Florida Current temperature cross-sections, were used to 
construct Figure 1(6) a contour chart of the 10° isotherm in which the edge of the 
continental shelf is represented by the y-axis and distances normal and parallel to 
the edge by the x- and y-co-ordinates, respectively. Because of the slope of the 
edge of the continental shelf and because of the complicated thermal structure near 
the surface at the inshore edge of the current, this chart must be considered as 
merely schematic. It will be seen, nevertheless, that the agreement in general 
between the two charts is good. The principal disagreement is in the Florida 
Straits. Here the lack of difluence in the 
* T T lee observed contours, which has no counter- 
part in the calculated contours, is probably 
due to the offshore constraint of the Great 
Bahama Bank. 
The v-velocity profiles across the current 








are shown in Figure 2 for various values of y. 
At a given latitude, the maximum veloci- 
ties are found at the inshore edge of the 
current and reach a value of 4.06 m sec™! 
at Cape Hatteras. These are, as would be 
expected in the absence of boundary fric- 





tion, on the upward side of the observed 
range of 2.00-2.50 m sec~! normally found 
in the swiftest part of the current (Iselin 
and Fuglister!!). The location of the high- 




















- | | est velocities agrees with the observation 
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7 xine) wa 20 that they are “confined toa relatively nar- 


Tet. 2Caelieel weal ote. 5 band 10 to 15 miles wide situated near 
across the Florida Current at various lati- the left or inshore edge of the current.” 
a LA many neti Sh ‘Florida Seralte Application of the Geostrophic Approxt- 
and Cape Hatteras. mation; Concluding Remarks.—The general 

theory which has been presented takes a 
particularly simple mathematical form when the quasi-geostrophic approximation 
is made, i.e., when both wu and v are assumed to be geostrophic and when h varies 
linearly with y: h = ho + u(y — yo). In this case equation (11) becomes redun- 
dant, and equation (10) takes the form 


1] 0 /q’ *) 
of. = Fh). (24) 
h E (f Or/ s| ), 


where 
3 0 rat 0 0 - ho = 
F(h) = ui bo f + Bly — yo) ie fo + Bh 10) a (25) 
h hy + ply — yo) h 
We thus obtain 
07} fB 4 Bho fo8 - 
‘ — iv" h= - (i —f— 6B ) = ; h(y), 
Ox* J b 4] ue gb 
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the solution of which, subject to the condition that h = hatx = ©, is 
h = Ae-*/> + h(y), (26) 


where A? = g’u/f8. The constant A is determined by specifying the volume trans- 
port Yo at y = yo. Wehave, geostrophically, 


. = g’ Oh Lr : 
a } ——_ (h? — h,2), 
. J, ‘fdr af 


where A, is the (necessarily constant) value of h at « = 0. With the condition 


h = h,atx = 0, equation (26) becomes 
h = h(y) — (h — h.je~*™”. 


The width of the current is of the order \ = (g’u/f8)’*. The quantity g’u/f is the 
magnitude |a@ = (g’/f) dh/dy of the lateral velocity at the outer boundary. 
Hence \ = (a@/8)*. It is thus seen that the Gulf Stream is a phenomenon that 
depends in an essential way on the variation of the Coriolis parameter with latitude. 
Substituting observed values for @ and 8, we find that \ is about 50 km. In general, 
it may be said that the geostrophic approximation predicts the rough character of 
the current but does not predict all details. Thus in the geostrophic theory h is 
constant along the coast, and therefore no phenomenon such as separation can occur. 

In conclusion, we mention two extensions of the theory which suggest them- 
selves. First, the boundary-layer technique which has been presented here can 
readily be extended to a continuously stratified ocean. In this case, one prescribes, 
in place of the depth of the upper homogeneous layer, the continuous density 
distribution at the outer edge of the boundary current. Second, an internal fric- 
tional boundary layer should be introduced to reduce the velocities to zero on the 
inshore edge of the currents. This can be done only if one knows something of the 
nature of the eddy viscosity. 

The author wishes to acknowledge his debt to Mr. Henry Stommel, of the Woods 
Hole Oceanographic Institution, from whom he gained much of his knowledge of 
the ocean circulations and in whose inspiring company the Gulf Stream problem 
was discussed over and over again until eventually it may be said to have solved 
itself. 

* This study was sponsored jointly by the Office of Naval Research and the Geophysics Re- 
search Directorate of the Air Force Cambridge Research Center, under Contract No. N-6-ori-139 
Task Order I with the Office of Naval Research. 
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REGULAR FUNCTIONS WITH PRESCRIBED MEASURABLE 
BOUNDARY VALUES ALMOST EVERYWHERE 


By F. BAGEMIHL AND W. SEIDEL* 
INSTITUTE FOR ADVANCED STUDY AND UNIVERSITY OF NOTRE DAME 
Communicated by J. L. Walsh, August 3, 1955 


1. Let Z denote the finite complex plane, and let q be a natural number. We 


shall use k throughout to denote an arbitrary one of the integers 1, 2,...,q, and n 
to denote an arbitrary natural number. Consider a region G ¢ Z bounded by q 
mutually exclusive Jordan curves B),..., B, in Z, with Bs, ..., B, lying in the in- 


terior of B;. For every k and every n, let C,, be a Jordan curve in G, with C,, 
..., Cy mutually exclusive and lying in the interior of C)™, such that (i) B, lies 
in the exterior of C,,™, and C,,4; lies in the exterior of C,,; (ii) for k > 1, B, lies 
in the interior of C,,, and C,,,; lies in the interior of C,,; and (iii) the curves 
CO Pre md. a's q) constitute the boundary of a region G, such that uG, = G. 
Let H, be the region bounded by C,,” and B,. For every k and every ¢ « C,”, 
let J, be a simple, continuous curve whose initial point is ¢, such that (a) J; — 
{¢} ¢ H,™; (b) for every n, J, 1 C, is a single point; and (c) for every z « Hi, 
there is precisely one ¢ « C; withz €J;.! 

THEOREM 1. For every k, let M, be an F, of first category on Cy, and let aj(¢) 
and 8,(¢) be real-valued? functions of Batre class © 1 defined for ¢ « My. Then there 
exists a nonconstant, single-valued, regular f(z) in G such that, for every k and every 
6 « M,, Ri(z) > a(¢) and Sf(z) > B.(¢) as z > B, along J;.° 

Proof: There exist closed, nowhere dense subsets M“ of C: such that M, = 


@o 

u M,,. For every k, m, and ¢, let J, be that subare of J, whose initial 
m=1 
and which, except for this point, lies in H,. Set E, = v 


m=1 


point is J, NC, 
q 

( uv J-™) and FE = uv E,. Then E£ is obviously nowhere dense and closed rela- 

teMm(k) — __ k=l 4 

tivetoG. IfW, =G, u(E G,,,. — G,), then, clearly, for every k, every z¢ W,, 


Y (k) 


that lies in the region bounded by C,, and C,,,»™ can be joined to B, by a Jordan 
arcinZ — W,,. 
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Since a,(¢) is of Baire class <1 on M,, there exist bounded, continuous, real 
functions a, ,(¢) on M, such that, for every ¢ « My, ax ,(O) > arf) asn > @; 
let ¢, = J, 9 C,™, and set ax, ,* (&, )= ax,,(f). Letty,” (7 = 1, 2) be a 
continuous function defined on Z such that y,?(z) = 2 — 7 for z € Gy, vn (2). = 
j — lforzeZ — G,4,, and 0 < y,(z) < 1 for every other z¢ Z. Forze E and 
lying in the closure of the region bounded by C,,” and C,,,,“, let ¢ be the unique 
point on C; such that z eJ;, and define 


* (Kk (1 e (kh 9) / ( 9 
a(z) — fap n(E sn) n (2) + Ak, n+1 (¢ met (Ay (2), (y,(?(z)+ Vn‘ (2)). 


Then a(z) is a finite, continuous, real function on FE such that, for every k and 
every £ € Ch, a(z) > a,(£) as z > B, along J;, and we can obtain an analogous 
function 6(z) for the functions 6,(¢). 

The proof of our theorem can now be completed along the lines of a proof of a 
theorem of Keldysh and Lavrentiev.® 

2. If we make the additional assumptions® that the curves B, are rectifiable and 
that for every w ¢« B, there exists precisely one ¢ « Ci: such that w is the terminal 
point of J, (it will be convenient to denote this /, also by J.), we obtain 

THEOREM 2. For every k, let \,.(w) and u,(w) be real-valued,? measurable functions 
defined for w € B,. Then there exists a nonconstant, single-valued, regular f(z) im 
G such that, for every k and almost every w € By, Rf(z) > Ax(w) and Yf(z) > ux(w) as 
z—>walong J. 

Proof: There exists’ a subset B,* of B,, with meas B,* equal to the length of 
B,., such that A,(w) and yu;(w) are of Baire class <1 on B,*. There exists* a subset 
M, of B,* such that M, is an F, of first category, with meas M, = meas B,*. 
Theorem 2 now follows’ from Theorem 1. 

3. We now make some additional assumptions concerning the curves B,, 
C,™, and J. We suppose that each B, has continuous curvature. Then it can 
be shown that some Green’s function of G has a level line consisting of q com- 
ponents which we call Ci, ..., C,, such that, for every k and every w e B,, the 
inner normal to B, at w intersects C;“ (call this point of intersection ¢) before it 
intersects B, again (if at all), and that, if we let J. denote the segment whose end 
points are w and ¢, the curves J/,, satisfy the assumptions made in sections | and 2. 
Define ¢, to be the point on J, satisfying |¢ — ¢,!/|F — w| = (n — 1)/n. 
Then, as w describes B,, ¢,“ (n fixed) describes a Jordan curve which we call 
C,. Theorem 2 now immediately yields 

THEOREM 3. For every k, let \,(w) and u,(w) be real-valued,? measurable functions 
defined for w € B,. Then there exists a nonconstant, single-valued, regular f(z) in G 
such that, for every k and almost every w € By, Rf(z) > Ay(w) and Yf(z) —> ux(w) as 
z—w along the normal J..." 

!. Let us assume, further, that n = 1, B; is the unit circle, and, for every w e By, 
J. is a radial segment terminating in w.!' Then Theorem 3 implies 

THEOREM 4. Given real-valued,? measurable functions \(0) and u(@) (0 S 0 < 
2x), there exist two real, conjugate trigonometric series in 6 that are Abel summable 
almost everywhere to (0), u(0), respectively.'? 

5. Under the assumptions that are in force in section 2, Theorem 2 for n = 1 
leads to 

TuHeoreEM 5. Given real-valued,? measurable functions p(w) and o(w) on By, there 








~J3 
cae 
~~) 
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exists a nonconstant, harmonic function u(x, y) in G such that, for almost every w € By, ‘ 
Ou/Ox —> p(w) and Ou/Oy > o(w) as x + ty > w along J.. 
Proof: There exists a nonconstant, regular g(z) in G such that, for almost every 
w € By, Rg(z) > p(w) and Yg(z) > — o(w) asz—~walong/.. If f(z) isa primitive 
of g(z) in G, then u(z) = Rf(z) satisfies the conclusion of our theorem. 
6. Suppose once more that the assumptions in section 5 hold and that, for 
every w € B;, J. is a rectilinear segment terminating in w. If we let r(w) be the vec- 
tor ¢ — w issuing from w, 7(w) is a continuous function of w. 
THEOREM 6. Given real-valued,? measurable functions p(w) and oa(w) on By, there 
exists a nonconstant, regular f(z) 7n G such that, for almost every w € Bi, ORf(z)/Or—> 
p(w) and OX f(z) /Or > o(w) as z—> walong J." 
Proof: Denote the length of B; by L. Fix a point wy on Bi, and denote by s the 
length of the are of B, extending in the counterclockwise sense from this point to 


a variable point we B;. The function w = h(z) = |z — w|//|/f — w| + ts (2 € Ja, 
w * wo) is a homeomorphism between H,‘? — J, and the rectangle 0 < Nw < 1, 
0 < Xw < L, preserving measure on B,, and p(h~'(7s)) and o(h~'(7+s)) are meas- 


urable functions on 0 << s < L. According to a slightly modified Theorem 3, there 
exists a continuous (even regular) H(w) in the interior of the aforementioned rec- 
tangle, such that, for almost all s with 0 < s < L, RH(w) > p(h-' (as8)) and 
YH(w) > o(h—! (ts)) as z > &s along the segment Yw = s,0 S Nw < 1. The 
direction cosines, y:(w) and ye(w), of J. are continuous functions of w. For every 
w ~ wo and every z €./., define 


@(z) = { v1(w) RH (h(z)) + ¥2(w) XH (h(z)) | + 
t{yi(w) SH (h(z)) — yeo(w)RA(A(z)) }, 


a continuous function in G — J.,.. Obvious modifications of the proofs of Theorems 
1 and 2 enable us to infer the existence of a nonconstant, regular g(z) in G such 
that, for almost every w « B,, Rg(z) — R&(z) > O and Yg(z) — J&(z) ~ 0 asz> w 
along J. If f(z) is a primitive of g(z) in G, then f(z) satisfies the conclusion of our 
theorem. 


* W. Seidel’s contribution to this paper was made under the sponsorship of the National Science 
Foundation. 

' Note that /¢ may be a spiral tending to By. 

2 The values + © are not excluded. 

’ For a related result ef. F. Bagemihl and W. Seidel, “Some Boundary Properties of Analytic 
Functions,” Math. Z., 61, 186-199, 1954, Theorem 7. 

‘Cf. C. Kuratowski, Topologie I (2d ed.; Warsaw and Wroctaw, 1948), p. 116. 

5 See S. N. Mergelyan, “Uniform Approximations to Functions of a Complex Variable,’’ Am. 
Math. Soc. Translation No. 101, pp. 35-37, 1954. Mergelyan’s polynomial p(z, 2, €) has to be 
replaced by a rational function whose poles lie in the exterior of our B, and in the interiors of our 
By, ..., Bq (ef. Mergelyan’s remark in the middle of p. 24 concerning Theorem 2.3). In case 
our a;(¢) and 8;(¢) are identically constant, we can insure that f(z) is not identically constant as 
in the proof of Theorem 6 in our paper referred to in n. 3. 

6 These assumptions do not imply the existence of a homeomorphism between H,” and a ring, 
under which the curves J, correspond to the radial segments. 

7 Cf. H. Hahn, Theorie der reellen Funktionen (Berlin, 1921), p. 567. 

8 This follows from Lusin’s theorem; cf. E. W. Hobson, The Theory of Functions of a Real Vari- 
able and the Theory of Fourier’s Series, I (3d ed.; Cambridge, 1927), p. 192. 

® Note that F, sets of first category on Bx go over into F, sets of first category on Cc, under 
the homeomorphism between B; and C, induced by the curves Jw. 
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0 This theorem generalizes a result proved by O. Lokki, “Uber das Randwertproblem der analy- 
tischen Funktionen,’ Ann. Acad. Sct. Fennicae, Ser. A-I, 144, 1-8, 1952. 

‘1 Under the present assumptions, the result for harmonic functions corresponding to Theorem 
3 was conjectured by W. Kaplan in a conversation with one of the present authors; a proof of it 
was obtained jointly, and the result appears in W. Kaplan, “Approximation by Entire Functions,”’ 
Mich. Math. J., 3, 1955 (in press). 

12 This theorem extends a result due to Lusin, who showed that if \(@) is finite almost every- 
where, and measurable, then there exists a trigonometric series in @ that is both Abel and Riemann 
summable to \(@) almost everywhere. It is not possible to replace “Abel’’ by “Riemann” in our 
theorem, as follows from a result proved by Germeier. Cf. N. N. Lusin, Integral and Trigono- 
metric Series (in Russian ), (Moscow and Leningrad, 1951), pp. 236, 488 ff. 

138 Tf, in particular, B,; is assumed to have continuous curvature, then Jw may be taken along the 


inner normal to B, at w. 


FUNCTIONAL EQUATIONS IN THE THEORY OF DYNAMIC 
PROGRAMMING. V. POSITIVITY AND QUASI-LINEARITY 


By RicHARD BELLMAN 
RAND CORPORATION, SANTA MONICA, CALIFORNIA 


Communicated by S. Bochner, July 24, 1955 


1. Introduction.—-Given the functional equation 





u(p) = Max [L(u, p, g) + a(p, ¢)], (1.1) 
q 


where L(u, p, g) is a linear operator for every fixed value of g, we wish to determine 
conditions upon the operator which will permit us to write the solution in the form 


u(p) = Max v(p, q), (1.2) 
q 
where v(p, q) is the solution of the linear equation 
v(p) = Liv, p, g) + a(p, q). (1.3) 


We shall also investigate the corresponding problem for the equation 
u(p) = Mm [L(u, p, q, g’) + a(p, 9, g’)] = mM [L(u, p, g, g’) + ap, q,q’)]. 4) 


Here we use the notation Mm = Max Min, mM = Min Max. 
q q’ q’ q 

Special cases of this quasi-linearity of the functional equations of dynamic pro- 
gramming have been noted previously by 8. Karlin in unpublished notes, by L. 
Shapley,! and by the author in connection with the calculus of variations? and in 
unpublished notes. 

After an abstract formulation, we shall present a number of examples involving 
differential and partial differential operators. There is a close connection between 
these questions and the concept of approximation in policy space,*® but space does 
not permit us to discuss this here. A detailed exposition will appear presently. 

2. Abstract Formulation.—Let u(p) be a scalar function of a vector p, belonging 
toaregion R. Let q be a vector variable, an element of a set S, which may depend 
upon p. Similar results to those we state below hold for systems of equations, or 
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vector equations, provided that we define inequality between vectors in the usual 
way, componentwise. The first result we wish to state is 
THEOREM |. Assume that 
a) There is a unique solution of equation (1.3) for any fixed q = q(p), 
denoted by v(p, q), for p in R. 
b) There is a unique solution of equation (1.1) for pin R. 


(2.1) 
c) (Majorization Property) If w(p) = L(w, p, gq) + a(p, q) for a 
fixed q = q(p), all p in R, then w(p) = v(p). 
Under these conditions the solution of equation (1.1) is given by equation (1.2). 
The second result concerns equation (1.4): 
THEOREM 2. Assume that 
a) There is a unique solution of the analogue of equation (1.3) for all 
pin R for any fixed q = q(p), q' = q'(p), denoted bi v(p, q, 7’). 
b) There is a unique solution of Equation (1.4) for all p in R. (2.2) 
eee ae ; 
c) (Majorization Property) If w(p) SL, P, 9, V') + a(pyq, 9’) 
7 . ! , . < 
for fixed q = q(p), g’ = q'(p) and all p in R, then w(p) v(p). 
Under these conditions, the solution of equation (1.4) is given by 
u(p) = Mmv(p, q, g’) = mM v(p, q, 9’). 2.3) 


The proofs are readily obtained as consequences of the assumptions we have 
made. With reference to any particular functional equation, the difficulties arise in 
verifying the majorization properties and in establishing existence and uniqueness. 
Results of this type may be found in other papers by the author.*5 

Examples.—Let us now give some applications of the general principle enunciated 
in the foregoing paragraph. 

Differential equations: It is easy to see that the majorization property is valid for 
the first-order linear equation dv/dt = b(q, )v + a(q, t), v(0) = c. Hence the solu- 
tion of 

du 


— Max [b(q, thu + a(q, 0], u(0) =e u(0) = 2 (2.4) 
q 


(see an earlier paper for sufficient conditions for existence and uniqueness), may be 
written u = Max v(q, 0). 


qg 
Since the Riccati equation, dv/dt = —v? + a(t), v(0) = c, may be written dv/dt 
= Min [q? + 2q + a(d)], o(0) = c, we have an interesting explicit solution of the 
qg 


ticcati equation, namely, 
at t t 
v = Min [ce2 4 4 [eres gee + as) |ds]. (2.5) 
q 0 


Using the equation 


u = Mm [(Aq, q’) + Cl — (Bq, q’)u] = mM [.. .], (2.6) 
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Shapley! obtained a short proof of the result due to Von Neumann that 
(Aq, q’) (Aq, q’) 
Mm ~— = mM : i (2.7) 
(Bq, 7’) (Bq, 9’) 


where g and q’ are probability vectors, provided that (Bq, q’) > d > 0 for all ad- 
missible g and q’. 

We can obtain an extension of these results by considering the differential equa- 
tion 


du 


Be Mm [(Aq, g’) — (Bq, q’)u] = mM [.. .], u(O) =e. (2.8) 
( 


Turning to systems of differential equations, we are confronted by the difficult 
question of verifying majorization properties for systems of the form dz/dt = 
A(q,t)x. The problem was encountered in another connection® and was resolved 
for the case of constant coefficients. A necessary and sufficient condition that e4! > 
0 for all ¢ > 0 is that a;; > Ofor? ¥ 7. For variable matrices it suffices to have 
a,(t) > 0 fora ¥ J. 

Calculus of variations: In an earlier paper? we showed that the problem of maxi- 
mizing J(y) = f{6 F(x, y)dt, subject to dx/dt = G(x, y), 2(0) = ¢, leads to the 
functional equation 


fr = Max [F(e, v) + G(e, v)f.], f(c, 0) = 0 (2.9) 


The corresponding linear equation possesses the requisite majorization property. 
This means that, under suitable conditions in F and G, the local and global varia- 
tional problems are equivalent. 

Parabolic equations: Equations of the form 


u, = Max f(z, t, u, uz, Uzz; Q), (2.10) 
q 


with suitable initial and boundary conditions, arise in the study of various classes of 
multistage decision processes of stochastic type. The problem of majorization for 


the equation 
vy, = f(z, t, v, v2, Vez; @) (2.11) 


arises in connection with stability of solutions of the nonlinear heatequati on™ * and 
is resolved by results due to H. Westphal® and G. Prodi.? 


'L. Shapley, “Stochastic Games,”’ these PRocEEDINGs, 39, 1095-1100, 1953. 

2R. Bellman, “Dynamic Programming and a New Formalism in the Calculus of Variations,” 
these PROCEEDINGS, 40, 231-235, 1955. 

3R. Bellman, ‘“Monotone Convergence in Dynamic Programming and the Calculus of Varia- 
tions,’ zbid., pp. 1073-1075. 

‘R. Bellman, “Some Functional Equations in the Theory of Dynamic Programming. I. Func- 
tions of Points and Point Transformations,’’ 7'’rans. Am. Math. Soc. (to appear). 

5 R. Bellman, ‘Functional Equations in the Theory of Dynamic Programming. II. Nonlinear 
Differential Equations,’’ these ProcEEDINGs, 41, 482-485, 1955. 

6 R. Bellman, I. Glicksberg, and O. Gross, ““‘On Some Variational Problems Occurring in the 
Theory of Dynamic Programming,’ Rend. Circolo Mat. Palermo, Ser. 11, 3, 1-35, 1954. 

7G. Prodi, “Questioni di stabilit’ per equazioni nonlineari alle derivate parziali di tipo para- 
balieo,”’ Att? Accad. nazl. Lincet Rend., Classe sci. fis. mat. e nat., 10, 365-370, 1951. 
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* R. Bellman, “On the Isxistence and Boundedness of Solutions of Nonlinear Partial Differential 
Equations of Parabolic Type,’ Trans. Am. Math. Soc., 64, 21-44, 1948. 

* H. Westphal, “Zur Abschiitzung der Lésungen nichtlinearer parabolischer Differentialgleichun- 
gen,” Math. Z., 51, 690-695, 1949. 


EXISTENCE OF FUNCTIONALLY INDEPENDENT AUTOMORPHIC 
FUNCTIONS 
By 8S. BocHNER AND it. C. GUNNING 
DEPARTMENT OF MATHEMATICS, PRINCETON UNIVERSITY 
Communicated July 11, 1955 


In the theory of automorphic functions of n complex variables, the existence of n 
functionally independent meromorphic functions is demonstrated, separately, in 
the multiperiodic case by means of theta functions, in the case of a bounded domain 
by means of Poincaré series, and in the modular case by means of Eisenstein series. 
We have found that these constructions differ but little from one another, and we 
are going to established the existence theorem in a unified version which will sub- 
sume all those known to us. 

1. We take a complex analytic manifold D of complex dimension n, and a 
countable group TI’ of analytic homeomorphisms of D onto itself, transformations 
with fixed points being permitted. A meromorphic function f(z) on D is auto- 
morphie if f(7z) = f(z) for al! T «I. 


Definition 1: A set of factors of automorphy' * is a set of meromorphic functions 

nr(z) on D indexed by the elements 7’ ¢ T and satisfying 
nst(z) = ns(T'z)nr(z2). 

Definition 2: A countable collection of sets of factors of automorphy {n,. 7(2)}, 
a = 1, 2, ..., is called a complete system of factors if the following conditions are 
fulfilled : 

i) For each a there is a holomorphic function h,(z) # 0 such that the functions 
h,(2)|n,. 7(z)]~! are holomorphic. 

ii) For each a there is a subgroup [, of T such that (a) 9, -(z) = 1 whenever 


T « T,; (b), letting A, be a set of representatives of the right cosets of T, in I, 
the series 
DL ha(2) [na, r(2)]~ (1) 
Te Ae 


is absolutely uniformly convergent on any compact subset of D; (c) the sum function 
(1) is not identically zero. 

iii) There isa fixed set M c¢ D (M being a countable union of closed subsets of D 
‘ach of which is of topological dimension at most 2n — 1) such that, for any z° « D — 
M, each function n, ,(z) is holomorphic in a neighborhood of z°, and the point set 


V(z°) = tz] ng, r(2) = ma, r(2°), all T ¢ T and all a} (2) 


is zero-dimensional in a neighborhood of 2°. 
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Note that in definition 2 the convergence of the series (1) is independent of the 
choice of the representatives of the various cosets. The set (2) may also be defined 
by 


V(z2°) = {2 Na, r\2) = Na. 7(2°), all T « A. and all e}. 


Assuming that {»,, ;(z)} form a complete system of factors of automorphy, we 

may introduce the functions 

r A = vA —k 

Pa, K\) a >» (Ne, r(2) | 

T € Ag 

for all integers / > 1. These functions are meromorphic on D, and ¢,. ,(z) # 0. 
It is also clear that they are independent of the choice of representatives of the 
various cosets and that they satisfy the functional equations 


Pa,~(T2) 2 (Ne, (2) }* eq, x(2) 
forall 7 ¢« T. Therefore, 
fa, (2) = Pa, (2) (Gq, (2)]~* (3) 
are well-defined automorphic functions. 
Analytic functions w;(z), ..., W»(z), m <n, are functionally independent on D 
if the Jacobian matrix 
Oe. i ae) 
: (4) 
One 2) 


has its maximal rank m at some point of D. 

THEOREM |. Jf the manifold D and the group T admit a complete system of factors 
of automorphy, then there exist at least n functionally independent automorphic func- 
tions. 

Proof: From the collection of automorphic functions (3) we select a maximal set 
of functionally independent functions 


Wi(z) = fas, x, (2), + -) Wm(Z) = fam, km(2); (5) 


Ja, 


that is to say, the functions (5)are functionally independent, but, by adjoining to 
the set (5) any additional functions from (3), we secure a set of functions which are 
not functionally independent. We shall assume that m < n, and from this derive a 
contradiction. 

The set M’ of all points of D at which the Jacobian matrix (4) is of rank less than 
m is a closed subset of D of topological dimension at most 2n — 2. The sets 

M, = {z| ¢,. ,(2) = 0} 
and 
a 


M,' = {z|h,(z) = 0}, 


where h,(z) are the functions introduced in part i of definition 2, are also closed sub- 
sets of D of topological dimension at most 2n — 2. Letting be the set introduced 
in part iil of definition 2, the set 


N = ™M uM’ vu (U.M,) v (U.M,’) 
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is of topological dimension at most 2n — 1, so that there isa point 2° «D—N. We 
shall then consider the subvariety 
V(z°) = jz | w,(z) = w,(2z°), ..., Wm(Z) = w,(2°)}. 


Since the matrix (4) is of rank m at 2° 


, we may introduce the new co-ordinates 


is Sales Rea scape 
in a neighborhood of 2°; then locally V(z°) will be a linear subspace of topological 
dimension at least 2. For any other function f, ,(z), the Jacobian matrix 
O( fa, e(Z), Wi...) Wm) 
O(w, » ey Wm Smt) -- + 2n) 
is by hypothesis of rank m, which implies immediately that 
Of, x(2) Of, (2) 


OZm+1 Oz 


= () 


ll 


n 


in a neighborhood of 2°; thus all the functions (3) are constant on the subvariety 
V(z°). 
Following Siegel,* we introduce for each a the function 
¥.,2) = DD OF — e0,:)me, r(2)]-", (6) 
T ¢ Ag 
noting that it is a meromorphic function of the point z on D and of the variable ¢. 
As a function of ¢, for fixed z, V,(¢, z) has simple poles at the points 


= 


£ = Ga, 1(2) Ma, r(2) (7) 


° ° >" 4° . 1 . 

for the various T’ « A,. Since, by our construction, the numbers | ¢,, ,(2)ng, 7(2) 
are bounded from below by some positive number 6 for all z in a neighborhood of z°, 
the function (6) may be expanded as a power series, 





V.(f, 2) ss Cheat > Su, Mae" is” (8) 
(=1 


for z near z° and [¢| < 6. 

Restricting the variable z to the subvariety V(z°), we have proved that the 
coefficients in the power-series expansion (8) become constants, so that the function 
(6) is independent of z « V(z°). In terms of the poles (7), this means that 

- sara ° ° 
La, 1(Z) Ne, 7(z) De. (2 )Ne, r(z ) 
forallze V(z°). If T €T,, this equation becomes simply 


*} 


Pa, (2) — Pu, (2 


so that we must have indeed 
4? 
Na, 7 (2) = Na, (2 ) 
for all z e V(z°). This contradicts the crucial condition iii of definition 2, and the 


theorem follows from this contradiction. 
2. Our first application shall be to the multiperiodic case. We assume that the 








- 
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period matrix (the nm X 2n matrix whose columns represent vectors generating the 
invariant lattice) is a Riemann matrix;* * © by a standard normalization, we may 
assume that this matrix is of the form 


(IQ), 


where / is a square unit matrix, Q = (w,,) is a square matrix, and there is a diagonal 
matrix A = (6;;) with integer entries such that AQ is symmetric and 7A(Q — Q) is 
positive definite. The transformations of the group [ may then be written ex- 
plicitly in the form 


n 

ryy 4 U 3 ti = / 

T'(m, ..., Mg; My, ..., My): 2, 2, + Mm, + > WM, 
k=1 


and as factors of automorphy we shall take the functions 


NT (m, mn; m,’, — 


imns) (2) 


n n 
T 
exp ; +> M 6 ;;(22; — w5;) + i 6 jj j:,M My, t. (9) 
j=1 


; j,k=1 

THEOREM 2. [fin the multiperiodic case the period matrix is a normalized Riemann 
matrix, the factors of automorphy (9) form a complete system of factors. 

Proof: (i) ‘This is an example of a case in which a single set of holomorphic 
factors of automorphy is sufficient. The first condition of definition 2 is trivially 
satisfied by taking h(z) = 1. 

ii) Let Tl, (1, with a = 1) be the subgroup of T consisting of the transformations 

TO ov OP te 
and select as the coset representatives A, the transformations 
To(m,, ..., Mn) = T(m, ..., ma; O, ..., 0). 


On any compact subset of D the series 


LD [nr(z))-! (10) 


Te At 


is eventually dominated by the series 


sided ” 
te exp !/,77 ) Dd §;0,.mm, A 
mi, mn=—« ' 7k=1 f 


which converges as a consequence of the hypothesis that the period matrix is a 
normalized Riemann matrix. Since (10) is the Fourier expansion of the sum 
function, the latter is not identically zero. This verifies the second condition. 


ili) Suppose that for some points 2° and z, 


a aaa! ~O 
NT y(m, mn) (2) or NT (my, _my) (2 ) (11) 


for all To(m, ..., m,). Using the explicit form (9) of these factors of automorphy, 
condition (11) means that 


n n 

= ate 
>. M652; = tS Mj;2; + g(m, ..., My) 
1=1 =] 
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for all mj, ..., m, and some integers g depending upon the indices m, ..., m,; 
in particular, by setting m, = | and m, = ... ie he SM, SM, 


we secure 
Yr 
° ‘ 
eee (12) 


for some integers g,. Whatever the choice of 2°, the set of all points z satisfying a 
condition of the form (12) with any arbitrary integers g is surely zero-dimensional at 
z°. This then completes the proof of the theorem. 

3. We next consider the case in which D is a bounded subdomain of the complex 
affine space and T is any properly discontinuous group of analytic homeomorphisms 
of D onto itself. Letting J7(z) be the complex Jacobian determinant of the trans- 
formation 7 and h,(z) be the set of all not identically vanishing polynomials with 
rational coefficients, we introduce the factors of automorphy 


nap. (2) = [dr(2)]-2*[h, (12) ] hy (2) (13) 


for all integers R > 1. 

THEOREM 3. J/f D is a bounded subdomain of a complex affine space and T is a 
properly discontinuous group of analytic homeomorphisms of D onto itself, then a sub- 
collection of the sets of factors of automorphy (13) form a complete system of factors. 

Proof: (i) As the functions h,(z) of part i of definition 2, we may take the 
polynomials h,(z) appearing in (13). 

ii) Taking l, = land A, = PF for each a, it is well known* 7 that the series 

¥ [J r(z) |? 

Ter 
is absolutely uniformly convergent on every compact subset of D; since each poly- 
nomial h,(z) is bounded on D, the series (1) likewise converges in this case. 


° 


Now select a pointz° « D with rational co-ordinates, and let 


I’ = {7 | T €f and |J7(2°)| < 1}; 


there will be only finitely many transformations J, 7), ..., 7, in T — I’, including 
the identity transformation 7. Considering only those polynomials h,(z) for which 
h(2°) = 1, h,(7:z°) = ... =h,(T,2z°) = 0, (14) 


we have 
a RW2°) = 1+ DY h,(T2°)[Jr(z°)]* ¥ 0 
rar’ 

whenever R > R(a) for some suitable index R(a) sufficiently large. We shall con- 
sider only those factors of automorphy (13) for which the component polynomials 
h,(z) satisfy (14) and for which R > R(a): condition ii is thus verified for this 
subcollection. 

iii) Let N = {z|z2 = T2° forsome Te T}, N, r = (2|h,(Tz) = 0}, Nr be 
the set of all fixed points of the transformation 7’, and 

M=Nvu U NrvUN,.-. 
Tel, T¥1 Tel,a 


We may of course always assume that the group T is nontrivial. Suppose, then, 








{ 
' 
} 
} 





; 
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that for points 2° « D — M and z' « D, 


Na, R, r(2') = a, R. st (15) 


for all factors of automorphy ,. x. 7(z) admissible under the restrictions imposed 
in part ii above. Since (15) must hold for all R suitably large, we have 


h(2)h,(T2°) = ha(z°)h,(T2') (16) 


for all polynomials h,(z) satisfying (14). If z' is not congruent to z° under TP, 
select any T ¢ Tl’ other than the identity transformation; we may then construct a 
polynomial h,(z) satisfying (14) and such that |h,(z°)| < 1, |h,(T2°)| > 2, and 
Vy < |h,(2")|, |h,(T2')| <1. This clearly contradicts (16) and is thus sufficient to 
demonstrate condition iii of definition 2 and complete the proof of the theorem. 

4. For our final example we shall consider Siegel’s modular group of degree p, 
where n = '/p p(p + 1). If A is a matrix, then A > 0 means that it is positive 
definite, and ‘A denotes its transpose. The manifold D is the space of all symmetric 
p X p complex matrices Z such that 


ha = 
(Z — Z)>0, 
22 
$ 


and the group T is the set of all transformations 
T: Z—~ (AZ + B)(CZ + D)—, 


where A, B, C, D, are p X p integer matrices satisfying 


‘(A B 0 / A B\ _ 0] (17) 
CD ie CD A=- oe ; 


The transformations 7 thus correspond to the 2p X 2p integer matrices satisfying 
(17). We then introduce the factors of automorphy 


nr(z) = det (CZ + D)?*?, (18) 


where ‘det’? denotes the determinant of the matrix; these functions are holo- 
morphic and nonvanishing on D.*% 

TueoreM 4. For the modular growp of degree p, the factors of automorphy (18) 
form a complete system of factors. 

Proof: (i) The first condition is trivial, as in Theorem 2. 

ii) Let T, (0, with a = 1) be the subgroup of I consisting of all transformations 


whose corresponding matrices are of the form 


‘U SU 
0U-) 


where S isa p X p symmetric integer matrix and U is an integer matrix with det 
U’ = 1. Then, for any system <A, of representatives of the right cosets of T; in I, 
it is known® ® that the series 

>» (nr(z)]-! = by det (CZ + D)-?-? 

Te Ai T ¢€ Ai 
is absolutely uniformly convergent on any compact subset of D. A Fourier ex- 
pansion’ shows that the sum function is not identically zero. 
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iii) Let us note that for an arbitrary p X p symmetric matrix A, the matrix 


O07 ; 
9) 
oe 4 (19) 


satisfies (17) and therefore corresponds to a transformation 7’ « [. Suppose, now, 
that for two points 2° and 2! of D, nr(z'!) = nr(z°) for all T ¢ T. For the particular 
transformations (19) this may be written explicitly as 


det (A — Z°) = edet (A — Z'), (20) 
where «?*+? = 1, and A is an arbitrary symmetric integer matrix. Since (20) is a 


polynomial relation in the variables A which is true for all integer values of its 
arguments, it must be true for all symmetric matrices, and, further, « = 1. This, 
in turn, is only possible when z° = z', so that condition iii of definition 2 is verified 
in this case also, and the proof is complete. 
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SOME SPECIAL EQUATIONS IN A FINITE FIELD 
By L. Caruirz anp H. H. Corson 
DUKE UNIVERSITY, DURHAM, NORTH CAROLINA 
Communicated by H. S. Vandiver, August 25, 1955 
Let N(a1, dz, ..., a,) be the number of solutions in a finite field GF(q), ¢ = p”, of 
Q42,"' + Aote” +... + 4,2," =. (1) 
tecently, investigati fN has claimed particular attention in papers 
vecently, investigation of V(a;,..., a,) has claimed particular attention in papers 
by L. K. Hua and H. 8. Vandiver,! H. Davenport and H. Hasse,? N. C. Ankeny,’ 
and others. We shall state some further theorems concerning this same problem: 
in particular, dealing with the variance of N(a, ..., @,) over the coefficients a, .. . , 
a, of equation (1). 
We shall need the following result: 
Lemma. I[f Ao ts the average of N(a,,..., a,) over all nonzero a, in equation (1) with 


c = 0, and if A, ts the average of N(a, ..., a,) over all nonzero a, in equation (1) with 
c = |, then 












a aE SR oe 


a 
= Se 
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Ay = A, = 7" 
First, consider the case where c = 0 in equation (1). Also, let m, me, . . . , m, be 
positive integers which, it may be assumed, divide g — 1, and let 
ee MM, . . . Mi 

A(124, te, 2 0 yp OE) = i ~ ’ 

[m;,, Ma, ---, My] 

where the denominator is the least common multiple of the m,. 
THEOREM |. Jf c¢ = 0 in equations (1), then the variance of No = N(a, dz, .. . , 


a,) over all nonzero a; is 


2B (No — Ao)? = q—(q — 1)'*? ys > (-1)’ | areas 5 3 
ai#0 k=0 (i) 
The inner sum is over the r indices taken k at a time, and (0) = 1. 

The determination of Ay and >> N» follows from computing >> No, summed over 
all a;, and then using the inversion formula which is applicable. To find >> N,?, 


consider 


1 m 
an" +...+42, = 0, 
m 
ayy +...+ ay,’ = 0. 


If M is the number of solutions of equations (2) in GF(q), then > M = DON,?. 
Treating equations (2) as a linear system in the unknowns q, a, . . . , a,, the prob- 
lem is then to find the number of matrices in the x;"" and y,”' of rank 1, x; and y; 
in GF(q). 

The next theorem is an immediate corollary of Theorem 1: 

THEOREM 2. A necessary and sufficient condition that No, the number of solutions 


(2) 


of equation (1) in GF(q) with « = 0, be independent of a;, ... , a, (a; # O) ts that 
os 7 Sa ix) = 0. (3) 
k=0 (i) 


As several solutions for equation (3) are known,':* the determination of all 
solutions may be of interest. This result is contained in the following theorem: 

THeoreM 3. Let S be the set of m;, the exponents in equation (1). Then a neces- 
sary and sufficient condition that No, the number of solutions of equation (1) with ¢ = 
0, be independent of a, ..., a, (a; # 0) 2s that S, or some subset R of S with each 
member of R prime to each member of S — R, has either of the following properties: 

a) One member of the set is prime to all other members. 

b) The set has an odd number of elments and each member of the set is of the form 2g, 
and the g; are prime in pairs. 

These results have suggested various problems. We have proved, for example, 
the following theorem: 

THEOREM 4. Let m; = vs; in equation (1) with c = 0, where (s;, 8;) = 1, % #7. 
Then No equals ihe number of solutions of 


ar; + at.” + ...+ 4,2,” = 0. 


Hua and Vandiver' proved an equivalent result when v = 1. 
Theorem 3 implies the following extension of a result of Carlitz: 
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THreoreM 5. If N ts the number of solutions of 


9 9 . ws 
ay + ...+42," = by" +... + Oy,”’, 


with (2u;,w;)) =1(i@=1,..., r;j=1,...,8), and (u;,, uj) = 1,7 4 J, and r odd, 
then 
‘uot 
Let us now assume that ¢ # 0 in equation (1); we may suppose that ¢ = | with- 


out loss of generality. In the same way as in Theorem 1, one can prove the follow- 
ing: 


TuHeoreM 6. /fc¢ = 1 in equation (1), then the variance of Ny = N(a,..., al;) 


over all nonzero a; ts 


Wig ss Mh, ib. 
} 


> (Ni — Ai)? = q’? (q — 1) dom, .. 2m, — (q — 1) — 
ie sss 5 yet 


aixO0 \ 
A, ts the average N, defined in the lemma, and |m,..., m,| ts the least common multi- 
ple of the m,. 

Corresponding to Theorem 3 is Theorem 7: 

THeorEM 7. For N, to be independeni of ay, a2, .. . , a, (a; # O) tn equation (1), 
all m; must be 1. 

These results will be developed more fully elsewhere. 

! These PROCEEDINGS, 34, 258-263, 1948; ibid., 35, 94-99, 1949. 

2 J. reine u. angew. Math., 172, 151-182, 1935. 

3’ These PRocEEDINGs, 40, 1072-1073, 1954. 

' Pacific J. Math., 5, 177-181, 1955. 


THE LOCAL STRUCTURE OF THE ARTIN ROOT NUMBER* 
By B. Dwork 
ey AOE ane ae 
Communicated by O. Zariski, July 25, 1955 


Let k be an algebraic number field, K a normal overfield of finite degree, X a 
character of the Galois group, G(K /k), of K over k, and let W(X) be the root number 
appearing in the functional equation of the Artin L-series, L(s, X, K/k).' We 
shall, in this preliminary paper, outline a proof of 

THEOREM ©. W(X) has (modulo +1) a decomposition into local terms which 
depend only upon the local behavior of X. Viewed as functions of characters of local 
Galois groups,” these local terms have (modulo +1) the usual linearity, factor group, 
and induced character properties. The qualification “modulo +1” may be dropped 
if K is of odd degree over k. 

1. A set of finite groups is said to be a family if it is closed under the process of 
forming subgroups and factor groups, groups being identified by the second iso- 
morphism law whenever it is applicable. A function defined on the set of linear 
characters of the groups in a family, A, and taking its values in some fixed (multi- 





aS 
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plicative) group is said to be extendable with respect to A (or A-extendable) if it can 
be extended to a function on the set of all characters of groups in A with linearity, 
factor group, and induced character properties’ and the property of invariance 
under inner automorphisms of groups in A. 

Extendability can be characterized in an important case: 

THEOREM 1. A function, F, detined on the linear characters of a family, A, of 
nilpotent* groups is A-extendable if and only if it has the factor group property, ts in- 
variant under inner automorphisms of groups in A, and satisfies the following conditions: 

a) If @; (¢ = 1, 2) ts a linear character of a proper subgroup, Gi, of G « A which 
induces a given irreducible character of G and if both G, and Gs properly contain H, 
an invariant Abelian subgroup of G, and G/H is Abelian of type (p, p) (p prime), then 
F(0,) = F(6.). 

b) If@%s a linear character of G ¢ A with restriction, 0, to H, an invariant Abelian 


p 
subgroup of prime index, p, then F(0@) = Il F(@X,), Xi, ..., X, being the linear 
t | 
characters of G which are trivial on H. 

2. The known, analytically derived properties of the Artin root number are 
applied by means of 

THrorEM 2. Let F be a function defined on all linear characters of local Galois 
groups with the following properties: 

a) F is invariant under character transformations corresponding to those group 
isomorphisms which are derived from topological field isomorphisms. 

b) If X is a linear character of the global Galois group, G(K/k), then F(X,) = | 
for almost all primes, y, of k, X, being the character of the local Galois group, G(Kk,/ky), 
obtained by restricting X to a y decomposition subgroup of G(K/k). 

c) If, in the notation of property b, M is defined on the set of linear characters of 
global Galois groups by M(X) = Il F(X,), the product being over all primes of k, then 


¥ 
M is extendable with respect to every family of global Galois groups. 

d) F ts extendable with respect to every family of nilpotent local Galois groups. 

Then it may be concluded that F is extendable with respect to every family of local 
Galois groups. 

This result is a consequence of the following existence theorem. 

THEOREM 3. Let A be a normal noncyclic extension of a local number field, B. 
If m (+1) ts the degree over B of an intermediate field, then there exists an algebraic 
number field, k, with normal overfield, K, such that B is isomorphic to the completion of 
k at each prime which has just one prime divisor in K, the corresponding completion of 
K is isomorphic to A, and the number of these primes of k is either m or one. 

3. If @is a character (i.e., a bounded, continuous® homomorphism into the com- 
plex numbers) of the multiplicative group, B*, of a local number field, B, let R(@) 
be the root number of @. Specifically: R(@) = 1 if B is Archimedean, with the 
exception that R(@) = —7 if (—1) = —1, Breal. If B is y-adic, then R(@) = 
7(6)/|7(0)|, where 7(@), the gauss sum in 8, is as given by Hasse (op. cit., sec. 6, eq. 
[16]). The arithmetic methods described by Hasse® give’ 

TuHeoreM 4. Let B be a local number field; 

a) If A tsa cyclic overfield of prime degree and ® is a character of A* obtained from 
a character of B by composition with the relative norm, then [ I] RO) | R(®) = V(A/B), 


46.0 
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a fourth root of unity, independent of 8, whose square is | Tl 6(—1 | 0(— 1), the prod- 
6\e 


ucts being over all characters of B* from which 8 may be obtained by composition with 
the relative norm. V(A/B) = 1 tf degree A/B ts odd. 

b) If A isan Abelian overfield of B with Galois group of type (p,p) and if 0; (7 = 
1, 2) «s a multiplicative character of A;, an intermediate field of degree p over B, such 
that both 6, and 6. divide (in the sense of part [a]) the same character of A*, while 6, is 
not divisible by any character of B*, then 


R(0,)V(A,/B) = R(@.)V(A2/B) (1) 
6:(—1) [R(A,) |? = 62(—1) [R(2) ]*. (2) 


t. Theorem O follows from Theorems 1, 2, and 4: Let F’ be the function on 
linear characters of local Galois groups obtained from R by the identifications of 
local class field theory. Let F be the function X — [F’(X) ]?X(—1) defined on the 
same set of characters. It follows from Theorems 1 and 4 that F is extendable 
with respect to all families of nilpotent local Galois groups. The function, M, in 
hypothesis (c) of Theorem 2 is now the mapping X —~ [W(X)]*, X being a linear 
character of a global Galois group. It follows from the analytically derived proper- 
ties of W(X) and from Theorem 2 that F is extendable with respect to all families of 
local Galois groups. The last sentence of Theorem O follows in the same way, 
using F’ instead of F. 

A purely local study of the arithmetic structure of the Artin root number is now 
possible (cf. Hasse, ep. cit., sec. 8). 

* This paper is the result of an investigation suggested by E. Artin. It was carried out with 
the advice and encouragement of John Tate. 

' KE. Artin, “Zur Theorie der L-Reihen mit allgemeinen Gruppencharakteren,’’ Abhandl. Math. 
Sem. Hamburg, 8, 292-306, 1931. 

2 “Global” and “local’’ are used to distinguish between algebraic number fields and their com- 
pletions under a valuation (Archimedean or non-Archimedean ). 

’ Compare with eqs. (8)—-(10) in Artin, op. cit. 

‘ A similar result for supersolvable groups will not be needed. 

> The topology is as given by the valuation. For our purposes it may be further assumed that 
the kernel of @ is of finite index. 

° H. Hasse, “Artinsche Fiihrer, Artinsche Z-Funktionen.. .,’’ Acta Salmanticensia, 1954. 

’ Part (a) of Theorem 4 was stated independently and paritally verified by Hasse (zbid.). Fora 
full proof see B. Dwork, “On the Artin-Weil Root Number” (unpublished doctoral dissertation, 
Columbia University, 1954). 


THE DIVISION PROBLEM FOR DISTRIBUTIONS* 
By Leon EHRENPREIS 
Communicated by C. Chevalley, August 2, 1954 

In a previous paper! we gave a partial solution of the following problem: Let 
D be a partial differential operator, and let 7 be a distribution ;? can we find a dis- 
tribution S such that DS = T? We had previously shown the existence of such an 
S in the case that 7’ is a distribution of finite order. In the present paper we shall 

complete the solution and show that we can find an S for an arbitrary 7’. 





=] 
or 
“J 
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The method of proof is to give explicitly the topology of the space D which is the 
Fourier transform of the space ® of L. Schwartz.?:* Once this is done, the main re- 
sult is proved by essentially the same methods as in our previous paper. We give 
below an outline of the proof; the details will be published elsewhere. 

For the description of the topology of D, we consider only the case where D is the 
space of indefinitely differentiable functions of compact carrier on the real line; the 
case of higher-dimensional Euclidean space is slightly more complicated but is 
handled in exactly the same manner. 

THEOREM |. The topology of D can be described as follows: We choose, first, post- 
tive numbersc,n. Next, for each integer 1 = 0, we choose 

a) Two positive integers at, dj with } at}, \di} monotonically increasing to infinity 


(aj, di are detined for t = 1 only.) 
b) A positive integer by = 5 such that }bi} is monotonically increasing. 
We consider the N of all G ¢ D which satisfy, for each t, 


max jexp (—3di \y\)z" G(z)| Sn 


zel'y™ 


(1) 


jor eo Gol 2... bi, where z = x + ty, ot = o> ts the set of z for which ly} < 
c, and the regions Ty~ for t > 0 are defined as follows: Let yi~ be the curves defined by 


the equations 


Hi t= eM! for ti = I, (2) 
Dh y = U for yy ia cad 
i i io a 
th) that is, 
HH 
i bi log | z 
i y ) 4 for Yt SF (3) 
| \ { d} 
ah 
At 0 for ae hal 
ai (yit ts the upper curve, yt~ the lower curve). Then Ty~ is the set consisting of all com- 
Hy plex numbers of the form 
‘ 
ue 2=2' + iat + 2%, Zev =» —ce<it<e. (4) 


Assertion: ‘The sets N described above form a fundamental system of neighbor- 
hoods of zero in D. 

Outline of Proof: It is not difficult to see that the sets N form a fundame tal 
system of neighborhoods of zero for a locally convex topological space and that 
each N is a neighborhood of zero in D. The main tool used in proving Theorem | is 
the following lemma, the proof of which is not difficult: 

LemMa. With yi~ defined as in the statement of Theorem 1, let 6¢~ be the curves 
yt + iat. Then, for any G ¢ D, we have, for any real t, 


‘ LG) exp (itz) dz = f°. G(z) exp (itz) dz, (5a) 
Sf, G(z) exp (itz) dz = f2. G(z) exp (itz) dz. (5d) 


Once this lemma is established, Theorem | then follows by a method similar to 


that used in a previous paper.® 
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Using Theorem 1 (or, rather, the analogous result for higher-dimensional Euclid- 
an spaces) we deduce* ® 

THEOREM 2. Let D be any partial differential operator with constant coefficients, 
and let T be any distribution. Then we can find a distribution S such that DS = T. 
This result is also valid if D is a partial differential-difference operator with constant 
coefficients. 

Remark: We can also give a direct description of the topology of the space E’ 
which is the Fourier transform of the space of distributions of compact carrier. * 
This leads to another proof of the result® that, if D is a partial differential-difference 
operator with constant coefficients, and f is an indefinitely differentiable function, 
than we can find an indefinitely differentiable function g such that Dg = f. By 
similar methods we can also show that, if f is any entire function, there exists an 
entire function g with Dg = f. 

* Work supported by National Science Foundation Grant No. NSF5-G 1010 at the Institute 
for Advanced Study. The author is on leave of absence from the Johns Hopkins University. 

! “Solution of Some Problems of Division. I,” Am. J. Math., pp. 883-908, 1954. 

2 L. Schwartz, Théorie des distributions, Vols. 1 and 2. 

81. Ehrenpreis, ‘‘Analytic Functions and the Fourier Transform of Distributions. I,’’ Ann. 
Math. (to appear). 

4 See “Solution of Some Problems of Division. I,’’ especially the proof of Theorem 8, p. 895. 

> See ibid., especially the proof of Theorem 7, p. 893. 

6 ],, Ehrenpreis, ‘‘Solution of Some Problems of Division. II,’’ Am. J. Math., pp. 286-292, 
1955. See especially the proof of Theorem 1, p. 288. 


DEGENERATE SCALAR INVARIANTS AND THE GROUPS OF MOTION 
OF A RIEMANN SPACE 


By ArTHUR KoMAR 
PALMER PHYSICAL LABORATORY, PRINCETON UNIVERSITY 


Communicated by John A. Wheeler, July 18, 19565 


Introduction: Functionally Independent Scalars and Groups of Motions.—If an 
n-dimensional Riemannian manifold admits n functionally independent scalars— 
constructed from the metric—it clearly cannot admit a continuous group of mo- 
tions. For if A‘ (7 = 1, ...,) are such a set of scalars, we must have §(0A‘/Oz") 
= 0 for the infinitesimal transformations 2‘ ~ x' + &€(a', ..., x”")ét of the group 
(the summation convention is assumed for Greek indices), while by hypothesis the 
rank of the matrix 0At/Ox’ is n. We see, further, that if a space admits groups 
of motions in at most r independent directions, it may admit at most n — r func- 
tionally independent scalar invariants. (Example: Schwarzschild metric: n = 
4,r = 3.) The purpose of this paper is to establish the converse of this result, 
namely, that if a space admits at most n — r functionally independent scalars, 
it will admit groups of motions in r independent directions. 

Notation: Orthogonal Ennuples; Coefficients of Rotation; Intrinsic Components of 
Tensors, Derivatives, and the Riemann Tensor; Integrability Condition; the Ricct 
Ennuple; Asymmetric Connection; Colon-T ype Covariant Differentiation.—We shall 
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denote an orthonormal ennuple by \;,;, where the indices to the left of the bar label 
the vectors and the indices to the right of the bar label the components of the 
vectors. We call a vector “spacelike” if gag\y"\i" > 0 and “timelike” if gag\,)*- 
\,° <0. For spacelike vectors we define A"; = A for timelike vectors we define 
"; = —~A, 


iljs 
; In this notation the ennuple is normalized, so that 
AaliA”; = Qi); (1) 
or, equivalently, 
{ 0 for a+ J, 
h,* = ¢,, = ‘’ for i = j = spacelike, (2) 
—] 


"* = 6,/'; nN ; 
for i =) = timelike. 


U 


r 


ta ta 


We shall denote covariant differentiation with respect to the metric of the space 
by means of a semicolon. Thus 
= OA, J ,a ‘ 
d; ike or* — I iE @’ (3) 
where I';, are the Christoffel symbols. 
The “‘coefficients of rotation” of the ennuple are defined by Eisenhart :! 


Yijk| = Aijaspr7W*ei- (4) 


We see that, as defined, they are scalars, and from equation (1) we deduce that 
Yijk| + Viik| = 0. 

Similarly, using the ennuple, we can define the “intrinsic components” of any 
tensor, which will then be a set of scalars, e.g., 


A jin) = Aap rr Sire, (5) 


and we can define the “intrinsic derivative” 


a) a fe) 6) 
=) (0 
Os' " Oe* 
Using the above notation, we can readily deduce for the intrinsic components of 
the Riemann tensor! 
am OV ijk Ovij1 a @ a@ rd 
R int) = 7H je OT Ve-| (Vast) — Vos) TY -uYast — V -wYast (4) 
Os Os 
and obtain analogous expressions for the intrinsic components of the higher co- 
variant derivatives of the Riemann tensor. 
If f is any sealar function, it satisfies the identity 
o Of Oo Of - of (8) 
= (7 6g: ee 
Os" Os’ Os’ Os" ? om Oe 
otherwise known as ‘‘the condition of integrability of intrinsic derivatives.” 
Throughout the remainder of this paper we shall take \,,; to be the ennuple de- 
fined by the eigenvalues of the Ricci tensor R,;. For this choice we find that not 
all the intrinsic components of the Riemann tensor are independent (apart from the 
evident symmetries due to those of the Riemann tensor), but instead there exist 
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n(n — 1)/2 additional identities: 
2 on a Bal: ’ 
Rix, $s Riek: = Pi\Cik}, 


where p, are the eigenvalues of the Ricci tensor. 
Using this Ricci ennuple, we may now define the asymmetric connection? 
dd" On", 
(10) 


I i = Ye " — = 
4 kj “A ¢ y / + 
Oat ’ Ox" 


We denote covariant differentiation with respect to this connection by a colon, 


thus: 
_ 
am. 


i 


Be ae he (11) 


- 
c 
Ss 


Use of the L';, notation allows one to state conditions of integrability in a covariant 
form. 
By comparing equation (10) with equation (4), we find that 


1 


i i 
x) — 0 ki| = Y -k - jk|: 


In addition, we define the tensor? 


hey! a 
= ye fr, TL, ~ 1%; (13) 


. Ox? oa* 


Summary of Results: Conditions for a Group of Motions; Integrability Condi- 
tions; Rank of Equations; Translation of Conditions into Intrinsic Form; Expression 
of Conditions in Terms of Coefficients of Rotation and Their Intrinsic Derivatives; 
Proof of Existence of Groups of Motions in r Independent Directions.—A sufficient 
condition that a Riemann space admit a group of motions is? 
£.. = 0. (14) 
This condition is also necessary if the Ricci directions are not degenerate. 
The integrability conditions of these equations are? 
ear i 
. L ‘ajk =e 0, 
x U ~ 
EL" ikl = 0), (15) 


ay. co 
g L -ajkl ge 0, 


A necessary and sufficient condition that equation (14) admit a solution is that 
there exist a number qg such that the rank of the matrix of the first ¢ sets of equations 
(15) ben — r => 0 and that the rank of the augmented matrix when the g + 1 set 
of equations are added remain unchanged. When this is the case, equation (14) 

,r) such that the rank of the 


.., «") then acts as one of the generators 


will then admit r independent solutions ,£' (a = 1, ... 
matrix ,é' isr. Each solution ,é'(2!, 
of an r-parameter group of motions. 

To gain some insight into the meaning of these integrability conditions, we will 
re-express them in terms of the coefficients of rotation and their intrinsic deriva- 
tives. Since equations (15) are tensor relations, we can consider them in the 
equivalent intrinsic form: 
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al Dane 0 

a) 0 (16) 
g* L' -jkil m ‘ 0 


If we substitute equation (10) in equation (13) and employ the identity R';;,., = 0 


and equations (7) and (12), we find 
ae (y" a) — Ysa) (17) 


Analogous to equation (11), we have 


OL ix: 


O es | Pye Le fa L* m Diack! i, ha L, al ~ Lim Lika: (18) 
av 


L, kl:m 


Again employing equations (10) and (12), we find 


/ OL, kl a ‘ui / Pe an ] 
4ijklom Pa m + ( Y um eal T mi\ ) 4ajkl ot gh, T m em T m) } 4iakl + 
Ss 


a a a a 
\Y -kml 9 VL CN ake ee gl Like: (19) 


mk jal 


However, from equation (17) and the integrability condition (eq. [8]), we find 


OL 5x1 a ‘a i 
Pa) m r (¥ oO Riel mj\) 44iakl 
8 
Se ei ee 
Os” Os (Yilk| — Yak . m|:— Yo -myji) Os* (Vitk| — Yikij) (20) 
o os 
CYits V iki|)- 


pi Os’ Os” 
If we replace this in equation (19) and employ the first of equations (16), we can 
express the first two equations directly in terms of the coefficients of rotation and 


their intrinsic derivatives: 


ahi 
. Os" (Vis Yijk|) 0) 
s 
ig (OT 
, ds"! Os! (sail — Yaya|) = O.- (21) 
s“' Os 


Similarly, we obtain, for the remaining equations (16), 


= O Go © 
é ds"! Oa” Ds (Yiaji — Yasui) Q, 
pe Re O O 
£ (Yieji — Vix) = Y, (22) 


] 
Os” Os? Os m Os’ 


The coefficients in these equations are scalars. But the number of functionally 


independent scalars by hypothesis is equal ton — r > Consequently, the rank 
of equations (21) and (22) is less than or equal ton — r. Actually, the equality 
holds, for if it did not, equations (21) and (22) would admit more than r linearly 
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independent vectors é", in which case the number of linearly independent scalars A' 
(cf. Introduction) would be less than n — r, contrary to hypothesis. 

Thus the theorem is proved. 

We may add that if the Ricci directions are not degenerate, this r-parameter 
group is the full group for the space. When the Ricci directions are degenerate, 
then the space may admit a larger group of motions, as described by Eisenhart.’ 

It. is of interest to note that the “‘constants of composition” of the group, C4, 


which are defined by 
OnE OnE ee poeeoe 
Cat =f —- i= (23) 
Ox" Ox" “0, em |]: 2.4, 


can be expressed in terms of the particular combination of the coefficients of rotation 
which occur in the integrability equations (21) and (22), via 


1o a a| B i t| 
( ab E - a 3 (¥Y “Bal —~ Y - af) 


$a,8=1,...," 
( ) (24) 
Se oe eS 
'L. P. Eisenhart, Riemannian Geometry (Princeton, N.J.: Princeton University Press, 1926). 
2 L. P. Eisenhart, ‘“Groups of Motions and Ricci Directions,’ Ann. Math., 36, 823, 1935. 


KINEMATICALLY PREFERRED CO-ORDINATE SYSTEMS 
By T. Y. THomas 


GRADUATE INSTITUTE FOR MATHEMATICS AND MECHANICS, INDIANA UNIVERSITY, BLOOMINGTON, 
INDIANA 


Communicated August 9, 1955 


|. Jntreduction.—In a recent communication dealing with the structure of 
the stress-strain relations for a material medium in a state of motion, we made use 
of a rectangular co-ordinate system relative to which the instantaneous rotation of 
the medium vanished at the point and time under consideration.!. Such a co-ordi- 
nate system was referred to as a kinematically preferred co-ordinate system. We 
shall now give a definitive discussion of these co-ordinate systems, and we shall, 
moreover, show how they can be used for the construction of invariants under the 
group G of co-ordinate transformations? 


La = Agg(tts oo b(t), (1) 


relating the co-ordinates x, and %, of moving rectangular systems.’ As indicated, 
the coefficients a,, and 6, in equation (1) depend on the time ft. The condition of 
orthogonality is expressed by either of the following two sets of equations, which 
are in fact completely equivalent, namely, 


Aag(t)Aay(t) se bg; Gaa(t)a,a(t) = Biy- (2) 


The theory of extension, including absolute time and covariant differentiation, 
which is based on the kinematically preferred co-ordinate system, depends on the 











Hi 
| 
' 
i 
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fact that a transformation of the type of equation (1) induces a linear homogeneous 
transformation of the co-ordinates y, and 7, of the preferred systems which are 
associated with the x and # systems, respectively. Hence the kinematically pre- 
ferred systems are completely analogous to the normal co-ordinate systems which 
are introduced in the discussion of Riemann and other generalized spaces. This 
fact makes possible the application of the replacement theorem,’ which provides 
immediate information concerning the structure of sets of invariant differential 
relations, e.g., the stress-strain relations for a correct dynamical theory of a material 
medium. 

2. A Basie Relation.—Defining the velocity components by v, = dx,/dt and 
d, = dé,/dt in the x and & co-ordinate systems, we have, from equation (1), 


Ve = Agg(t)ig + dag(t)%s + 5,(t), (3) 


where the dot is used to denote the time derivative. Differentiating equation 
(3) partially with respect to the co-ordinate J, the resulting equations can be writ- 
ten 

Vas = Age splat eis gp cay (4) 
where the comma is employed in the usual manner to denote partial differentiation. 


If we now interchange the indices a,8 in equation (4), subtract, and then make use 
of the relations obtained by differentiation of equations (2), we can show that 


Gas = 168% ae + QasP Bo; (9) 


where 


1 . : RNs e 
Pap = 2\%ap — YBa); Pap = 2(\Vag — VBa)- 

The simple relation in equation (5) will have a basic application in the following 
discussion. 

3. Kinematically Preferred Systems.—Let us consider, for definiteness, that 
the velocity field is known, and let us represent by x, = g,(t) the trajectory of a mov- 
ing point P of the medium referred to the time ¢ as parameter. Let us also con- 
sider a moving rectangular co-ordinate system y related to the x system by the 
equations 

La = Cash) Ys — Jail). (6) 
This transformation belongs to the group G (see sec. 1) and 1s, moreover, such that the 
y system, at any time t, will have its origin at the moving point P. 

Denote by w, the components of the velocity of points in the medium relative 

to the y system, and put 


1 
Vas a 2(Wa a, We a)- 


We now make the following assumption. The quantities ~,3 vanish at the origin of 
the moving y system. To see the consequence of this condition, let us consider rela- 
tions (5) corresponding to the transformation (6), namely, 


Cap = CopP aa + Cac¥ Bo: (7) 
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Referring equation (7) to the origin of the y system, we have 
Cap “a CopPao(l), (8) 


in which the quantities ¢,, can be regarded, as we have indicated, as functions of 
the time ¢. Assuming the continuity of these functions, it follows that the coeffi- 
cients c,(t) will be determined as a solution of equation (8) when we assign their 
values at some specified time to. We now impose the condition that cas(t) = bas 
att = t. It is to be understood in this connection that we are concerned with the 
behavior of the medium at the particular time fy), which can, of course, be selected 
arbitrarily. This determination of the coefficients c,3(¢) by means of equation (8) 
is such that 


d 


yea _ Caplay + Caplay = CopCay\Pae os Poa) = 0, 
( 


in view of the fact that the quantities ¢,, are skew-symmetric. Hence the ex- 
pression C,Ca, is independent of the time ¢, and, since this expression has the value 
5g, at t = ft, by the above assumption, it follows that the required orthogonality 
condition is satisfied. Moreover, by combining equations (7) and (8), we see that 
the quantities Y,, vanish at the origin of the moving y system. 

The above determination of the transformation (6) depends essentially on the 
velocity field within the medium. We shall refer to any of the co-ordinate systems 
y which are given by equation (6) as a kinematically preferred co-ordinate system. 
In particular, the kinematically preferred system corresponding to the above value 
of the time ¢ = éo will be called the initial system, and the time fy will be called the 
initial time. 

For use in the following discussion, let us here observe that at the initial time fp 
we have 


Cap = OaB; Cap _ Pas; Cap _ (dbag dt) + PacP ap: (9) 


The first of these relations is the result of the above assumption. The second fol- 
lows by evaluation of equation (8) at ¢ = t. The third is obtained by differenti- 
ating equation (8) and evaluating at ¢ = ¢. Higher time derivatives of the c,(¢) 
can readily be determined by a continuation of this process. 

4. Transformations of Kinematically Preferred Systems.—Let us now make a 
transformation (1) of the co-ordinates x,, and let us then consider the kinematically 
preferred co-ordinate system 7 which is associated with the system £ in the same 
way that the system y of section 3 is associated with the x co-ordinate system. We 
suppose in this connection that the y and 7 systems have their origins at the same 
moving point P and also that the initial y and 7 systems are determined for the 
same value of the time ¢ = to, i.e., the initial time is the same for the two systems. 

The relations between the y, and §, co-ordinates will now be deduced. We ob- 
serve immediately that the transformation y <—— ¥ must have the form 


Ya = Capt) Ve, (10) 


since the y and ¥ systems have the moving point P as their common origin. From 
equations (5) we see that 
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€ ap ag CopV ae + Cas¥ Bo 


for the transformation equation (10). But ¥,g and Y,, vanish at the origin of the 
y and ¥ systems. Hence é,, = 0, i.e., the coefficients e,3 in equation (10) must be 
constant. To determine these constants, we consider the relations 

Oa Wa OX, OF, 


= 11 
OV, OX, OF, OVs a2) 


Cap = 


with reference to the initial y and 7 systems. But at the common origin of these 
systems equation (11) reduces to egg = @,3(t0). We have now proved the following 
result. 


When the co-ordinates x, are subjected to a transformation (1), the kinematically 


a 


preferred system y associated with the x system undergoes a linear homogenous trans- 
formation 


Ya = asl to )\Vg ( l 2 ) 


in which the coefficients a,, are constant and equal to the corresponding coefficients 
Ga3(t) in equation (1) at the initial time ty. After the proper evaluations have been 
made, the designation ¢ can be used in place of é, since this time is arbitrary. In the 
following sections we shall show how the kinematically preferred systems can be 
used for the construction of tensors under the transformation group G. 

5. Absolute Time Derivatives.—The discussion in this and the following sections 
will be based on a representative tensor o having components o,,(2, ¢) relative to 
the x co-ordinate system. In addition, certain formulas will be derived which are 
connected with the velocity v, but these will be of a special nature, owing to the 
relationship between the kinematically preferred system and the velocity field in 
the medium. 

Consider the transformation 


T ap a T pvyey_, (13) 


relating the components 7, and 7,, of the tensor o in the y and 7 systems, respec- 


D" oa = (==) . D"b a8 ~ (=) (14) 
De \ a /” Di = \ a /” 


in which the right-hand members are considered to be evaluated at the origin of 
the initial y and 7 systems. By these equations the quantities D"o,,/Dt" and 
D"é,3/ Dt" are defined at the point P and time fy relative to the x and Z systems, re- 
spectively. Now, since the a’s in equation (13) are constant, it follows, by total 
time differentiation of these equations and evaluation at the origin of the initial 


tively. Let us put 


y and ¥ systems, that 
t Pe 


= a 
Dt" Dt" 


pal yp, 

in which the a,, can be regarded as the partial derivatives 0x,/O#, at the time fo. 
Hence the quantities D"o,g/Dt" defined by equations (14) are the components of a 
tensor under the group G. We shall call this tensor the n-th absolute time derivative 
of the tensor o. 
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To determine the formulas for these time derivatives, we have merely to differen- 

tiate the relations 
Tap = TyCpal rp (15) 


and evaluate att = tf. Thus we find 


Deas do ag 


Dt 3 dt + TupP uc + F apP up) (16) 
and 
D?a., dees do,s AG ay ddbva 
ee f 9 f » 
Dt? eo ae ee ee 
do,s 9 
2 T ay dt + “Oye wa? vB + (FpsPvaF auP v8) Puvs 


when use is made of equations (9). Formulas for higher absolute time derivatives 
are obtained in a similar manner by repeated differentiation of equation (15) and 
evaluation at the origin of the initial y system. 

The process of absolute differentiation can also be applied to certain quantities 
of nontensorial character under the group G, provided that the tensor law of trans- 
formation is valid for the linear transformations (12) relating the kinematically 
preferred systems. For example, we have 


Wa = Beas} Dan = BelbeiWec: Vas = Ugh pVor (17) 


va 


under the transformations (12). Hence, if we define quantities D"v,/Dt", 
D"v,.g/ Dt", and D"$,/Dt", corresponding to the definition (14) of the absolute 
time derivative of the tensor ¢, these quantities will be the components of time 
derivatives of tensorial character under the group G. However, owing to the fact 
that w, andy,, vanish at the origin of the moving y system, it follows that 


Dv, D"¢, 
= 0; Pap 


= UV; =0,n = 1,2,3,... (18) 
Dt" Dt" 


The formulas for the derivatives D"v, 3/ Dt" can be obtained by differentiating the 
relations 


Wap ae Ve, eal rp om Coal apy 


by which the quantities v,., are transformed from the x to the y system, and then 
evaluating at the origin of the initial y system. Thus, using equations (9), we 
find 


Dvag Wag Ibapg 


; 19 
Dt dt dt TF Vepboa + Yaebes pink 


6. Co-ordinate Extensions —Let us now differentiate relations (13) repeatedly 
with respect to the co-ordinates ¥ and then evaluate the resulting equations at the 
common origin of the initial preferred systems. We thus show that the quantities 
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Tas/7...6 defined by the equations 


O"rag 
Tap ree - (20) 
Oa.» s/o 


are the components of a tensor, which we may call the n-th co-ordinate extension 
of the tensor o, under the group G. However, the formula for the components 
Tos/,...6 Of this extension, which is obtained by differentiation of equation (15) 
and evaluation at the origin of the initial y system, is merely 


TaB/y...8 ~ TFaB,y... 5) 


and it is of course immediately obvious that the partial derivatives oag,...5 are 
the components of a tensor under the group G. 

Less obvious results can sometimes be obtained if we apply the process of co- 
ordinate extension to a quantity of nontensorial character, which we can do, pro- 
vided that the components of the quantity transform by the tensor law under the 
linear transformation (12). From the first set of equations (17) we see that the 
velocity conponents satisfy this condition. Hence quantities vg/g,.,., can be 
defined by equations of the type of equation (20), and these quantities will be the 
components of a tensor under the group G. To obtain the formula for the com- 
ponents ¥,/g... ,, We must differentiate the relations 


WwW, = C 


a va 


fo > Ceabedle 7 Cus (21) 


with respect to the variables y, and evaluate at the origin of the initial preferred 
system. For the case of first extension or covariant differentiation, this gives 


1 *pe 
Va . Va.B it Pas ” 2(Vag +} Vg, a)- (22) 


It follows that the quantities appearing on the right-hand side of equation (22) are 
the components of a tensor under the group G. This tensor is usually called the 
rate-of-strain tensor. We shall employ the symbol €,, to denote its components in 
the following discussion. 

Now that we have before us the fact that the quantities e,, are the components 
of a tensor, we can simplify the derivation of the formulas for the time derivatives 
D"v,.3/ Dt" considered in section 5. Thus, from the identity 


Va. = €as SE Pas (23) 
and the second set of relations (18), we can immediately infer 
D"vag D" € as 


= Fp 1 28..0 24) 
Dt" Dt" ; 


Since the e,, are of tensor character, the formulas for the time derivatives in the 
right-hand members of equation (24) are similar to those in section 5 for the tensor 
o. Hence, taking n = 1, we have 


Dogg de 


ag or 
Dt dt F fupPua F Cau Pups 7) 


in place of equation (19). 
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Finally, we see from equation (21) that 


l 


’ = 2 
@/B... a@B... ¥) 


provided that there are two or more indices in the set 8... y. However, it is im- 
mediately evident from equation (3) that these partial derivatives v3... , enjoy 
the tensor transformation under the group G. 

7. Mixed Time and Co-ordinate Extension.—Quite generally, the quantities 
Ta3/y......¢ defined in the x co-ordinate system by the equations 


O*r. 
Vep/y ..6%i..t = ( e ) (26) 
OY yess Oy 30 55 ee 0 


where the derivatives in the right-hand member are evaluated at the origin of the 
initial y system, will be the components of a tensor under the group G. When 
partial differentiation with respect to the time is not involved, formula (26) reduces 
to the definition of the co-ordinate extension considered in section 6. 

In particular, equations (26) define the partial time extension with components 
Tas/t...u The relation between the components of this extension and the com- 
ponents of the absolute time derivative in section 5 can be derived from the equa- 
tions 


AT ag at OT as 


+ Tap,y WV, 
dt ol mn 


and the equations which result from them by total differentiation with respect to 
the time. Evaluating these equations at the origin of the initial y system, it is seen 
from the first set of relations (18) and the fact that w, = 0 at the origin that 
Dogs 
Di 


= Tag; Z — = Cag/us ss (27) 
De? 
In other words, the absolute time derivatives of a tensor are identical with its correspond- 
ing partial time extensions. 
By transforming equation (23) to the kinematically preferred system y, differ- 
entiating totally with respect to the time, and evaluating at the origin of the initial 
system, we can deduce that 


Dvas Deas Dv... Dea (28 
= Catan = - " = 0, = rm 20) 
Dt Ps oe De wet 
In particular, from equation (25) and the first of relations (28) we have 
déag } 
Va Bt a5 It + €usP ua + €auP up: (29) 
( 


It may be of some interest to observe that, in general, repeated covariant dif- 
ferentiation is not commutative when both co-ordinate and time variables are in- 
volved. Thus we see immediately that v,,,, equals zero, since v,,, is equal to zero. 
However, the quantities v,/,, are equal to the components of the absolute time 
derivative of the rate-of-strain tensor. 
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8. The Replacement Theorem.—As an illustration of the replacement theorem,® 
consider a scalar or tensor invariant whose components are functions 


’ O43 Wag wv, We 
F\ aaa} 7 a Pas ; (30) 
ot Ox, Ol Otz 


of the components of the tensor o, the velocity v, and their first partial derivatives 
with respect to the co-ordinates and the time. The replacement theorem then 
states that the components / of this invariant can be expressed by replacing the 
arguments in the above functions (30) as follows: 


00 a8 Deas OG ag 
o ——'e - ae = ; = @ 3 
ap ap» ap ’ apy; 
; of Dt ox, Ps 
Ove Oe 
ve -> 0); — 0; ie Va oe €ap- 
ot Org 


This result is obtained immediately by transforming the components (30) to kine- 
matically preferred co-ordinates and evaluating at the origin of the initial system. 
Hence the components F can be expressed by functions of the type® 


H (e483 Tap t; OC ap,y> €ap)- (31) 


This procedure can evidently be extended to invariants whose components in- 
volve derivatives of any order. Thus, if the components (30) are assumed to depend 
also on second co-ordinate and time derivatives of the v,, the following substitutions 
for these derivatives are to be made, namely, 


OV Ove re Deas Ova 


—+ 0): “ae _ 
Dt Ox 02, 


9 =? 0/9. = l 
or Ox Ot ' 


a,By* 


In this case the functions H will involve as arguments the quantities De,3/Dt and 
'4,8y, in addition to the arguments appearing in (31). 
Remark: Consider the following two sets of stress-strain relations: 


U 


do Oa 


ap 


= F (oc, €); = F .a(e, €), (32) 


dt ot 
where the F’s are the components of a tensor invariant of the stress tensor o and 
the rate-of-strain tensor «. If we adopt the point of view which we have recently 
proposed, namely, that such relations are valid only at the origin of a kinematically 
preferred system where the rotation of the medium vanishes, then it follows from 
the first of relations (27) that the invariant formulation of each set of equations 
(32) is the same and is given by! 


Dogg 


Det = F.,(¢, €). (33) 


Hence either set of equations (32), when considered apart from the origin of a kine- 
matically preferred system, can be regarded only as an approximation to the dy- 
namically correct relations (33). 
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'T. Y. Thomas, “On the Structure of the Stress-Strain Relations,’’ these PRocEEDINGs, 41, 
716-720, 1955. See also C. Truesdell, ‘“The Simplest Rate Theory of Pure Elasticity,’’ Communs. 
Pure and Appl. Math., 8, 123-132 1955; and W. Noll, “On the Continuity of the Solid and Fluid 
States,” J. Rat. Mech. Anal., 4, 3-81, 1955. 

2 Since rectangular co-ordinate systems will be used exclusively, there will be no distinction be- 
tween covariant and contravariant indices. Hence indices will be placed in the lower or covariant 
position, including the index on the symbol for the co-ordinates which have therefore been denoted 
by tq and £, in relations (1). 

* Prepared for the Applied Mathematics Branch, Naval Research Laboratory, Washington, 
D.C. This paper was written while visiting the University of California, Numerical Analysis 
Research, Los Angeles, California. 

4 The completely analogous definition of extension based on normal co-ordinates was introduced 
in the paper by O. Veblen and T. Y. Thomas, ‘“‘The Geometry of Paths,’’ Trans. Am. Math. Soc., 
25, 551-608, 1923. For further discussion of this subject see O. Veblen, Jnvariants of Quadratic 
Differential Forms (London: Cambridge University Press, 1927); also T. Y. Thomas, The Dif- 
ferential Invariants of Generalized Spaces, (London: Cambridge University Press, 1934). 

5 The replacement theorem first appeared in an article by the present writer entitled ‘“‘A Pro- 
jectory Theory of Affinely Connected Manifolds,’’ Math. Z., 25, 723-733, 1926. For references 
to earlier work employing this theorem see T. Y. Thomas, The Differential Invariants of Generalized 
Spaces, p. 111. 

6 The procedure involved in the proof of the replacement theorem shows that the functions F’ 
cannot depend explicitly on the co-ordinates, since these variables cannot enter as arguments of 
the functions H. 


RELATION OF THE THEORY OF CERTAIN TRINOMIAL EQUATIONS IN 
A FINITE FIELD TO FERMAT’S LAST THEOREM 
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1. In several previous papers" * * 4 the writer brought out a relation between 
a certain theory of trinomial congruences and the second case of Fermat’s Last 
Theorem. In this work the aim was to employ powerful tools from each of two 
deep theories in number theory and to unify the methods employed. Of course, 
it would be nice to prove the Fermat theorem, and by any means possible; but 
failing this, any by-products derived from such a failure are more likely to be novel 
and of value if we employ initially what appear to be some of the important results 
in algebraic number theory than if we stay with elementary results and method. 

The depth of the first theory mentioned above is well known; it includes the theory 
of the laws of reciprocity, class numbers, and class fields in the theory of alge- 
braic fields, and stemmed originally from the consideration of the conditional bi- 
nomial congruence 


n 


"= (mod py), 


where p is a prime ideal in an algebraic field AK and @ an integer in K, not an nth 
power in K. This led to laws of reciprocity for special values of n (generalizing 
quadratic reciprocity) and particularly for /th powers, where / is a regular prime, 





') — _ —_ _ — —_—— _ —- — Es 
i 

1 

i} 
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] , ' >" 
! and the other topics mentioned above followed in its wake. However, the depth 
it of the second topic is not so well known, and I should like lo explain my viewpoint 


here. The latter theory had its beginnings in the work of Lagrange® in his con- 
tributions to cyclotomy, mainly his use of the Lagrange resolvent. Then Gauss’s 






















\ i 
|) work in exponential sums (cyclotomic periods) led, for special values of m, to the 
t} consideration of the congruence 
ax™ + by™"=d (mod p), 
nL 
a's where all the symbols represent rational integers. Then, based on the Lagrange 
ait resolvent 
it 
Yh $ 9 2 ~l —2 2 
if ria, f) = atta’ + fa +...+ 0770" (A) 
i where ¢ = e7'"/" @ = e”'*’”, » = 1 + lm, risa primitive root of p, and I is an odd 
i prime, Cauchy,® Jacobi, and Eisenstein studied the sum 
ik 
bh+(a+b) ind (Gh + 1) 
V(a) = > a y 7 ; (B) 
h 
where A ranges over all the integers 0, 1,..., p — 2, excepting (p — 1)/2, and 
M) 
, gint* == zy (mod p). 
A These numbers (A) and (B) were further studied and extended by Kummer, 


Stickelberger, H. H. Mitchell, and, particularly since 1920, many others, and the 
literature is now vast. Many applications have been worked out, mainly to the 
theory of conditional equations in finite fields. And, in particular, we shall point 
out some work that has been done which perhaps illustrates best how far-reaching 


our topic is. Fueter? used the complex multiplication of elliptic functions to gener- 
alize the Lagrange resolvent (A) and the Cauchy-Jacobi sum (B) and obtained ex- 
ponential sums involving the use of a quadratic field and then the field obtained by 


- <= 


adjoining a cyclotomic field to the former. This resulted in proofs of certain impor- 
tant known laws of reciprocity involving this extended field, so that here, in a sense, 
the two theories we mentioned at the beginning of our paper merge into one. It 
is this development and several others that led me to use the term ‘‘depth”’ in con- 
nection with our second general topic. The theory of the symbol (7, 7),; for ¢ and 1 
arbitrary (as defined just below our relation [7] for / prime), which is used so much 
in the present paper and in a number of our previous ones, has yielded applica- 
tions to the theory of quadratic forms, the theories of congruences in algebraic 
fields, finite geometries, algebraic varieties in geometry, difference sets, and Dio- 
phantine analysis, as well as to the theory of cyclotomic fields. 


—— 


= 


ee 


2. Consider 
a + y! +2) =0, (1) 


| being an odd prime, xyz # 0 (mod J), the so-called first case of equation (1). 
Let ¢ = e”*' and let a be a primary integer in the field K(¢) or K with 
(a, 2(¢ — 1)) = 1; then, from equation (1), 


(x + gy) = ai, (2) ° 


0,1,...,2—1,a,;an idealin K. Then 
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Qa 
(. 4 <) : 


i= 0,1,...,/— 1, and, since a is primary, we have by the law of reciprocity for 


r a sot 
( re) = 1. (3) 
a 


We now note also that if (a) = 6”, with 6 an ideal eX with (h, 1) = 1, then equation 


r+ ¢ 
(: “ ’) = ], (3a) 


Note the difference in our discussion? of the second case of the Fermat problem. 


Here we found that if y == 0 (mod J), 


r+ ¢' 
( ey) = 1, (4) 
Dd 


. . . - h . . . . h . 
where d is any ideal in K such that d" is a principal ideal such that >” = (a), with 
ain K and (h,/) = 1, whereas a was limited to be primary, as we explained above for 


“JI 
~] 
bo 


Ith powers in K, 


(3) gives 


the first case of equation (1). 
If, in equation (3a), we set 6 = p, a prime ideal in K, then it follows as elsewhere 
(relation (4) of an earlier paper*), with p; « K, that 


a+ ty = p; (mod p), (5) 
i= 0,1,...,1— 1, and p divides (p) with p an odd prime >-zyz. 
Consider the equation in the finite field F(p") 
gt 4 git tabs 1: (6) 
p" = 1+ cl = 1 + el*, since yp’ is primary; 7 is in the set 0, 1,...,¢ — 1;/ is in 
the set 0, 1,..., 1 — 1; 8 is in the set 0, 1,...,12 — 1;¢isin the set 0,1,..., 


ec — 1; and g isa multiplicative generator of the nonzero elements of F(p"). Now 
equation (5) may be written in the form (6), if we first write (5), using s in place of 


2, as 
Y ue i a fae ~ 
— +27 ‘pj =1 (mod p) (7) 


and set ¢ =i‘ (mod p), —y/x =gi', 2 | =gv’, p: =’, where g; is a primitive root 
of pin K. The residue classes modulo p form a finite field of order p", with p” 
also the norm of pin K. Hence g corresponds to g;, and x and y have corresponding 
integral elements in the F(p"). Hence relation (7) becomes relation (6) in F(p"). 
In equation (6), represent the number of solutions s, ¢ as (7, j).,; then it follows from 
relation (7) that (7, 7) = J, an example? of case A of equation (6). 

3. Our work so far in examining equation (1) has made use of powerful tools in 
the first general theory mentioned in our introduction. We now apply some re- 
sults which belong to the second topic mentioned at the beginning. In a previous 
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paper’ we obtained the relation, since p” | + el? and m is replaced by the odd 
prime 1, 


ja] — 
>> (b, 2)... (6b —a,t—ad)y = bf (d,7)(d — a,j — b),, (8) 
l=0 7=0 
d 
with (7, 7). = (2, j)ccs d ranges over the set of distinct integers in the set 0, 1,..., 


c — 1 of the form d,; + kl, with d; some integer and with 0 < d; < l. 
If we take a = 0 in equation (8) and d # 0 (mod J) and use 6 in place of 7 in g' = 

y/x, then we find’ that we have 
(dy) .(d, 7 — b), 0, (9) 


for any d # O (mod 2), b # O (mod c), which we may use if we assume y/zx ¥ 1 
| in F(p"), which we may do when p is sufficiently large; ¢ is any integer in the set 


ee Sa c — 1. In another paper‘ the writer gave the result that if p" = 1 + el, 
i with (c, 1) = 1,v = el, 
m—1 c-1 
> (b, Dea (b — a,j — adder > (di, jeld — a,j — b)., (10) 
1=0 j re 
where d ranges over d,, d; + 1, d; + 2/,..., d, + (e — 1)l, a and d, are any in- 


tegers such that we do not have both the relations d; = 0 (mod J) and a = 0 (mod 
c) holding. Now, if we examine the argument we gave to prove relation (10), we 
find that it will also hold if we assume that the only restriction on ¢ and lis p" = 1 + 











\\ cl. If we use it in this more general form, we may employ it to obtain relation (8). 
iy For, if we take relation (10) and substitute in turn d;, di, + ¢,d; + 2¢e,...,d;+ 
| ae . ° . . P e 
\ (1 — 1)e, add the resulting equations, and employ the formula (relation (13) of an 
ial earlier paper’), 

HI 

t) H m | 

> (d, + re, ] — b),, (di, 7 — b),, (11) 
{ 0 


Waa 


we have equation (8). 


In a previous paper? we found a relation of the type (7) applying to the second 
case of Fermat’s Last Theorem, but as a relation in the finite field F(p), which 
corresponds to a congruence modulo p, as then we had in mind the application of 
rapid computing machines to some of our related problems, and this was relevant 
is since g in this case is a rational integer and p = | + el, that is, the special case 
when n = 1. 


a re a 





Relation (1) shows that in relation (7) we may replace y/x by x/y, 2/2, x/z, y/z, 
z/y, and obtain relations similar to (6) involving equations in a finite field F(p"), 
and in the same way we shall find that for these different pairs there exist an 7 

and 7 in relation (6) such that (b, h).; = land (6, t)., = Ofor # hmodulol. Hence 
we may employ equation (10) and find, as elsewhere (relation(11) of another paper‘), 


(d, 2),.(d, 7 — b). = 0 (11) 


(where, if we replace l by v in eq. [6], then (j, 2),- is the number of solutions ¢, s in 
eq. [6]), which we have already noted holds (in the present paper) for any c and m 
such that p” = 1 + em, and we here set m = /, an odd prime. 
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If we use some previous results? (p. 249) with relations (6) and (7) of the present 
paper, we may state 

THEOREM I (FIRST CASE). Assume the following: 

1. Relation (1) holds with xyz = 0 (mod 1), with | a given odd prime. 

2. In the algebraic field defined by e?'*’' = ¢(say, K), let there be a prime ideal » 
such that p" = (6), with h # 0 (mod 1) and @ a primary integer in K; also, (p) = 0 
(mod »), where p ts an odd rational prime >xyz. 

Then there exists a finite field F(p") of order p" = 1 + cl, ¢ = el; and the norm of 
p is p" for any p satisfying the last assumption and such that relations (9) hold, where 
dis any integer of the form d; + kl,k = 0,1,...,¢: — 1, with d; # 0 (mod 1); 7 is 
any element of the set 0, 1,...,¢ — 1; and tf g is a multiplicative generator of the 
nonzero elements of F(p"), then g” takes on any one of the values 

we Ee He Sy tet, ee oe in F(p"). 
x y Z x z y 


The symbol (h, k), for any h and k, used in formula (9), ts defined just below our rela- 


tion (8). 
THEOREM II (SECOND CASE). Assume the following: 
1. Relation (1) holds with y = 0 (mod 1) (xz, l) = 1, y ¥ 0, la given odd prime. 
2. Suppose that p is an ideal prime dividing (p), where p is a rational odd prime 


>axyz in K(¢?"*) and yp" is a principal ideal in K with h # 0 (mod 1), with p" the norm 
of p. 

Then there exists a finite field of order p" for every p satisfying the last assumption 
such that relations (11) hold in the finite field F(p"), where g’ has either of the values 
—y/x, —y/z, in F(p"). 


' Erna H. Pearson and H. 8. Vandiver, “On a New Problem concerning Trinomial Congru- 
ences Involving Rational Integers,”’ these PROCEEDINGS, 39, 1278-1285, 1953. 

2H. S. Vandiver, ‘‘New Types of Trinomial Congruence Criteria Applying to Fermat’s Last 
Theorem,” these PROCEEDINGs, 40, 248-252, 1954. 

3H. 8. Vandiver, ‘On Trinomial Equations in a Finite Field,” ibid., pp. 1008-1010. 

*H.S. Vandiver, “On the Properties of Certain Trinomial Equations in a Finite Field,”’ cbid., 
41, 651-653, 1955. 

5 (Euvres, 2, 375, 1868. 

6 It seems that the earliest published paper in which the number (B) appears is an article by 
Cauchy: Bull. sci. math. astr. phys. chim. (ed. Ferussac), 12, 205, 1829; Mem. Acad. Sci. Inst. 
France (Ser. 2), 17, 1840; In his Quores, Ser. 1, 3, 112, he gave a wide generalization of the 
number (B), which was only rediscovered, apparently, several years ago. 

7 Nachr. Ges. Wiss. Gottingen, Math.-physik. Kl., 1-11, 427-455, 1927. A historical note ap- 
pears in order here. I spent several months in Zurich in 1928 and visited with Professor Fueter 
a number of times. He talked to me a great deal about the contents of the above papers and how 
these ideas generalized a large part of the material given by Hilbert in his Jahresber. Deut. Math.- 
Vereinitgung “Report on Algebraic Numbers’”’ ((2uvres, 1, 195-248). There is a connection here 
with the remarks made by Hilbert at the end of the third chapter (p. 194), which treated the 
theory of quadratic fields, and at the end of this chapter he refers to a ““Héheren Theorie’’ of 
quadratic fields obtained by adjoining other fields to the quadratic. He then says: ‘Die 
theorie der zu einem imaginiiren quadratischen K6érper gehérigen Klassenkérper sowie der dazu 
gehorigen relativ-Abelschen Kérper erfordert jedoch zu ihrem Aufbau die Methode der komplexen 
Multiplikation der elliptischen Funktionen, und dies ist ein Gegenstand, welcher von der Auf- 
nahme in diesen Bericht ausgeschlossen werden miisste.”’ 

Fueter told me that Hilbert had been to see him in Zurich and had remarked that the reason 
he had omitted the use of the complex multiplication of elliptic functions in the report was that 
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this subject had not been sufficiently developed, as far as possible applications to number theory 
were concerned, at the time he wrote his report. Hilbert also said that he was very glad to hear 
that Fueter had fashioned the necessary tools from the theory of elliptic functions to yield an ad- 
vanced theory of quadratic fields. 

We hope to be able to discuss in a later paper the methods of Fueter in their relation to the class 
field theory as it has been developed up to the present. As to the work of Fueter’s predecessors in 
this line, cf. H. Hasse, Jahresber. Deut. Math.-, Vereinigung, 35, 39-44, 1926. 


ON CYCLOTOMIC RELATIONS AND TRINOMIAL EQUATIONS IN A 
FINITE FIELD 
By H.S. VANDIVER 
DEPARTMENT OF MATHEMATICS AND ASTRONOMY, UNIVERSITY OF TEXAS 


Communicated July 22, 1955 


Let p be an odd prime and F|[p”| denote a finite field of order p” and consider 
the solutions s, ¢ of 


git + git S. 1. (1) 
where g is a multiplicative generator of the nonzero elements of F[p”|; 7 is a fixed 
integer in the set 0, 1, ..., ¢ — 1; 7 is a fixed integer in the set 0,1, ...,/ — 1; 
p" — 1 = cl; sis in the range 0, 1, ...,/ — 1; ¢ 1s in the range 0, 1, ...,¢ — 1. 


Let the number of sets of solutions s, ¢, of the type above be (7, 7).:. Then (noting 
relation [2] of the present paper) the maximum value of the (2, 7). is life > lL. 
If there exists a 7 for a fixed 7 such that (7, 7).: = 1, we shall call this “case A”’ of (1). 

In several previous papers various consequences of this assumption were ob- 
tained,'~‘ and the relation of some of them to Fermat’s Last Theorem was brought 
out.2 Other results of this type will be derived in the present paper. 

To obtain the desired results concerning (1), it will be convenient to first 
write equation (1) in another form, and we shall also set 1 = m, with / an odd prime, 
and assume that c # 0 (mod /). Let [7, 7] denote the number of solutions of the 
equation 


P ioe ote | # git® 
in a finite field F(p") of order p” = le + 1, s’ being in the range 0, 1, ...,/ — 1, and 
t’ in the range 0, 1, ...,¢ — 1, for givenz andj. It is known that 


- 
jo if a gf - (mod ec), 


0 ny c 
l if 1=>5 (mod ec). 


— 


Evidently [7,7] = [¢ + ca, j + ly]. Let g be a primitive root in the finite field of 
order p” and g* = y, with x an integer and y an element of the F(p"). We write 


ind y=2z (mod p” — 1). 


It is known that 
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pd (a, Am tre Tt ba > Mita “eda (ght 4 1) (2) 
ti,J hy 
where a, b, and ¢ are any integers in the set 0, 1, ,' — 1, and 
Sailr Pag (1) mee) lee pees. 
>> [%2, 7218 e's js(t) a ~ tind (g +1) (2a) 
12,J ki 
and where 7; and 7 range over 0, 1, ...,¢ — 1; j; and je over 0,1, ...,1 — 1; and 
h, over 0, 1, ..., p" — 1, excepting (p" — 1)/2. Also, k; ranges over this latter set, 
and 
—; Z ‘ é * rhal z h+( jo a ‘ 
Bk >. (1, jr] (te, Je gal i ¢? le) ee = rt aa (3) 
ti, Ji 42,9 hii 
where A = a(h, — k;) and B = bind (g" + 1) + ¢ (ind (g”" + 1) — ind (g" + 1)). 
Set lo = iT v, qe —_ yn = d, 
4 | 
4 h 4 h Y - k 
g q ae out 
g + 1 
Conversely, these give 
eae h Sel i+¥ 
g ea nh g v k—h gy 7 k—h 
i paket ‘| 1 — qd g — |] 


so that if we substitute the (A, /) sets for the (h;, k;) sets, then we obtain all possible 
= 0, ork, = h. Also, k = h givesh, = ky. Also, the 


(hy, ky) sets unless k = 0, h 
0 gives h = 0 and conversely, so that these exceptions are included 


condition k = 
under the exception inh = k. Then (3) becomes, if we insert on the right the values 
d over 


ky, 1 = di, 7 conde 


hi, letting v range over 0, 1, ,ce- 1; 


such that h; = 


O23, /— 1] 
>> 7. AiG ta j vi d|g wes ge ee +z gua a yi a (4) 
i,7,v,d h,k h 
where C = b(ind (g~" — 1) — ind (g*~" — 1)) — tk and D = bind (g” + 1). Write 
k — hfor s’, and set h in lieu of (—h); this gives, if s’ ¥ h, 
Dy. st +5 + A sP-" = Fee + 3, (5) 
h,s hi 


1,.jv,a 


b(ind (g" — 1) — ind (g* — 1)) — t(s’ — h) and F = bind (g” + 1) 


where FE = 
We now collect the terms involving distinct 8°’s on the right, so that the above may 


be written in the form, if s, + H,, 
( l } 
tO Se Be (5a) 
e=0 Hs 
where G = b(ind (g”* — 1) — ind (g* — 1)) — t(s, — H.) and where the H,’s range 
In view of this, we see that (5) 


over all the integers <c/l of the form e (mod ¢). 


may be written in the form 
A, 180 ~ 9 + A, 28°"? + 2. + Ae = 0, 

= (0,1,...,¢—1. Hence, in the complex field, the equation 

A,-ot "+... + Ap = 0 


which is true for a = 








en 


—— 
- 


—= 








4 
i 
i 
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has ¢ distinct roots; consequently, 


A, =0; f=6@,1,...,.¢e—1. (5b) 
We now assume that in 
ca eg  shebge (6) 
for a fixed m # O (mod c) as s; ranges from 1 tol — 1, we get a fixed r on the right, 
so that, fO <r <i, 
ind (q”* — ]) r (mod l). (6a) 
for m fixed. Now set ¢ m in (5a); then, from (5) and (5b), A,, = 0 gives, if we 
set 7 = H,.. 
> [i 7]fi — mj + dle” aD ww (7) 
i,j, d H 
where / b(ind (yg! — 1) — ind (g* — 1 i(s — H), where H ranges over each of 
the integers in the set 0, 1, ,¢ — | of the form m (mod c), s # H. As in the 


left-hand member, we can substitute 7 + ¢ in lieu of 7; then we may write, if « = 
(p” — 1)/2, 


~ & AG —mj+aer-* = de, (8) 
i,j,d H 
where J b(ind (g” — 1) ind (g° — 1)) — &(s — H) and where (7, 7) is the 


number of solutions in the F'(p”) of 
g +g" l. (9) 


This is the same as the number of solutions of (6), with m replacing 7 and r replacing 


j. Also, this is equivalent to (1), with s», é replacing s, ¢, of (1), and we write 


(7, 7) as an abbreviation for (7, 7).:, since we keep ¢ and / fixed throughout our 
argument. Now introduce the terms s on the right such that H = s. The number 
of these is 1. 


Assume that ¢ + 0. and use (6a): then we have 


oe ° pe ~_jb—dt | bt d (2° 1) ts ~ —tH 
> GG — mj + dj + | pay Be 2 1 
ij.d H 
where s now ranges over all the integers in the set 1, 2, , le — 1. If mis prime 
to l, then the integers </c of the form m modulo ¢ are congruent in some order to 0, 
1,2, ...,l - 1 (modulol). Hence, since m # 0 (mod ce), 
/ 


er a Se and l+ } @A)G-— mi t+ d)j¢?-" = 0, (11) 
H i,j,d 


fort #0,0 <b<l. However, if ¢ = 0, then we obtain 


>» @&ADG—mj7 +a” = Fee -mMe- YT - 1. (12) 


Lnre 


For b ¥ 0 we have 


br < .—bds _b 
6 4 ei ee 
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Si 


where s; ranges over 0, 1, ..., le — 1, except (p" — 1)/2, using g’ — 1 = g". This 
gives, for b ¥ 0, 


» (J — mj + dye” = —-"} — I: (13) 


t,J 
and for b = 0 we find from (12) 


> @&|)Gi—mj+d) =Ule-1) -1 


i,j, d 
whence 


(2, )(@ — m,j + d) = Ule — 2) for i= 0b = 0. (14) 
‘, 3,4 


Then, in view of (13) and (14), we have, if 


A,=DLD)G- m,j + d), (15) 


i,d 


by eliminating ¢ and its powers by known methods from these / relations, if r # 0, 


Ao = l(c bea 1), | 
A, = l(c — 1), (16) 
A; =Ie if J # 0, I# es 
If, however, r = 0, we find 
A, = el, J #03; (17) 


Ao = l(c — 2). { 


Now, since [7,7] = (7 + ¢,j), then (la) gives }> (¢ — m,j +d) = land 


d 


> (i, (i — mj +d) = Ui, j) for i zm (mod c), (18) 


d 


> (m, j)(0, 7 + a) = I(m, J) — (m, J). (18a) 
a 


Hence 


fol — (m, 7) for J # O, 


i EA a FE Vc — 1)l — (m, 3) for 7 = 0. (19) 
Now set 
B, = ‘P (2, pa —_ m, J =} d)j¢~”: t ~ 0. (20) 
i,d 


Using (11), we may then write 


1-1 


D Bg” = -l, (21) 
=0 
6 =0,1,...,1— 1. Eliminating the quantities ¢” from these relations, we obtain 
By = —-l, | 


] 9s 
B.=0, jj 05 -) 
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We have 


bai 
> Coa = Ue — 1) 
i=0 
where 
pee 
Cyu= >, (77)(« — m,7 + d) (23) 
1=0 
and 
1-1 
ey C ja¢ ; —l, (23a) 
d=0 
forj = 1, 2, ,/— 1. Elimination of the quantities ¢~’ yields 
Coo 2c. l, Coa = ¢, d = Q. (24) 
In a similar manner we obtain 
C3 = ¢— I, ry © d arbitrary, (25) 
C34 = ¢, ) # 9, 1 #7, r £0, d arbitrary, (26) 
Coo = Cct-— l ate I. Coa =Cct- l, ad F Q). (27) 
Cya = C, ) # 9, d arbitrary. (27a) 


In view of (24), (25), (26), and (27), we may summarize these by the following: 
THEOREM I. Let (7, j) represent the number of solutions in F(p") of equation (1), 
and p" — 1 = el; lis an odd prime with c # 0 (mod l): 


Cja = > (1, 7)(« —m,j7 +d); 


1=0 


also, 
‘a 


ind (g™**' — 1) =r (mod 1), (27b) 


for each s,; with m, ce, and r fixed; thenifO<n<ld<d<lO<j<l, 


Co =c-l r #0; Cos = ¢, d # 0, r #0, 

Ca=z=e- i] for any d with r ¥ 0, 

Cy =, 3 #0, j YT, r #0, for any d; (28) 
Co =ec—-1-—-], r = 0; Cos = € — 1, d # 0, r = 0, 

Cya = 6, 7# 0, r = 0, for any d. 


The assumption in (27) is equivalent to case A in connection with equation (1). 

In a previous paper‘ we found the following results: If case A holds as defined 
in the beginning of the present paper and if c # 0 (mod J), if d is any integer # 0 
(mod /) in the range 0, 1, ...,v — 1, then 


(d. t),.(d. 2 — b).. = 0. 
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Also, 7 is any integer in the set 0, 1, ...,¢ — 1, (6, h)c = U for some h with v = 
p" — 1. Note the differences between this result and Theorem I of the present 
paper. The latter gives criteria depending on the value of h, which is not the case 
for the other criteria. Also, our present Theorem I involves only symbols of the 
type (2, j)em instead of (2, J)rc and (2, J) ce. 

| E, H. Pearson and H. 8. Vandiver, “On a New Problem concerning Trinomial Congruences 
Involving Rational Integers,’”’ these PROCEEDINGS, 39, 1278-1285, 1953. 

2H. S. Vandiver, ‘“New Types of Trinomial Congruence Criteria Applying to Fermat’s Last 
Theorem,”’ ibid., 40, 248-252, 1954. 

3 H.S. Vandiver, “On Trinomial Equations in a Finite Field,’’ ibid., pp. 1008-1010. 

4H. 8S. Vandiver, “On the Properties of Certain Trinomial Equations in a Finite Field,” zbid., 
41, 651-653, 1955. 


ON THE DIVISORS OF THE SECOND FACTOR OF THE CLASS NUMBER 
OF A CYCLOTOMIC FIELD 
By H. 8S. VANDIVER 
DEPARTMENT OF MATHEMATICS AND ASTRONOMY, UNIVERSITY OF TEXAS 
Communicated July 18, 1925 

Let ¢ = e”*/" where J is an odd prime, and K(¢) or K be the cyclotomic field 
generated by ¢. Let A‘ be the class number of the field. It is generally written 
in the form 


AO = hh, (1) 


where hy‘? and h.“” are each positive integers which were determined by Kummer 
in algebraic form.' A key problem in the theory of cyclotomic fields! * is the ques- 


tion of the possible existence of fields AK such that ho =0(modl/). Let k; = ¢" 


ace 2 and e; = k;41/k;. Set 
a! —2ta 
Et) =E= 1 ¢. (2) 
1=0 
with » = (1 — 1)/2. Then it is known* that if ho“ = 0 (mod J), then F,(¢) is 


the [th power of an integer in K for at least one value of a. 
This criterion was shown! to be equivalent to the statement that if h, = 0 (mod 
l), then 


u 
oa /S és pl—1—2a ; 
Q;.= ¢°°" TI GQ? 1) = ay’, (3) 
b= 
where ap « K, with d = 1'~7* + 2'-*2 4 |... + ,'~*4 for some value of a in the 
set 1,2, ...,4— 1. From this it is evident that 
_-d/2 u re pt! —1—2a 1 ; 
§ lt (¢"* — 1) = a, (4) 
i=1 
with a e K, and where };, bo, ... b, is a set of integers included in the set 1, 2, 


1 — 1, having the property that if 6; is in the set, then / — 6; is in the set. 











| 
j 
} 
| 
] 
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The problem as to A‘? = 0 (mod J) enters into all the eriteria® which have been 
given for case II of Fermat’s Last Theorem which have so far yielded any proofs 
of the theorem for special exponents / > 7. In fact, if ho“? = 0 (mod 1) for some 1, 
then all such criteria fail for said /.  Also,® if the Fermat equation holds in case I 
of the theorem for some /, then h.‘” = 0 (mod /) in K, for said 1. 

From (4) we find easily that 


“ 

—?), ail Bi " 

¢ II (¢° — 1)" = ay! (0) 
#=1 


2a 


where B, is the least positive residue of 6; , modulo / and az e K. It was 
noted‘ (see pages 28-29) that ap is either an integer in K or may be expressed as a 
fraction in K whose denominator is of the form (1 — ¢)", h > 0, and whose 
numerator @ is such that (6, 1 — ¢) = 1, with @ an integer in K. This is obviously 
true also for @ in (4) and a in (5). Hence, from (5), we have 
“ 
aQ— M4? Og — 1)® = a, (6) 
+= 1 

and this gives h = 0 because of the unique factorization theorem for the integral 
ideals in K. Thus a» in (5) is an integer eK. 

We shall now describe a direct method for the examination of the properties of 
the left-hand member of (5). It may be written in the form 


Co+ Ce +... + Cat, (7) 


with the C’s rational integers. We shall show how the C’s may be determined 
in any particular case without multiplying out the binomials in (5) explicitly. Set 


P(z) = 2 ?? TI (x — 1)*. (8) 


Nicol and the writer’ gave the recursion formula 


m 


2 H,,~:(Z) » ka,Z’ = (m + 1)H,41(Z), (9) 


h=0 rk=h+1 
where the inner sum is taken over all integral solutions for which r > 0 and n > 
k > 0, and the H’s are defined as follows: Let x and z be indeterminates over the 
complex field, and set 


nN 
F = F(Z, x, n, a, _ I (1 — ar')”, (10) 
i=1 
where a), a, ..., @, are nonnegative integers not all of which are zero, and write 
n 
F= >> H,(Z)z", (10a) 
s=0 
where 
ul 
N= ) UA; 
j=1 


Relation (9) may be used to calculate the coefficients of 
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ll (x** — 1)* (11) 


i=! 


when it is written in explicit polynomial form. Write (10) as 


“ Me 
MT (a — 1) = 1 (—1)84(1 — a?) Fe (11a) 
t=1 i=1 

and then set z = 1, B; = A,in (9). Now in the proof’ of (9) we might have sub- 

stituted z = 1 originally and assumed that x was any nonnegative number <1; 


hence it follows that (9) is true for all values of x, since otherwise we have a poly- 
nomial with integral coefficients with more zeros than its degree. 

Using (10a) with (9) for z = 1, it is clear that we can calculate Ho(1), H,(1), 
H,(1), ete., in turn from (9) by taking m = 0, 1, ..., in turn in (9). This gives 
the coefficients of the polynomial in (11a), and hence the coefficients of (8), where 
u now has any positive integral value. We may then set 2 = ¢, fory = (1 — 1)/2, 
in our polynomial form of (8) and reduce the expression, using 


FD © a ae 


then the left-hand member of (5) may be written in the form (7), and we may 
examine the possibility of the existence of (5) by using its reduced form. We 
have a wide choice for the b’s, so that perhaps it would be best to start with the 
form (3). Also, we might obtain simpler forms by not reducing the exponents in 
(3) or (5), modulo J. Clearly, however, we could not proceed very far, because 
of the length of the computations, unless we employed a rapid computing machine. 

In a previous paper‘ by D. H. Lehmer, Emma Lehmer, and myself, we gave a 
criterion in connection with Fermat’s Last Theorem as follows: 

Let ¢ be any integer such that t* # 1 (mod p), where p is a prime of the form 
p=kl+1<l-—lU. Define Q, by 


pi —2—2a 


u 
Q,= 1" fe? — 1) (12) 


where 


Let 
Ba, . a ee | &.. 


be the only Bernoulli numbers, with indices less than (1 — 1)/2, which are divisible 
by l. Then if 


OF (mod p) (13) 


holds for 7 = 1, 2, ..., s, Fermat’s Last Theorem holds for the exponent /. 
Obviously the expression in (12) is the same as the one in (3), with ¢* replacing 
¢ in the latter. Pursuing the methods of the present paper, it is also clear that 
we cannot expect to find out much about the residue* of Q,,” unless we know a great 
deal about the coefficients in (8) for the case where the 6,’s constitute the set 1, 


2, 3; : B. 
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INHERITANCE AND REGENERATION OF CYTOPLASMIC DAMAGE IN 
PARAMECIUM AURELIA 


By Eart D. Hanson* 


OSBORN ZOOLOGICAL LABORATORY, YALE UNIVERSITY 


Communicated by J. S. Nicholas, July 29, 1956 


The gullet of Paramecium aurelia is a complicated organelle which duplicates itself 
at each fission. This process involves the formation of an anlage by the outpocket- 
ing of the right, dorsal side of the pharynx,' which then differentiates into a new 
gullet in the posterior fission product (opisthe). The gullet in the anterior fission 
product (proter) dedifferentiates partially during formation of the new gullet bud 
and then redifferentiates into a functional gullet. Studies on P. caudatum and P. 
multimicronucleatum* have shown that feeding is necessary for the regeneration of 
missing or damaged structures. No regeneration of gullet injuries was found. This 
has been interpreted as being due either to the inability of Paramecium to form feed- 
ing organelles de novo? or to gullet damage precluding feeding which prevented re- 
generation.* The term “gullet” includes mouth, pharynx, esophagus, and the food- 
vacuole-forming region. These latter four terms‘ refer to parts of the gullet ar- 
ranged in an essentially anterior-posterior alignment, the mouth being the most an- 
terior structure. 

Reported here are the results of irradiating the gullet, or its parts only, by means 
of a fine beam of ultraviolet light. _Homopolar doublets of P. aurelia were studied. 
The use of double animals, a natural parabiosis, permitted damage to one gullet 
while allowing the animal to feed from the other, uninjured gullet. The gullet can 
be selectively damaged, and the damage can be restituted. Analysis of this re- 
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generation shows that formation of a new gullet does not depend on the normality of 
a pre-existing gullet. 

Methods.—The ultraviolet microbeam is, in principle, the same as that developed 
by Uretz, Bloom, and Zirkle. A Leitz Ortholux microscope with a special nose- 
piece carrying a side arm holding the primary aperture, a beam-splitter, and the 
Bausch and Lomb UV reflecting objective (design V, 2.8 mm., 0.72 N.A. 53X) was 
used for irradiation. The microspot is 5 «square. The energy of the microbeam has 
not been measured; its source is a Hanovia high-pressure mercury arc. Animals 
were immobilized by a brief immersion in a freshly made 10 per cent solution of 
methyl alcohol in 0.1 N Ringer’s solution. The motionless animals were slightly 
flattened under a quartz coverslip for the irradiation. To observe and photograph 
animals no longer immobile, a microcompression chamber (Biological Institute, 
Philadelphia, Pennsylvania) was used to hold the animals still. The effects of the 
irradiation could be examined under phase-contrast microscopy seconds after 
irradiation simply by removing the irradiation nosepiece and substituting for it a 
nosepiece carrying phase objectives. The versatile Heine condenser was used for 
both the phase and the UV objectives. The general methods of culturing the ani- 
mals are those given by Sonneborn.® 

Materials.—The clone of double enimals used here arose spontaneously from the 
homopolar fusion of two conjugants of stock 51, variety 4, P. aurelia. Conspicuous 
morphological features of double animals are (1) two sets of oral structures, one set 
on each side of the body; (2) four contractile vacuoles, two on each side midway be- 
tween the feeding organelles; and (3) large size. 

Observations During and Shortly after Irradiation.—The ability of the microspot to 
cause localized damage is best illustrated by irradiation of a dividing animal. Of 
importance here is the close proximity of the dedifferentiated old gullet to the anlage 
of the new gullet which it has just produced (see Fig. 1, a). The anlage at this early 
stage shows convulsive movements. It was irradiated until all movement ceased, 
about a minute, and subsequently for another two minutes. Activity of the near-by 
old gullet was not visibly impaired during the course of the irradiation. Two hours 
later this animal had completed dividing, and observation revealed two normal gul- 
lets in the proter (Fig. 1, b) and only one in the opisthe (Fig. 1, c). The irradiation 
had eliminated the second gullet normally present in the opisthe. 

The gullet damage depends on (1) placement of the microspot, (2) the stage of de- 
velopment of the gullet, and (3) the length of time of irradiation. The average 
length of the gullet is about 20 uw in nondividing animals. This allows the small 
area of the microbeam to damage selectively parts of the gullet. At least three dif- 
ferent placements of the microspot on the long axis of the gullet are necessary to 
stop ciliary beating in the whole gullet. If only one placement is used, irradiation of 
the ventral side of the esophagus results in the most extensive damage for a given 
time of exposure. Newly formed gullet anlagen lose all activity about twice as fast 
as the matured gullets. Dedifferentiated gullets are not so sensitive to the irradiation 
as are the anlagen but are more sensitive than the mature gullets. Longer ex- 
posures, of course, result in greater irradiation effects. 

The gullet on the nonirradiated side is normal in appearance and function. On 
the irradiated side the damage can range from loss of a gullet, through a definite 
series of intermediates, to no detectable damage. The essence of these intermediate 
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stages is that a given part of the gullet is present only if all the more anterior parts 
are already in place. 

Other more generalized effects are also caused by the irradiation. Occasional 
erratic bursts of trichocyst discharge are often seen during the course of the irradia- 
tion. Also, gullet radiation can delay fission if administered prior to the onset of 
visible manifestations of division. As compared to controls, which divide within 5 
hours, the comparable irradiated animals can be delayed up to 24 hours. 

Restitution and Regression of Gullet Damage.—In confirmation of Tartar,? it has 
been found that no regeneration occurs if the damaged animals are starved. How- 
ever, there is the new finding that feeding animals can reform lost or damaged gullet 
structures. 

In two cases, animals with only one gullet, but two oral grooves, and four contrac- 
tile vacuoles gave rise to clones in which about one-fourth of the animals were nor- 
mal doublets. Also included in these clones were forms intermediate to singles and 
doubles, e.g., possessing one gullet and three contractile vacuoles. Reorganization 
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Fic. 1.—a:; Gullet region of dividing animal just prior to irradiation (Bausch and Lomb UV 
reflecting objective). 6: Proter from same animal two hours later (phase contrast). ¢: Opisthe 
from same animal two hours later (phase contrast). d = Dedifferentiated gullet; a = Anlage 
of new gullet; g = Mature gullet. 


into singles is the most common fate of doublets which have lost a gullet. Such re- 
organization has been reported for doublets of Leucophrys.’ 

Regression of a damaged gullet has also been observed. Animals with one normal 
and one damaged gullet have given rise to progeny with only one gullet. 

Inheritance of Gullet Damage.—In other cases, behavior of a gullet rudiment was 
followed fission by fission. Three major findings emerge from this study. First, if 
only a very small gullet rudiment, e.g., mouth and a small part of the pharynx, is 
present on the irradiated side, the opisthe will receive no gullet structures on that side 
but the proter will retain the rudiment. Apparently no anlage-forming region is 
present. Second, larger rudiments, e.g., mouth and pharynx, will pass on to the 
proter, as above, but also will provide the opisthe with a damaged rudiment. This 
indicates the presence of a damaged anlage region. Third, a damaged gullet, con- 
taining, perhaps, mouth, pharynx, and most of the esophagus, can give rise to a 
normal gullet in the opisthe while persisting as an injured structure in the proter. 
Here the anlage region is evidently undamaged. It has been observed that these 
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three categories represent the chronological order of events in recovery of the nor- 
mal doublet phenotype. After each succeeding fission, in such a case, less damage 
is present than in the preceding interfission animal. 

The data strongly suggest as an explanation for the regression and recovery of 
damaged gullets that the proter line of inheritance of an injured gullet will eventu- 
ally lose the damaged structure and that the recovery of the normal doublet condi- 
tion occurs only in an opisthe line of descent. 

Discussion.—The mode of gullet formation suggests that the production of new 
parts depends on pre-existing ones. This suggestion is in keeping with the expecta- 
tions of Lwoff’s* kinetosome theory, according to which all morphogenesis in ciliates 
is under the control of the kinetosomes and new kinetosomes arise from pre-existing 
ones. Thus regeneration of the gullet structures as reported here could conceivably 
depend on the replacement of destroyed kinetosomes by reproduction and dif- 
ferentiation of the remaining ones. However, the origin of the new kinetosomes in 
the case in which the whole gullet has been lost is enigmatic. Two alternative 
sources are possible: (1) the remaining normal gullet and (2) the pellicular kineto- 
somes. Resolution of this problem awaits further study. 

Finally, there should be mentioned an attempt to obtain damage in the opisthe 
and none in the proter by irradiating a nondividing animal. This was successful 
and was achieved by irradiating only the right wall of the pharynx, i.e., the anlage- 
forming region. The experiment was not an unqualified success. The irradiated 
animal was large, and though no formation of a new gullet was visible, it is possible 
that formation of a new gullet had already commenced and only the newly formed 
structures were irradiated. At the next division the normal-appearing proter gave 
rise to two normal animals, showing either that a damaged anlage region did not 
persist or, as already suggested, that the irradiation was delivered to a very early 
stage of formation of the opisthe’s gullet. 

Summary.—A microbeam of ultraviolet light has been used to damage one gullet 
of the two which exist in homopolar doublets of P. aurelia. Irradiation can result 
in varying amounts of damage to this structure, even causing complete loss of the 
irradiated gullet. Both regeneration and regression of the damaged structure have 
been observed. Double animals entirely lacking one gullet commonly reorganize 
as singles but can also regenerate the entire lost gullet. An injured gullet can pro- 
duce either no new gullet, a damaged gullet, or a normal gullet. These differences 
in behavior are explained by assuming different degrees of damage to the gullet 
anlage-forming region. 

* The author wishes to thank Drs. R. E. Zirkle, W. Bloom, and R. B. Uretz for the use of their 
laboratory, where this work was initiated. Special thanks are due Dr. R. Uretz for his patient 
counseling during the construction of the author’s microbeam apparatus. 
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DIFFERENTIATION IN CULTURE OF MIXED AGGREGATES OF 
DISSOCIATED TISSUE CELLS* 


By J. P. TRINKAUS AND Peccy W. GROVES 
OSBORN ZOOLOGICAL LABORATORY, YALE UNIVERSITY 


Communicated by J. S. Nicholas, July 1, 1955 


Our concepts of the process of determination in embryonic development have 
been based largely on results obtained by extirpation, transplantation, and explanta- 
tion of tissue masses consisting of many cells. These methods have revealed an 
apparent restriction of potencies during development. The varying results ob- 
tained with different methods, however, have led to serious questioning of the con- 
clusions. For, as Harrison! pointed out, “‘there always may be new conditions, 
not yet tested, under which other potencies might be revealed.” 

An important factor influencing the degree and variety of differentiation in 
tissue culture is the size of the tissue mass involved.?~> Thus small pieces from the 
node region of the chick blastoderm (definitive primitive streak stage) fail to dif- 
ferentiate in culture, while larger pieces invariably form brain tissue. Explants 
of presumptive head mesoderm from a single early gastrula of Triton taeneatus 
form only muscle in vitro, but larger masses consisting of identical cells, derived 
from the fusion of several of these same primordia, ofrm muscle, chorda brain 
sensory primordia, and ectoderm.’ It is evident from these observations that a 
description of the potencies of a given tissue mass composed of many cells does not 
necessarily reveal the capacities of the individual constituent cells. 

These studies raise the possibility that, when an intact tissue mass is judged to 
be determined by all the usual criteria, individual cells within the mass may retain 
wider potencies. Perhaps in transferring a mass of cells to a new site the immediate 
environment of the inner cells is but slightly altered. Such a change would be 
sufficient to cause a shift in the direction of cell differentiation in the earliest phases 
of development (e.g., early amphibian gastrula); but in more advanced stages, 
generally considered to be determined, a more acute change may be necessary. 
Accordingly, in the latter instance, fixity of the mass under the imposed conditions 
could not be taken as evidence of cell fixity. 

With the discovery that it is possible to dissociate vertebrate tissue cells by treat- 
ment with high pH,® trypsin,” * or versene® and, after aggregation, to obtain char- 
acteristic differentiation in vitro, we have a promising method for attacking this 
problem. By mixing dissociated cells derived from different primordia (or even 
from differentiated tissues) and combining them into complex tissue masses by 
aggregation, we may drastically alter the immediate environment of embryonic 
tissue cells. The differentiation of such complex aggregates may then provide 
important information on the potencies of these cells. 

It was with this hope that we began the present investigation. Our method has 
been to combine cells derived from two different organs at various advanced stages 
of organogenesis and study the differentiation tn vitro of the complex aggregates. 

METHODS 

The organs used in this study were mesonephros and wing bud from 3'/:- and 

5-day chick embryos (stage 19-20 and stage 27)."° White Wyandotte, Rhode Island 
787 
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Red, and White Leghorn eggs were used, with no differences in result. To obtain 
an adequate mass of aggregated cells, large numbers of these organs were required 
(e.g., 30-40 mesonephroi and wing buds from stage 19 embryos). 

Both trypsin and versene (disodium versenate) were used as dissociating agents 
during beginning phases of the investigation. Trypsin was soon found to be more 
satisfactory for our purposes, however, because after dissociation with this agent 
the cells tend to readhere to one another more readily. This yields a firmer ag- 
gregated mass. A 3 per cent trypsin solution was used (Difco 1:250 trypsin in 
Ca- and Mg-free Tyrode’s solution, buffered with phosphate and bicarbonate). 
The trypsin solution, Tyrode’s solution, and bicarbonate buffer were all sterilized 
by filtration through a fritted glass disk. This not only provides sterile media 
but also avoids the occasional cotton fibers that fall from plugs used in autoclaving, 
and which adhere firmly to dissociated tissues. Standard sterilization procedures 
were used for all other solutions and glassware. All solutions were prepared with 
demineralized water obtained by passing metal distilled water over an anionic- 
cationic exchange resin (Amberlite MB-1). In our experience this is quite as 
satisfactory for tissue cultures as Pyrex-distilled water. 

Suspensions of dissociated cells of both mesonephros and wing bud were obtained 
by a modification of the procedures devised by Moscona.’?. Entire wing buds and 
mesonephroi, cleaned of most of the extraneous tissue, were placed in trypsin I (pH 
raised to 7.8 with 1 per cent KOH). The tissues were cut into small pieces in the 
enzyme solution and then incubated at 38° C. for about 15 minutes. Upon re- 
moval from the incubator, the tissues were drawn into an orally controlled pipette 
of small bore (0.08—0.1 mm. inner diameter). This causes the tissue to fragment 
into a mixture of free cells and cell clusters, composed of ca. 20-100 cells. The 
partly dissociated tissue was then forced out of the pipette into trypsin IT (the same 
as trypsin I, except that the pH has been raised to 8.3 with the addition of 1 per cent 
KOH) and incubated at 38° C. for an additional 15 minutes. The suspension was 
then once again drawn into a very fine bore pipette (ca. 0.05 mm. inner diameter) 
and forced out of the pipette into a dish of Tyrode’s solution. This procedure yields 
a suspension consisting almost entirely of individual cells (Fig. 1). However, a few 
pieces of ectoderm, consisting of 5-10 cells, always remain undissociated in suspen- 
sions derived wholly or partly from 5-day wing buds. The entire process of dis- 
sociation, from the time the tissue is placed in trypsin I until the dissociated cells are 
suspended in Tyrode’s solution, requires approximately 1!/, hours. 

The cells were washed free of trypsin and aggregated by the following procedure: 
centrifugation (2-3 minutes in a clinical centrifuge at 1610 & g), removal of the 
supernatant, replacement with fresh Tyrode’s and centrifugation for an additional 
few minutes, removal of the supernatant, and replacement with a nutrient medium 
consisting of 2 parts horse serum, 2 parts Tyrode’s solution, and 1 part embryo ex- 
tract. The cells were centrifuged for 7-10 minutes in this medium, in order to 
force them into a more compact mass. 

This aggregated mass of cells was removed intact from the centrifuge tube with 
the aid of a small spatula and cut into pieces of the desired size. One of these was 
fixed (Fig. 2), and the others were cultured overnight in the fluid nutrient medium in 
deep depression slides. The following day the explants were removed from this 
medium and cultured by the Fell method on the surface of a plasma-embryo extract 
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clot. During the sojourn in fluid medium the explants acquire a firm enough con- 
sistency to make culturing on the surface of the clot a feasible procedure. While 
this method of aggregation has been satisfactory in providing us with firm aggre- 
gates which differentiate in culture, it has one distinct drawback. There is con- 
siderable loss of peripheral cells, both in cutting the aggregated mass into explants 
and during the several hours in nutrient medium. Because of the possibility that 
differences in specific gravity may cause cells to separate during centrifugation, 
care was taken to cut the aggregated mass in such a manner as to insure inclusion 
of cells of varying specific gravity in each explant. The subsequent complex dif- 
ferentiation characteristic of larger explants suggests that representatives of all 
cell types are present in each explant (Figs. 6, 7, 9) 

The clot was formed with 4 parts hen plasma, | part fresh embryo extract, and 
0.5 part 0.001 per cent protamine sulfate (to achieve a firm clot). Embryo extract 
was made from embryos which, as far as possible, are the same age as the cells of the 
explant.!! For example, aggregates from a 5-day embryo after one night in fluid 
medium were placed on a clot made with 7-8-day embryo extract, were transferred 
1 days later to a clot made with 10-day extract (the same age as the tissue cells), and 
{ days later to a clot made with 14-day extract. Explants were routinely shifted 
to a fresh place on the clot after 2 days and were transferred to a fresh clot every 4 
days. In order to prevent bacterial growth, penicillin (Potassium Penicillin-G, 
Lilly; 125 units/ec) and streptomycin (CaCl complex, Merck; 0.875 mg/cc) were 
regularly added to the Tyrode’s solution. These antibiotics have no detectable 
effect on the explants, as determined by comparison with parallel control cultures. 
Explants cultured by this method remain remarkably healthy for long periods in 
cultures (up to 16 days). Cytolysis, as judged by pyenosis, is at a minimum, being 
largely confined to the central area of the explant and to the first day or so in cul- 
ture. Mitoses occur throughout the period in vitro but are most frequent during 
the first few days, prior to the differentiation of tissues. Explants were fixed in 
Bouin’s solution after varying periods on the clot, from 1 to 12 days. Sections were 
cut at 6 w and stained with Delafield’s haematoxylin or Mallory’s triple stain. 

When mixed aggregates were desired, known quantities (by wet-weight determina- 
tion) of mesonephros and wing bud tissue were disaggregated together in the same 
dish. Thus mixing begins in trypsin I, continues in trypsin II, and is completed 
in Tyrode’s solution, where the cell suspension is drawn into a medium-bore pipette 
and flushed out several times. All other procedures were as described above. 

RESULTS 

Aggregates of Mesonephros Cells——All cell aggregates of dissociated 3'/.-day 
(stage 19-20) and 5-day (stage 27) mesonephros (a total of fifteen cases) formed 
characteristic mesonephric tubules after 3-6 days in watch-glass culture (Fig. 3). 
This confirms the results of Moscona and Moscona.°® 

Aggregates of Wing Bud Cells.—Aggregates of dissociated wing buds of embryos at 
stage 19-20 and stage 27 ( a total of ten cases) do not differentiate during the first 
few days in culture. By the fifth day on the clot, large areas of procartilage dif- 
ferentiation appear (Fig. 4). In no case, however, did cartilage development pro- 
ceed beyond this stage, even though explants were cultured for as long as 13 days. 
Furthermore, no other wing tissues appear (e.g., mesenchyme, muscle, or kerati- 
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nized epidermis), and all explants possessed large areas of undifferentiated cells. 
These results are at variance with those of Moscona and Moscona,’ who obtained 
excellent cartilage and wing mesenchyme formation in aggregates of limb mesoderm 
cells. In view of the results presented in the next paragraphs, the method of cul- 
turing cannot be held responsible. The only differences that appear significant 
are a longer exposure to trypsin and the presence of ectoderm cells in our explants. 

Mixed Aggregates of Mesonephros and Wing Bud Cells from Embryos of the Same 
Stage.—Aggregates composed of intermingled cells derived from dissociated mes- 
onephroi and wing buds invariably form tissues when cultured by the watch-glass 
method; however, the degree of differentiation and the number of tissue types 
formed vary considerably. 

The data on which the following observations are based are derived from study 
of twelve explants composed of dissociated cells from approximately equal quan 
tities of mesonephros and wing bud tissue of 3!/2-day embryos (stage 19-20). 
In one series 2.8 mg. of mesonephros and 2.0 mg. of wing buds yielded six explants. 
In another series equal quantities of mesonephros and wing bud tissue were used 
(5.5 mg. of each). This mass also yielded six explants. Comparable results were 
obtained in the two series. After 4 days in culture the beginnings of characteristic 
cartilage and mesonephric tissue are evident Areas of procartilage and of epi- 
thelial tissue appear (Fig. 5). Explants cultured for 7-9 days or more possess 
well-differentiated cartilage and mesonephric tubules (Fig. 6) and often sheets or 
whorls of keratinized material. In spite of the fact that the dissociated cells were 
thoroughly intermingled prior to aggregation, the individual tissues almost. in- 
variably occupy separate sectors of the explant. Moreover, the tissues generally 
have a characteristic regional relation to one another: cartilage in the central 
region of the explant, mesonephric tubules in a perichondrial position, and keratinized 
sheets and whorls near the periphery. No loose wing mesenchyme formed in 

Abbbreviations: c¢ = cartilage;/m = loose mesenchyme; mt = mesonephric tubule; pe = pro- 
cartilage; py = pycnotic nuclei. All photomicrographs except Figure 1 are of sectioned material. 

Fig. 1.—Dissociated cells of 5-day chick wing bud. Fixed directly on a slide immediately after 
dissociation. Stained with Delafield’s haematoxylin. 1,500. 

Fig. 2.—Mixed aggregate of dissociated cells of 3!/2-day mesonephroi and 3'/2-day wing buds. 
Fixed immediately after centrifugation. Stained with Delafield’s haematoxylin. 160. 

Fic. 3.—Explant derived from an aggregate of dissociated cells from 3'/:-day mesonephroi. 
Note well-formed mesonephric tubules. Seven days in culture. Stained with Delafield’s haema- 
toxylin. 160. 

Fig. 4.—Explant derived from an aggregate of dissociated cells from 5-day wing buds. Note 
absence of characteristic cartilage and loose mesenchyme. Only procartilage areas are present. 
Thirteen days in culture. Stained with Delafield’s haematoxylin. 160. 

Fig. 5.—Explant derived from a mixed aggregate of equal quantities of 3'/2-day mesonephros 
and 3'/.-day wing bud cells. Note the areas of beginning cartilage and mesonephric tubule for- 
mation. Four days in culture. Stained with Delafield’s haematoxylin. 160. 

Fig. 6.—Explant derived from a mixed aggregate of equal quantities of 3'/2-day mesonephros 
and 3'/:-day wing bud cells. Note regional differentiation of cartilage and mesonephric tubules. 
Eight days in culture. Stained with Delafield’s haematoxylin. 160. 

Fig. 7.—Explant derived from a mixed aggregate of equal quantities of 5-day mesonephros and 
5-day wing bud cells. Note regional differentiation of cartilage and mesonephric tubules. Ten 
daysinculture. Stained with Mallory’s triple stain. 160. 

Fig. 8.—Explant derived from a mixed aggregate of dissociated cells from 5 mg. of 3'/.-day 
mesonephroi and 8.6 mg. of 5-day wing buds. Note predominant differentiation of mesonephric 
tubules and absence of other tissues. Ten days in culture. Stained with Delafield’s haema- 


toxylin. 160. 

Fic. 9.—Explant derived from a mixed aggregate of dissociated cells from 8.3 mg. of 5-day 
mesonephroi and 3.6 mg. of 31/2-day wing buds. Note extensive formation of cartilage and loose 
mesenchyme, in contrast to the small number of mesonephric tubules. Ten days in culture. 
Stained with Delafield’s haematoxylin. 160. 
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explants of this series. Undifferentiated cells are always present but usually 
comprise but a small part of the explant. 

The differentiation of a similar series of four explants derived from an aggregate 
of mesonephros and wing bud cells of 5-day embryos (stage 27) also was studied. 
Cells from 21.6 mg. of mesonephros were mixed with cells from 22.3 mg. of wing 
bud tissue. These mixed aggregates differentiate in a manner similar to that 
just described for the 3'/.-day—3! .-day combination (Fig. 7). In these cases, 
however, advanced differentiation of mesonephric tubules and cartilage occurs 
sooner-—by the fifth day in culture. 

It is clear from these results that characteristic mesonephric tubules, cartilage, 
and stratified epithelium may differentiate with consistent regional arrangement 
in mixed aggregates of dissociated cells from mesonephros and wing bud, when 
these primordia are derived from embryos of the same developmental stage. It 
is now of interest to consider what happens when cells from an early stage of one 
primordium are mixed with cells from a later stage of the other. For, if the deter- 
mination of the constituent cells and the mass of tissue are equal at any given time, 
then cells from a more advanced stage would have less capacity to regulate. Accord- 
ingly, if any tissues predominate in an aggregate composed of cells from embryos of 
different ages, they should be those characteristic of the more advanced organ 
(i.e., in a mixture of cells from early wing buds and later mesonephroi we should 
expect mesonephric tubules to predominate). 

Mixed Aggregates of 3'/s-Day Mesonephros and 5-Day Wing Bud Cells.—An 
experiment using approximately equal quantities of 3'/s-day mesonephros and 5- 
day wing bud tissue (7.0 mg. each) yielded three mixed aggregates. After 12 days 
in watch-glass culture, all these explants possessed well-differentiated mesonephric 
tubules, a small amount of epidermal tissue (whorls of stratified epithelium), and 
small areas of undifferentiated cells. None formed cartilage or loose mesenchyme. 
Furthermore, in each instance the mesonephric tissue clearly comprises the great 
bulk of the explant! The predominant presence of mesonephric tissue in these 
explants was of such interest that we felt a more rigorous test to be desirable. 
Thus, in repeating the experiment, we combined a smaller quantity of 3'/s-day 
mesonephros cells (5 mg.) with a larger quantity of 5-day wing bud cells (8.6 mg.). 
This aggregate yielded four explants, which were cultured 7-9 days. In all of these, 
mesonephric tubules formed half or more of the bulk of the explant, the rest being 
largely composed of undifferentiated cells (Fig. 8). In two explants, however, 
there is evidence of beginning cartilage differentiation. 

Mixed Aggregates of 3'/.-Day Wing Bud and 5-Day Mesonephros Cells—When 
cells from a small quantity of 3'/s-day wing bud tissue (3.6 mg.) were mixed with 
cells from a much larger quantity of 5-day mesonephros (8.3 mg.), we obtained 
five explants which became differentiated during 7-9 days in culture. Wing bud 
tissues predominated in all these cases. The two largest explants formed ex- 
tensive areas of well-differentiated cartilage and loose mesenchyme and a small 
amount of poorly differentiated mesonephric tissue. Regional differentiation is 
characteristic—cartilage in the central areas, mesonephric tissue largely peri- 
chondrial, and mesenchyme most cortical (Fig. 9). The other three explants of 
this series possess both cartilage and mesonephric tubules, the former being 
definitely predominant. 
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From these results it is abundantly evident that the younger tissues predominate 
in mixed aggregates composed of mesonephros and wing bud cells from older and 
younger embryos. Not only do they comprise the greater portion of the differen- 
tiated explant, but they also differentiate to a more advanced degree. Moreover, 
these relationships persist even when a relatively smaller quantity of the 3'/2.-day 
tissue is used in forming the aggregate. 


DISCUSSION 

The remarkable differentiation of these mixed aggregates raises a number of 
interesting questions (e.g., why do cartilage and loose mesenchyme differentiate so 
much better in the presence of mesonephros cells?); but the central question 
concerns the cellular origin of the differentiated tissues. Basically, there are two 
possible answers to this question. Either each tissue is wholly derived from cells 
of the same tissue type (mesonephric tubules from mesonephric cells, cartilage from 
prospective cartilage cells, etc.), or each tissue is in part derived from cells of 
diverse origin, by a process of cell transformation (e.g., cartilage from mesonephros 
cells). We shall refer to these possibilities as the redifferentiation and transforma- 
tion hypotheses, respectively. In both these hypotheses we must take into account 
the import of critical mass in tissue differentiation. Several lines of evidence 
indicate that with standard culture conditions a mass of cells must exceed a certain 
minimal size before differentiation may occur in vitro.*4 It therefore seems 
justifiable for us to assume that differentiation of any particular tissue in our 
mixed aggregates will not occur unless the required number of cells of that tissue 
type are clustered. 

Redifferentiation Hypothesis—There appear to be three possible mechanisms 
whereby the regionally differentiated tissues of our explants may originate by 
redifferentiation : 

1. Redifferentiation in situ of the cells of each tissue type. This mechanism 
assumes there to be an occasional chance association of several cells of one tissue 
type, which increases in size by mitosis, if it is located in the proper region of the 
explant. Once the necessary size is reached, redifferentiation occurs. The cells 
of different tissue origin which are located in the midst of this differentiating mass 
will remain undifferentiated. In the case of an unmixed tissue, like cartilage, 
they must be eliminated. This may occur as a result of cytolysis or migration to 
another part of the explant or into the zone of outgrowth (negative cell movement). 

2. Redifferentiation of the cells of each tissue type after they have been sorted 
out by directed cell movements. Thus procartilage cells would generally move 
toward the center of the explant, mesonephric cells to the perichondrial region, 
and wing mesenchyme and ectoderm cells to the periphery. 

3. Redifferentiation of the cells of each tissue type after they become sorted 
out by random cell movements. This mechanism assumes that each tissue cell 
has a specific affinity for others of its type, so that upon contact they adhere.!® 
In this manner, masses of cells of each type will form. Redifferentiation will occur 
in those that meet the requirements of mass and regional differentiation (cartilage 
in the center, mesonephros in the perichondrial region, ete.). The remaining cells 
either will remain undifferentiated or will be eliminated by cytolysis or negative 
cell movement. 
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Cell Transformation Hypothesis—This hypothesis also assumes that differentia- 
tion of a particular tissue begins by redifferentiation (differentiation in accord with 
prospective fate) of a small mass of cells derived from that tissue. The manner 
of formation of such a cell cluster may be by one or a combination of the mechanisms 
just discussed. When this cell mass reaches the necessary size, it will begin to 
differentiate, if it is located in the proper region of the explant. It is assumed 
that with the onset of differentiation this mass becomes an active center, specifically 
influencing surrounding cells to follow suit. As these surrounding cells begin to 
differentiate, they, in turn, will influence cells peripheral to them. This process 
continues until a large area composed of many cells is involved. It ceases when 
cells are reached which have already passed under the influence of another active 
center. This hypothetical phenomenon resembles an embryonic field. If most 
of the cells of the explant remain randomly distributed during this process, it follows 
that some cells derived from one tissue must become transformed into cells of 
another tissue. This hypothesis requires neither differential mitosis, differential 
cytolysis, nor negative cell movement; nor are extensive directional cell movements 
required. 

In view of the comprehensive studies of Holtfreter® and Townes and Holtfreter'® 
on mixed aggregates of amphibian gastrula and neurula cells, redifferentiation 
from cells of the same tissue type after directed cell movements emerges as a likely 
possibility. It must be emphasized, however, that even though there is a general 
tendency in the amphibian material for prospective epidermis cells to move out- 
ward and for mesoderm and neural plate cells to move inward, the specific fate of 
the inner mesodern cells of such aggregates is not known. Thus the possibility of 
cell transformation has not been eliminated. 

A number of considerations suggest that cell transformations may occur in the 
differentiation of mixed aggregates of chick cells: 

1. The striking differentiation of complex aggregates composed of mesonephros 
and wing bud cells from embryos of different ages into tissues that are almost 
entirely of one organ type (e.g., mesonephric tubules [Fig. 8]) is more readily ex- 
plained by assuming a transformation of the cells of the older tissue. In addition, 
the predominance of tissues from the younger organ suggests that properties such 
as greater migratory activity and adhesiveness may be more significant in the dif- 
ferentiation of the aggregates than the presumed determination of the cells from 
the more advanced organ. 

2. Prospective wing mesenchyme cells are present in abundance in all aggre- 
gates of wing bud and mesonephros cells, yet loose mesenchyme differentiates only 
in large explants (Fig. 9). The absence of this tissue in small explants is most 
directly explained by assuming that some prospective mesenchyme cells have been 
transformed. 


3. Our evidence does not favor differential cytolysis. Pyenotie nuclei (our 
criterion of cytolysis!) are frequent during the first 2 days in culture, but their 
central location suggests that the cytolysis is due to general factors such as nutri- 
tional or respiratory deficiency and is not related to the regional differentiation of 
tissues. In general, explants are virtually free of cytolyzing cells after the first 


3 days on the clot. 
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CONCLUSIONS 
While it is clear from these observations that we must seriously consider the 
possibility that cell transformations are occurring during the differentiation of 
these mixed aggregates, the evidence is, in fact, equivocal at this time. We have 
as yet no information on cell movements in the differentiating aggregates and no 
direct information on cell origin. Evidence on this critical matter can be gained 
by marking one type of cells and following their fate in the explant. An effort 
to gain evidence of this sort, using S® as a marker, is now in progress. A decision 
as to which of the proposed hypotheses, or other possible views, is correct must 
be deferred until the completion of these and other marking experiments. 


* Aided in part by a grant from the National Science Foundation (G1024). Disodium versenate 
was generously supplied by Versenes, Inc. We are indebted to Mr. Thomas D. Sauerwein for 
technical assistance. 
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ERRATA: DIGITAL COMPUTATIONAL METHODS IN SYMBOLIC 
LOGIC, WITH EXAMPLES IN BIOCHEMISTRY 


In the article of the above title appearing in these PRocEEDINGs, 41, 498-511, 
1955, the following corrections should be made: 


On page 502 diagram at lower left should appear as follows: 


ANTECEDENCE 
(F jx) ® (E..:) cS (R;:) 


. ; 110 
On page 504, second sentence of second paragraph, instead of “#(R;;) -( Bs 


101 
110\ ,, 
read “(Rj;) = borg, 


On page 511, note 9, read “...G. Gamow, A. Rich, and M. Yéas, article in Ad- 
vances of Biological and Medical Physics... .” 


On page 501 the automorphism H might more customarily be called a homo- 
morphism of B into itself. 
R. S. LEpLEY 
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